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Preliminary Designing Method of External 
Pressure Vessels for Sea Subsurface Applications 


Jan P. Michalski, Assoc. Prof. 
Gdansk University of Technology 
Polish Naval University 


ABSTRACT 


This paper describes a unique engineering method intended for the preliminary designing 
of marine echo-location systems. The solved designing problem consists in determining 
geometrical parameters and selecting structural materials for a vessel of a buoyancy 
necessary to contain measuring instruments of a given mass, as well as its maximum 
gabarites and operational submersion depth. The assumed variability range of the system 5 
parameters determines a space of permissible applicability parameters of the method. 
Stress level related to permissible and critical stresses was assumed to be safety criterion 


for the vessel. The method may be applied also to solving converse problems consisting in 

determining maximum permissible submersion depth for a vessel made of assumed structural materials 

and having given geometrical parameters. The presented description of the method is illustrated by an 
example of its application. 


Keywords: subsurface buoy mooring systems; external pressure vessels 


INTRODUCTION 


Tight submerssible vessels intended for the fulfilling 
of assumed functions deep in water, capable of carrying 
large external hydrostatic loads find many civil and military 
applications such as submarines and submersible vehicles, 
bathyscaphes, crude-oil and gas storage tanks, or casings for 
measurement apparatus to conduct underwater oceanological 
experiments [ 3]. For example, the measuring, recording and 
investigating of physical fields of seas and oceans is conducted 
by means of echo-location systems (Subsurface buoy mooring 
systems) composed of such elements as: 

— floats (buoys of a required buoyancy) 

— vessels containing measuring and recording apparatus 
— vessels containing electric storage batteries 

— anchoring device in the form of ballast 

— tension members joining elements of a module. 


The system’s elements are connected to each other by 
means of elastic tension members, i.e. cables or chains. 
The vessel’s shell plating is loaded by hydrostatic pressure 
proportional to submersion depth, sometimes also by 
a hydrodynamic load due to water current or surface wave 
motion. Also, concentrated loads generated by tension 
members are applied to the vessel. For functional reasons 
the vessel’s structure should ensure its buoyancy, shape and 
watertightness to be maintained. 

This paper contains a description of the engineering 
method for the preliminary designing of echo-location system, 
which serves to determining design solutions satisfying a set 


of functional, technological and mechanical assumptions. 
Original algorithms for the determining of parameters of vessels 
subjected to action of external pressure — in the aspect of their 
buoyancy, shell structure strength and form stability — are also 
presented in it. 

Applicability of the method is determined by ranges of 
permissible values of vessel’s form parameters, external 
pressure values as well as a kind of structural material. 
The method may be useful in particular for the preliminary 
designing of echo-location systems (Subsurface buoy mooring 
systems) within the area of: 

— designing their geometrical configuration 

— designing geometrical parameters of their vessels 

— selection of their structural materials 

— assessment of quantity of their mass 

— determining optimum design solutions in an assumed 
sense. 


FORMULATION OF DESIGNING 
PROBLEM 


The designing problem consists in determining geometrical 
parameters and selection of a structural material for a vessel 
of a given effective buoyancy, capable of accommodating 
measurement apparatus of a given mass and dimensions, at 
a given operational submersion depth. Variability range of 
parameters for which calculation algorithms of the method have 
been elaborated, determines a space of permissible parameters 
for its application. The method makes it possible to solve also 
converse problems, namely determination of a maximum 
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permissible submersion depth of a vessel made of a given 
structural material and having given geometrical parameters. 
An axially symmetrical vessel is considered of empty 
interior (a float) or containing measurement apparatus of the 
mass m, and maximum linear dimension 1 ; a part of the vessel 
can be filled with a lightweight material, e.g. polyurethane 
foam in order to maintain buoyancy of the vessel in the case 
of flooding its interior by water due to damage or loss of 
watertightness of its shell plating. The effective displacement 
D, of the vessel of the volume V, plating surface area S and 
constant shell plating thickness h, immersed in water of the 
constant density p,, is described by the buoyancy equation in 
which: 
— o=0- in the case of keeping the vessel tight 
— a= l -in the case of flooding the vessel entirely 


m 
D, =(l—a.)-py:V+a-py: ll 
Pa Pf 
F (1) 
—(m,+m, +mr)+—4 

8 
where 
g - gravity acceleration 
m,- mass of vessel plating 
p, - density of vessel plating material 
m, mass of apparatus 
p, - average density of apparatus 
m,- mass of filler 
P, - density of filler (foam) 
F - value of external concentrated load 


Balance of components of mass and displacement of vessel 
in Eq. (1) can be described by the following important design 
parameters of the system: 

— the geometrical parameters of the vessel, V, S, h 
— the density of structural material of its shell plating, p, 
— the mass of apparatus inside the float, m,, and its average 

density p, 

— the maximum linear dimension of the apparatus, |, 
— the filler volume V,and its density p, 
— value of the external load F 


The balance equation is then expressed by the equivalent 
relation as follows: 


De = (1-2): 0 Vta pp tsn} 
(2) 


a 
E, 
-S-h-p,—m, Vf: ppt 


Positive value of the right-hand side of Eq.(2) as well as 
positive sign of the external load Fz stands for senses of forces 
opposite to gravity force. In Eq. (2) the belonging of parameters 
either to the set of the parameters to be determined or that 
of given ones constitutes a tag of variant of the method for 
preliminary designing the float. To achieve an effective solution 
of Eq. (2) it is necessary to add complementary equations 
whose set and genealogy constitute another tag of the designing 
method. In the case of the method formulated in the form of an 
optimization problem Eq. (2) should be complemented by: 

— acriterion of choice of an optimum solution 
— aset of inequalities which restrain allowable solutions. 


The elaborated calculation algorithms of the method 


concern the designing process of vessels having axially 
symmetrical shells of circular cross-section, made of: 
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— NW 210 steel of normal strength 
— AL 5085 aluminium alloy usually applied in shipbuilding. 


The range of permissible parameters of the method in 
question was so selected as to make it possible to obtain 
geometrical dimensions, form proportions, structural material 
as well as maximum operational submersion depth which would 
satisfy requirements for echo-location research systems for 
Baltic Sea environment. 

A hazard to structural safety of a vessel subjected to action 
of high external pressure may result either from exceedance 
of reduced permissible stresses or critical stresses leading to 
loss of form stability (buckling) of its shell plating [1, 2, 4], 
in consequence to loss of its functional qualities. The method 
in question covers the whole range of elastic deformations 
up to the material yield point R ; the vessel’s parameters are 
so determined as to obtain reduced stresses as well as critical 
stresses which do not exceed their permissible values at an 
external load enlarged by safety factor. 


ASSUMPTIONS AND GEOMETRICAL 
DESCRIPTION OF FORM OF VESSEL 


Vessel’s shell plating is composed of a cylindrical casing 
of circular cross-section closed by two bottoms of a specific 
form which ensures moment-less state of stresses in the 
shell plating. The choice of such form of the vessel results 
from its usefulness to accommodate measurement apparatus 
and a simpler process of its manufacturing in contrast to 
spherical vessels. In accordance with the linear theory of 
stability of thin-walled elastic shells permissible values of 
the vessel plating thickness h were assumed much smaller 
than the curvature radius R of shell casing. Effects of possible 
geometrical imperfections of shell form were taken into 
account by safety factor increasing calculation external loads. 
The shell surface area of the vessel is assumed the same 
as the middle surface area of the shell plating. Structural 
mathematical models for determining structural strength 
and stability of shells as well as reccommended values of 
safety factors were taken from the respective subject-matter 
literature sources [1, 2, 3, 4]. 


Form of bottom’s surface 


The bottoms constitute support for vessel casing loaded 
by constant pressure. The desired momentless state of internal 
forces in shell plating will occur at an appropriate form of 
the bottoms. In the case of convex bottoms such form can 
be described by an appropriate plane curve (directrix) whose 
rotation around cylinder’s axis generates, at vessel’s bottom/ 
shell casing contact, a surface of the circumferential curvature 
radius R, = R and axial curvature radius R, = oo, which are 
consistent with curvatures of the shell casing. The requirement 
is fulfilled by the bottom surface of the directrix described by 
the following relation [4]: 


1 
t 


cai (1) 0 


R 
where 
tER>2.0 - index exponent - a real number 
b= - dimensionless radius of vessel, 0 < & < 1 
a - dimensionless camber of bottom, 1 < & <0 
1 - camber of bottom. 


Area of bottom’s surface 


Bottom’s surface area constitutes a part of vessel’s shell area 
and affects quantity of shell plating mass. The bottom’s surface 
area S, can be expressed by using Guldin theorem: 


(4) 


R 
and after substituting dimensionless variables: 


dz Y 
Sy=2-n- fr- 1+) — | -dr 
dr 
0 
1 1 2 dí 2 
S,=2-7-R?.- lE. (a) [5] a- 
a52: R E ta E 


1 (5) 
=2-n-R? > [sGsRula.t)-aé 


0 
Inverse function of the integral is not known hence it cannot 
be expressed by an analytical formula. In order to elaborate 
a simple approximating formula useful in formulating Eq. (2), 
parametric investigations of numerical values of the integral 
were conducted within the following ranges of the parameters: 
0 <1/R <0.5 and 2.5 < t < 4.5, by using the quadrature [6]: 


la ia AR ` 
sR tt ]=22 Seot Seat 2D S an 6) 


The best obtained analytical approximation of discrete 
values of the quadrature is represented by the expression [6]: 


S Rt zn- R? | le Ei t3 
d R 1 R = 


1 1.411 
=n? f140.4281 1.) ae 


Within the considered variability range of the investigated 
parameters the relative percentage error of values of the integral 
determined by means of Eq. (7) and the quadrature (6), does 
not exceed 0.3%. 


Bottom’s volume 


Volume of the bottoms which constitutes a part of the vessel 
volume, affects value of its displacement. The bottom volume 
can be determined by using the Guldin theorem: 


1 
R 


2\t 
r 
Va = 2n-l, |r fE) dr 
R 
0 


In the case of the considered form of the bottom the integral 
is non-infinitesimal, as its primitive function is not known. 
By substituting the dimensionless variables r = R - & and 
dr = R - dé the following is achieved: 


(8) 


Vz =2n +1, -R’. bay dé = 
a (9) 
=2-n--R- [ve t) -dé 
R > 


0 


In order to elaborate a simple approximating formula useful 
in formulating Eq. (2), parametric investigations of numerical 
values of the integral were conducted within the following 
ranges of the parameters: 0 < 1/R < 0.5 and 2.5 < t < 4.5, by 
using the quadrature [6]: 


l AR ` 
d = 
Va RE a a Neo Mert 2D orale 


On the basis of the obtained numerical values,by testing 
various hypotheses on a form of approximating functions, the 
following analytical formula for approximate determination of 
values of bottom’s volume was elaborated [6]: 


l Med 
Va RA] mt 2.Ró= 


ly ,02339 p3 a 
= 0.5747-1 rae -R 


Within the considered variability range of the investigated 
parameters the relative percentage error of values of the integral 
determined by means of Eq. (11) and the quadrature (10), does 
not exceed 0.3%. 

The complementary part of the whole volume of the vessel, 
V, is the following volume of its cylindrical part, V „ of the 
length L: 

V_=2aRL*V=2V,+V, (12) 
CRITERION OF STRUCTURAL STRENGTH 
SAFETY 


As a measure of material effort of structure under external 


load the reduced stresses are assumed usually in compliance 
with Mises-Huber hypothesis expressed as follows: 


1 
oO d= >=: 
re J2 (13) 


2 2 2 
(6,-6,) +(0,-0,) +(0,0,) +6- G54 +T.) 


In the case of thin-walled shells (of R/h > 20) under 
momentless load their state of stress can be considered two- 
dimensional, in this case the reduced stresses can be expressed 
by the principal stresses: axial ones, o, and circumferential 


ones, O» [4]: 
= 216 6. 
Cg = E FOTO 


Structural strength safety is determined by the condition 
of maintaining the reduced stresses in the shell plating smaller 
than the permissible ones. Value of the permissible stresses 
related to the structural material yield point R,, was obtained 
by assuming the safety factor n, > 1 of a value determined 
according to a standard, rules, or on the basis of experimental 
knowledge or empirical practice: 


(14) 


R, 
O nax = o per a (15) 

S 
Depending on a required reliability of a device,value of the 
safety factor is contained within the range of: n = 1.1 +20, and 
usually: n = 2.5 + 4.0 [1, 2, 3, 4]. Hence the structural strength 


safety condition takes the following form: 


= 22 y = a < _R. 
Oved max = MAX (VO, oO, O, Oy ) $S per= a 


S 


(16) 
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In momentless state of stress vessel’s shell plating is 
compressed simultaneously by uniformly distributed axial and 
circumferential forces. Unit compression forces in vessel’s 
casing are determined by the following relations: 

— axial force per unit length of vessel’s circumference: 
2 
es -R-p (7) 
2-m-R 2-m-R 


— circumferential force per unit length of vessel’s directrix: 
F 2-R-L- 

1252 =R>p a) 

2-R 2°L 


The principal axial stresses are determined by the 
relation: 


2 
gaT e E 9) 
2:m-R-h h 
The principal circumferential stresses by the relation: 
5 <2 RPL _R. a 
2-L-h h 


If vessel load changes occur slowly then dynamic and 
fatigue loads can be neglected, and the structural strength 
safety criterion for vessel under uniform external pressure is 
determined by the relation: 


| 2 2 
Oed max = MAX ( o, +O¿70, ‘o )5 
BR R, 


(21) 


D : h Pmax SO per = 


Ss 


CRITERION OF FORM 
STABILITY SAFETY 


Static stability of shell form means that there is an 
equilibrium state in a position neighbouring to the considered. 
The state of static equilibrium of form of axially symmetrical 
vessel under constant external load at momentless state of stress 
and linear elastic strain is described by the deformation form 
function w = w(z, @) which is solution of shell form stability 
equation. Depending on parameters of a vessel, loss of form 
stability may take mode of axial or circumferential buckling. 
Form deformation depends on values of geometric and material 
parameters of a vessel and mode of its loading. The parameters 
constitute arguments of differential equation of shell stability. 
Under certain simplifying assumptions the form stability 
equation of cylindrical shell takes the following form [4]: 


E-h @&w 1 Ow 
D-V'w+——-——+—-p-R-V* | — |+ 
R? e 2 Oz 
6 (22) 
Ow 
+p-R-V* 0 
R’-d@° 
where 
3 
D= Ehi - shell stiffness 
121.0 — v^) 
4 4 4 
Vine yg Le 
oz R? dz’ -ðp? Rİ Ao! 


Eq. is solved by assuming a priori a form of allowable 
function which satisfies relevant boundary conditions of 
deformation: 
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= a DE Z 
Ww=ws:sin -cosn-@ (23) 
where: 

w  - shell deflection 
m - number of longitudinal half-waves of shell deformation 


n - number of circumferential waves of shell deformation. 


On substitution of the allowable function to Eq. (22) the 
algebraic equation of deformations is obtained as follows: 


= 2 4 
i E nme] +n’ +E-h(n-mé) + 
RË? R L 
2 2 
—p-R Dp me) +n) T me) 4 211, C9 
2 E L 
-sin 


T-M-Z 
————-cos(n:p)=0 
7 "08(2-9) 


The pressure at which the equation is fulfilled, is determined 
by the relation [4]: 


(s+n) A de 


2 
A an 
mn Mu i ) 2 «Be a (25) 
s? 2 R 
2 +n? |-@’+n’) 
2 
where: 
PES 
=m ms 


The critical pressure at which a shell under uniform external 
load loses its form stability within elastic deformation range 
(i.e. suffers buckling), is equal to the smallest value of the 
pressures p,, , determined at the successively assumed integer 
numbers m of longitudinal half-waves and n numbers of 
circumferential waves, at fixed R/L and R/h ratios and the 
parameters E, v, at which relevant modes of form deformations 
take place, namely the critical pressure p,,: 


P.¿= MiN{P,, a} (26) 
where: : 
manEN. 
MINIMUM OF CRITICAL PRESSURE 


FUNCTIONAL 


Value of the critical pressure p, is equal to the global 
minimum value of critical pressure relevant to all possible 
modes of shell form deformations within permissible variability 
range of design parameters: 


= min ia (27) 
Per rard ° "L'h” i 

The presented method concerns vessels of specific 
geometrical and material features which have been a priori 
assumed to be functionally useful [5], and the below given 
range of their variability determines the allowable applicability 
space of the method: 
— the range of allowable values of the ratio R/L: 

0.25<R/L<1 (0.25, 0.5, 1.0}; 


— the range of allowable values of the ratio R/h: 
25 <R/h< 150 (25, 50, 150} 

— Poisson number v (for NW-steel and AL-5083 alloy): 
v = 0.3 and v = 0.33 

— Young modulus E (for NW-steel and AL-5083 alloy): 
E =2.1 - 10° and E = 0.72 - 10° [Mpa] 

— the number of considered axial half-waves of deformation: 
m= 1 

— the number of considered circumferential waves of 
deformation: n = 3 + 12. 


For the above given combinations of parameters the 
following was elaborated: 
— discrete values of the critical pressure: 
o RR 
bs min Py ma. ME 


m=1 ne{l +12} 


— diagrams of (continuous) functions of critical pressures [6], 
exemplified in Fig. 1 
— from the diagrams was read a number of circumferential 
waves at which the critical pressure p „reaches its smallest 
value, at constant values of the remaining parameters. 
Critical Pressure 


R/h=150 Steel s210 


Number of circumference deformation waves n 


crl (R/L=0.25) —e— Perl (R/L=0.5) 


la. P 
me Per] (R/L=1.0) == minPer 


Fig. 1. Example of envelopes of the minimum critical pressure values [6]: 


AS nm = 1È = 150, E = 2.1-10° MPa,v= 03) 


APPROXIMATION OF CRITICAL 
PRESSURE 


On the basis of the set of discrete critical pressure values 
an analytical approximation formula minimizing the sum of 
squares of deviations, was elaborated. In the case of vessels 
made of NW-steel the best approximation of critical pressure 
is represented by the following expression [6]: 


R 


-2.66 1.11 
Pa = 0.0806 + 414330, =] (2) [MPa] 


Pcr [MPa] 


RA i 60 


02 40 


. . . R R 
Fig. 2. Diagram of the relation p.,= f WL 


— for structural material: NW-steel 


In the case of AL-5083 alloy the best approximation of 
critical pressure is represented by the following expression 


6]: 
16] R -2.66 1.11 
p„=0.0262 + 139640, È) 8 [MPa] (29) 


Per [MPa] 


Pcr [MPa] 


02 40 


; R R 
Fig. 3. Diagram of the relation p, = f TIL 
— for structural material:AL-5083 alloy 


APPROXIMATION OF MINIMUM VALUE 
OF SHELL THICKNESS 


On the basis of Eq. and was obtained the set of discrete 
values of the converse problem consisting in determining 
geometrical parameters of vessel at an assumed operational 
submersion depth: 


(2)- f f Per)» E) 


On the basis of the set of discrete values achieved from the 
relation (30) an analytical formula was elaborated by applying 
the method of least squares of deviations. In the case of steel 
vessels the best obtained approximation of critical pressure is 
represented by the following expression [6]: 


(30) 


z = 144745 — 5538. ao 4 


yo 


R Sie ) 
— 139075- E + 3.670: ES 
L L 


where: 
Pa y 


given in [MPa]. 


24, o 
240 
160 
00 = 
£ x 
ac 80 
20 
0 
o 


a 
RM er ? 
Fig. 4. Diagram of the relation 2 f (p..F) 
— for structural material: NW-steel 
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Tab. 1. 


PARAMETER Variant I Variant II Variant III | Variant IV 


Assumptions Unit of measure Value Value Value Value 
Structural material A1-5083 Al-5083 
0 - tight; 1 - flooded 0 1 
Material density 2.1 Za 
Submersion depth 200 200 


Mass of apparatus 


Average density of apparatus $ ; : 1.5 1.5 


Effective displacement | i : 0.1 0.1 
R/L ratio ] 1 1 
Id/R ratio ] ; ; 0.5 0.5 
Exponent t >2 ] : ; 2.5 25 
Young modulus 210000 210000 


Poisson ratio | 0.3 0.3 


Safety factor 


Life time [year] 
Water density [t/m*] 


Foam density [t/m] 


Results of calculations Unit of measure 


Design pressure [MPa] 


Vessel diameter [m] 


Casing height [m] 


Shell plating thickness [mm] 


Vessel height [m] 


Bottom height [m] 


Bottom surface area 


Casing surface area 


Vessel surface area 


Volume of two bottoms 


Vessel volume 


Vessel displacement 


Mass of vessel shell plating 


Mass of foam 


Total mass 


Effective displacement 


Axial stresses 


Circumferential stresses 


Reduced stresses 


Permissible stresses 


Minimum strength-based thickness 


Minimum permissible thickness 


Corrosion allowance 


Allowance for plate thickness manufacturing tolerance 


Allowance for other stresses 


Technological plating thickness 


Mass of vessel after thickness correction 


Total mass after thickness correction 


Effective displacement after thickness correction 
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In the case of vessels made of AL-5083 alloy the best 
approximation of critical pressure, within the range of 
considered design parameters, is represented by the following 
expression [6]: 


R > R 0.4533 
$ 75.682 + 82.78: Por om (R) (32) 


where: 
p,,- given in [MPa]. 
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420 
> 80 
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0 40 
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Fig. 5. Diagram of the relation ae f (p..E) 


— for structural material: AL-5083 alloy 


RECAPITULATION AND EXAMPLE OF 
APPLICATION OF THE METHOD 


By substituting, to the buoyancy equation , the approximation 
formulae for V(R,, Ly La), S(R,, Ly L), HR, Ly p,,) and 
p,,(R/h, R/L), the balance is obtained of mass and displacement 
components of vessel, expressed explicite by its design 
parameters and design assumptions, as follows: 


(1-a)- pp V (RL, L) a S(R,L, L): 


-h(R,L, Pa)’ pm, Pr +O: Pw’ 
(33) 


(2 +V,+S(R,L,L,)-H(R,L, Pa) -D, =0 
a 

Solution of the design problem in question is reached 
by numerical solving the buoyancy equation supplemented 
by relevant assumptions on requirements and conditions for 
a given case. 


Example of application of the method 


It has to be designed a vessel made of AL-5083 alloy or 
NW steel, capable of accommodating apparatus of the average 
density of 1.5 t/m*and 100 kg mass, intended for operation at 
200 m submersion depth in the water of 1.025 t/m*density, so 
as to obtain the vessel’s effective displacement of about 100 kg. 
Alternatively it should be considered the case of a vessel of 
foam-filled bottoms, so as to obtain, but after flooding its 
interior, the same effective displacement of about 100 kg. 

The supplementary assumptions and determined values of 
the design parameters are contained in Tab. 1. 


Note: This research has been performed in the frame of 
the R&D project No. PR NR O R00 0047 08 financed by the 
Ministry of Science and Higher Education from its resources 
for the years 2009-2011 


NOMENCLATURE 


- effective displacement 

Young modulus 

- value of external concentrated load 

- gravity acceleration 

shell plating thickness 

- length of casing (cylindrical part of a vessel) 
bottom camber 

- mass of apparatus 

- mass of filler (foam) 

- mass of vessel shell plating 

- circumferential force per unit of vessel directrix length 
- axial force per unit of vessel circumference length 
value of safety factor 

- casing curvature radius 

- material yield point 

- circumferential curvature radius 

- axial curvature radius 

- radial variable 

vessel surface area 

- bottom surface area 

- power exponent 

vessel volume 

- axial variable 

index of vessel’s variant (tight or flooded) 
quadrature step 

- Poisson ratio 

- reduced stresses 

axial stresses 

- circumferential stresses 

- average density of apparatus 

- filler density 

- vessel plating density 

water density 

- dimensionless radial variable 

- dimesionless axial variable. 
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ABSTRACT 


The possibly accurate numerical prediction of the detailed structure of vortices shed from the tips of 
hydrofoils is an important element of the design process of marine propellers. The concentrated tip vortices 
are responsible for the propeller cavitation erosion and acoustic emission. The purpose of the project 
described in this paper was to develop the numerical method for prediction of the tip vortex structure. 
In the course of the project the numerical calculations were confronted with the results of experimental 
measurements. This led to creation of the specific method of construction of the computational grid and 
to selection of the optimum turbulence model. As a result the reliable method for the accurate numerical 
prediction of the concentrated tip vortices for different hydrofoil geometry and flow conditions has been 
developed and validated. This method enables elimination of the unfavourable phenomena related to the 
tip vortices in the course of the propeller design calculations. 


Keywords: marine propulsors, vortex generation, numerical methods, experimental techniques 


INTRODUCTION 


The concentrated tip vortices shed from the tips of the 
marine propeller blades are responsible for many unfavourable 
phenomena. Due to the strongly reduced pressure in their 
centres they generate cavitation, which in turn leads to 
destructive cavitation erosion of the propeller and the rudder, 
generation of the intensive acoustic signals and induction of 
vibration in the elements of ship structure (cf. Fig. 1). These 
phenomena should be accurately predicted at the design stage 
of the propeller, in order to eliminate them completely or at 
least to limit their intensity to a harmless level. This can be done 
by a sufficiently accurate and widely experimentally validated 
numerical method, which responds correctly to even small 
modifications of the propeller blade geometry and to small 
changes in the flow conditions. Such a method may be used 
for numerical optimization of the marine propeller geometry 
and its operating conditions. 

The first step in numerical prediction of the tip vortex 
cavitation is the possibly accurate determination of the tip 
vortex structure, 1.e. the spatial distribution of the flow velocity 
and pressure in the vicinity of the vortex, including its centre. 
The accurate determination of location and magnitude of the 
minimum pressure is of particular importance, because this 
pressure is the decisive factor in cavitation inception. The 
region in question is dominated by the high velocity transverse 
flow generated by the vortex. Consequently, the resultant local 
velocity vectors deviate strongly from the direction of the 
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Fig. 1. Cavitating tip vortices shed from the blades of the marine propeller 
model and interacting with the ship rudder (courtesy of the Ship Design and 
Research Centre CTO S.A. in Gdansk, Poland) 


general flow. This creates special problems for the numerical 
methods, for example the most frequently used methods 
based on numerical solution of the Reynolds Averaged Navier 
Stokes equations. These problems concern the adequacy of the 
computational grid structure and the accuracy of the turbulence 
models, which have been calibrated for much simpler flows. 
The accumulated practical experience of the past shows, that 


in most cases the numerical methods tend to underestimate the 
degree of concentration and the minimum pressure of the tip 
vortices. Consequently, their ability to predict accurately the tip 
vortex cavitation inception was regarded as questionable. 

The purpose of the research project described in this paper 
was to develop a computational method, which makes use of the 
special structure of the computational grid and which employs 
the specific turbulence models, especially suited for the vortex- 
dominated flows. This method is based on the combined 
numerical and experimental research, presented in detail in 
[1, 2, 3] and described briefly in the following section. 


DESCRIPTION OF THE NUMERICAL AND 
EXPERIMENTAL RESEARCH 


The research forming the basis for the development of the 
new method included experiments [4] and calculations [5, 6, 7, 8, 
9, 10]. The experiments were conducted in the cavitation tunnel 
of CTO S.A. and they were concerned with the measurements 
of the velocity field in three cross-sections of the tip vortex 
generated by hydrofoils of different geometry in different flow 
conditions (cf. Fig. 2). The measurements were performed using 
the Laser Doppler Anemometer. The geometry of the hydrofoils 
was selected in such a way that it resembled the blades of 
contemporary marine propellers. Three variants of the hydrofoil 
were tested: one with optimum distribution of the hydrodynamic 
loading along its span, one with hydrodynamically unloaded 
tip and one with hydrodynamically loaded tip. Each of these 
hydrofoils was tested at three different angles of attack: +2.5 
degrees, 0 and -2.5 degrees. It was expected that the variation 
of the hydrodynamic loading along the span of the hydrofoils 
may influence the process of formation of the tip vortices and 
it may lead to generation of vortices of different strength. The 
sufficiently accurate numerical method should be able to detect 
these small differences. 


Ey e q Fa 
Fig. 2. Hydrofoil model mounted in the measuring section 
of the cavitation tunnel of CTO S.A. 


The numerical calculations were arranged in such a way 
that they reproduced the experiments as closely as possible. 
The details of the computational flow domain, based on the 
geometry of the measuring section of the cavitation tunnel, 
are shown in Fig. 3. 

The calculations were performed using three different 
CFD codes (commercial codes Fluent and Comet together 
with the code Solaga developed at CTO S.A.) and in all six 
turbulence models (Spalart-Allmaras, k-w, k-e, RSM, k-w SST, 
k-e RNG). The accuracy of the results of calculations for such 
a complicated, vortex-dominated flow, depends on the detailed 


structure of the computational grid and on the adequacy of 
the applied turbulence models. The purpose of calculations 
was to develop a method for construction of the optimum grid 
for this specific type of flow and to select a turbulence model 
which is best suited to simulate the process of turbulence 
generation, transport and dissipation in this type of flow. This 
was performed on the basis of comparison of the experimental 
and computational results in three cross-sections of the vortex, 
located respectively 10 mm, 70 mm and 300 mm behind the 
hydrofoil tip, as shown in Fig. 4. 


Fig. 3. The discrete mesh for numerical calculations on the hydrofoils 
surface (top) and on the external surfaces of the flow domain (bottom) 


Fig. 4. Location of the cross-sections of the tip vortex 
for computations and measurements 
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The process of development of the optimum computational 
grid was conducted in three stages: 

- thefirst, basic grid was constructed according to the general 
recommendations included in the manuals of the respective 
CFD codes (cf. Fig. 5). This was a block-structured grid, in 
which the O-type grid was employed around the hydrofoil 
and the H-type grid was used in the remaining part of the 
flow domain. The total number of finite volumes was around 
1.250.000 and the size of the smallest volume close to the 

about 0.02 mm. 
| 


Fig. 5. Basic computational grid in the vicinity of the hydrofoil 


- the second grid was refined inside the cylinder covering the 
anticipated location of the tip vortex (cf. Fig. 6). Inside this 
cylinder all initially defined finite volumes were divided 
into 8 smaller volumes, rising the total number of volumes 
to about 1.940.000. 


a à 
A E 
Fig. 6. Computational grid refined in an arbitrarily defined cylinder 
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- the third grid was refined in the semi-automatic way, 
using the criterion based on the local values of vorticity, 
as calculated using the basic grid (cf. Fig. 7). The values 
of vorticity were tested inside the blue region shown in 
Fig. 7 and the volumes inside the red region in which the 
vorticity exceeded the prescribed value were modified. 
Each finite volume in the modified region was divided into 
16 smaller volumes, rising the total number of volumes to 
about 1.638.300. In order to have the refined grid in the tip 
vortex (core and vicinity), the streamwise voriticity has been 
selected as the criteria. The range of the streamwise vorticity 
for the grid adaptation is case dependent. In general, one can 
apply the maximum of the vorticity as the upper limit and 
30% of the maximum with the opposite sign as the lower 
limit. The lower limit is also important, because the counter 
rotating vortices are generated downstream the trailing edge 
in some cases. They are much weaker than the main vortex, 
but their influence on the tip vortex development and the 
diffusion rate is important. In case of the Comet code, the 
streamwise velocity was applied as the criteria for the grid 
adaptation. 


Dy 
ms 
Fig. 7. Computational grid refined on the basis 
of the local values of vorticity 


The distributions of the x-component of the resultant flow 
velocity and the x-component of the local vorticity calculated 
using the first grid, which were later used in the semi-automatic 
grid refinement, are shown in Fig. 8. In the same Figure the 
resulting grid structures in the consecutive stages: basic grid, 
adaptation | and adaptation 2 are shown. 

It is easy noticeable that the basic grid is coarse. Such effect 
could be avoided if the tip vortex location could be foreseen 
during grid generation and the grid would be generated fine 
enough at the very beginning. It is difficult or even impossible 
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Fig. 8. Process of grid refinement in the tip vortex area 


for the full 3D propeller analysis. In such a case, the coarse 
grid can be created and then it can be refined according to the 
local flow conditions. The local grid refinement (shown in 
Fig. 8) enables much better prediction ofthe local velocity and 
pressure gradients, which are very important for the cavitation 
prediction. It is shown that grid adaptation gives higher cells 
concentration in the tip vortex and such approach seems to be 
the most efficient. 


ANALYSIS OF THE NUMERICAL AND 
EXPERIMENTAL RESULTS 


The selected results of measurements and calculations 
included in this section are supposed to illustrate the agreement 
between them for the three consecutive stages of the grid 
refinement. Fig. 9 presents the axial velocity component 
in the cross-section 10 mm behind the hydrofoil, obtained 


0.04 


0.02 


Fig. 9. Axial component of the velocity - 10 mm — L+2.5 — Fluent (SST) 
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experimentally and calculated using Fluent for the tip-loaded shows the analogical results for the cross-section 70 mm behind 
hydrofoil (L) at the angle of attack +2.5 degrees. This hydrofoil the hydrofoil. Figs. 11 and 12 show similar results obtained by 
geometry in combination with the flow conditions should the program Comet. Figs 13 and 14 include the comparison of 
generate the most intensive vortex in the entire project. Fig. 10 the experimental and computational results for the transverse 
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Fig. 11. Axial component of the velocity - 10 mm — L+2.5 — Comet (SST) 


14 POLISH MARITIME RESEARCH, No 2/2010 


(vertical) velocity component at the cross-section 10 mm behind 


the tip-loaded hydrofoil at +2. 


5 angle of attack. 


In all cases computed by Fluent, the best agreement with 
experimental data was obtained for the “adaptation-2”. In the 
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Comet code, because it was no possibility to set the vorticity 
as the grid refinement criterion, the longitudinal velocity was 
applied. For such a criterion, which is not so sensitive as 
the vorticity, the results obtained on the refined grids differ 
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Fig. 12. Axial component of the velocity - 70 mm — L+2.5 — Comet (SST) 
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Fig. 13. Vertical component of the velocity - 10 mm — L+2.5 — Fluent (SST) 
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Fig. 14. Vertical component of the velocity - 10 mm — L+2.5 — Comet (SST) 


from those obtained with Fluent. In some cases (Fig.11- 
12), the results obtained on the grid for the first adaptation 
(adaptation-1) are a bit closer to the measurements. 


METHOD FOR NUMERICAL PREDICTION 
OF TIP VORTICES 


The method developed and tested in the course of the 
described research project is in fact the method for effective and 
optimum use of the commercial CFD codes for the specific task 
of accurate prediction of the tip vortices generated by hydrofoils. 
This method is built on two basic recommendations: 

- the computational grid may be initially constructed 
according to general rules applicable to the CFD code used. 
After performing the initial calculations using this grid, the 
region of the anticipated location of the tip vortex must 
be defined and the grid in this region should be modified 
(subdivided at least by a factor of 16), using the criterion 
based on the calculated x (streamwise) component of the 
local flow vorticity. Then the calculation should be repeated 
using the modified grid, producing the final results. 

- in the RANS methods applied for the vortex-dominated 
flows the turbulence model k-w SST should be used, 
because it is best suited for such flows and it gives the best 
results as far as the prediction of the detailed vortex structure 
is concerned. 


These recommendations have been proven to be effective 
for two different commercial CFD codes used in the project. 
Therefore it may be safely assumed that they are valid also 
for other commercial codes. Calculations of the tip vortices 
generated by hydrofoils, which are conducted in accordance 
with the above recommendations, supply sufficiently accurate 
results for further analysis of the tip vortex cavitation. 
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CONCLUSIONS 


The following final conclusions may be formulated on the 
basis of the research presented above and in [3]: 


- contemporary Computational Fluid Dynamics methods 
are fully capable of predicting the vortex-dominated flows 
with accuracy sufficient for further analysis of the vortex 
cavitation inception and development 


- the accuracy of numerical prediction of the vortex- 
dominated flows is significantly improved by the semi- 
automatic refinement of the computational grid in the region 
of the vortex, based on the pre-computed values of the local 
axial (streamwise) component of vorticity 


- such a local refinement of the computational grid ensures 
acceptable accuracy of the results while keeping the 
calculation time to a reasonable minimum 


- inthe sufficiently dense grid (in both directions: crosswise 
and streamwise) the diagonal orientation of the resultant 
velocity vectors does not influence adversely the results of 
calculations of the vortex-dominated flow 


- among the most popular turbulence models the k-w SST 
model seems to be best suited for calculations of the vortex- 
dominated flows. 
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ABSTRACT 


During designing the ship its designer tends to obtain as high as possible efficiency of ship propulsion 
system. To this end on certain ships additional elements such as: nozzles, half-nozzles or suitably profiled 
fins attached to underwater part of ships hull before screw propeller, are applied (sometimes they are 
intended for the mitigating of vibration resulting from operation of screw propeller). Another device which 
affects operation of screw propeller is streamline rudder capable of improving its efficiency (most of the 
transport ships is fitted with single screw propeller and streamline rudder placed behind it). In this paper 
the influence of streamline rudder on screw propeller efficiency has been examined with the use of numerical 
methods of fluid dynamics (CFD).The obtained results indicate a very favourable influence of the rudder on 
screw propeller efficiency. Research in this area is continued and its results concerning impact of particular 
geometrical parameters of rudder on screw propeller efficiency, will be presented in the future. 


Keywords: Streamline rudder, improvement of screw propeller efficiency, 
computational fluid dynamics (CFD) 


INTRODUCTION 


During ship’s design process one of the most important 
tasks is a design of its propulsion system which should ensure 
reaching the assumed service speed by the designed ship at its 
as high as possible propulsion efficiency. The crucial element 
of ship propulsion system is a screw propeller whose efficiency 
is decisive of overall propulsion efficiency. Geometry of screw 
propeller decides on its efficiency, hence it should be optimum 
one in assumed design conditions, i.e. screw propeller should 
develop demanded thrust at an assumed ship’s service speed. 
Apart from the optimum design, screw propeller efficiency 
can be improved by using various additional devices such 
as: nozzles, half-nozzles, or suitably profiled fins attached 
to underwater part of ship’s hull before screw propeller [5] 
(sometimes they are intended for the mitigating of hull plating 
vibration resulting from screw propeller operation). The 
additional devices are aimed at the improving of distribution 
(direction and velocity) of water stream flowing from around 
the hull towards the screw propeller. The devices produce 
a positive effect only in strictly defined conditions. As results 
from practice, they are applied very rarely. 

Another device which influences screw propeller operation 
is a streamline rudder capable of improving its efficiency. 
Majority of transport ships is fitted with single screw propeller 
and streamline rudder placed behind it. Therefore it seems 
purposeful to examine which way and in which conditions 
streamline rudder affects screw propeller operation, and how 
much screw propeller efficiency can be increased. 
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SCREW PROPELLER — STREAMLINE 
RUDDER SYSTEM 


The screw propeller — streamline rudder system is 
exemplified in Fig. 1. 


3300 


3450 


Fig. 1. The semi-spade rudder behind the screw propeller of B 573 ship [8]. 
(All the dimensions given in mm) 


A few dozen years ago already it has been observed that 
a blade rudder of airfoil profile not only improves ship’s 
manoeuvrability but also increases screw propeller efficiency 
and — in consequence — overall propulsion efficiency of ship 
[6]. The improvement of screw propeller efficiency results from 
the following reasons: 


— part straightening the flow behind the propeller (which leads 
to lowering the induced circumferential velocity) 

— amore favourable distribution of water pressure behind the 
propeller and reduction of vortices flowing downstream 
from the propeller, mainly from its boss. 


However apart from knowledge of the facts there is a lack 
of results of systematic investigations, e.g. model tests, and 
in the subject-matter literature only a general information 
can be found, that does not provide any basis to state which 
way geometrical features of streamline rudder affect screw 
propeller efficiency and how large increase of the efficiency 
would be possible. Yet the application of numerical methods 
of computational fluid dynamics (CFD) to modelling water 
flow around propeller and rudder has made it possible to 
assess qantitively and qualitively influence of rudder on 
screw propeller efficiency. Most publications in which results 
obtained by means of CFD methods have been presented, deal 
mainly with influence of screw propeller on ruddder or with 
investigation of flow around the hull, propeller and rudder [1, 
2, 3, 4 and 7]. In this paper preliminary results of numerical 
investigations concerning influence of rudder on screw 
propeller efficiency, are presented. 


SCOPE OF THE INVESTIGATIONS 


Within the frame of the conducted research project 
a numerical analysis of flow around screw propeller was 
performed by using Fluent system to determine the thrust 
coefficient Kt and torque coeffcient Kq as well as screw 
propeller efficiency n, for: 
— a free propeller 
— ascrew propeller accompanied with a rudder of HS VAMP73- 
20 profile in the position of no deflection 
— for a screw propeller accompanied with a rudder of 
HSVAMP73-20 profile in the position of 15° deflection. 


The investigations were performed for the hull of a B573 
ship (built by Stocznia Szczecinska), its screw propeller and 
streamline rudder (Fig. 1). 

As results of model tests on rudder / propeller interaction for 
the B 573 ship have been lacking, the obtained computational 
results were compared, only at the advance coefficient 
J = 0.70, with the values calculated on the basis of the 
propulsion prediction for free screw propeller [8]. 


Geometrical parameters of the screw propeller of B 573 
ship (Fig. 1) [9]: 


Diameter D 5900 mm 
Pitch ratio P/D 0.739 
Expanded area ratio A,/A, 0.578 
Number of blades Z 4 


Geometrical parameters of the rudder of B 573 ship 
(Fig. 1) [8]: 

The rudder was designed on the basis of HSVAMP73-20 
profile. Its contour and size are in compliance with the rudder 
installed on B 573 ship. 

The computational model of the rudder is shown in Fig. 2, 
and that of the screw propeller — in Fig. 3. 

The computation domain consisted of three co-axial 
cylinders (Fig. 4). The screw propeller was placed inside 
the small cylinder put inside the fore large cylinder. Rear 
faces of both the cylinders were located in one common 
plane and adhered the fore face of the rear cylinder. Such 
arrangement makes it possible to apply the mesh for free- 


propeller calculations by using the sliding mesh method (the 
small cylinder is a rotating element and the large cylinders are 
motionless). For the free-propeller calculations no physical 
element is placed in the rear cylinder. For the calculations of 
propeller and rudder together in the rear cylinder is placed the 
rudder model either with or without deflection, respective of 
a modeled case. 


4 


—— 
WR 

r 2 
SSS a a a 
a A gr 


Fig. 2. Computational model of the rudder without any deflection 


Perspective 


Fig. 3. Computational model of the screw propeller 
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Fig. 4. Computation domain for the set of screw propeller and rudder 
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RESULTS OF THE NUMERICAL 
CALCULATIONS 


In advance of the actual calculations of rudder influence on 


screw propeller efficiency, the following has been calculated: 


20 


Hydrodynamic characteristics of the screw propeller model, 
which were compared with results of experimental model 
tests. The propeller model geometry acc. P 355 storage 
propeller and results of its model tests are contained in the 
report [8]. 

Description of the numerical calculations of the 
hydrodynamic characteristics and their results are 
presented in [10]. 


Drag face of propeller 


2.06e+04 
1806+04 
9.45e+03 
8.80e+02 
-7.090+03 
-1636+04 
-2.48e+04 
-3.34e+04 
-4.20e+04 
-5.06er04 
-5.91e+04 
6.77 e+04 
-7.63e+04 
8.400404 


-9.34e+04 
-1026+05 


Driving face of propeller 


Fig. 5. Pressure distribution and streamlines for the set of propeller 
and rudder without any deflection, the advance ratio J = 0.7 
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e Hydrodynamic characteristics of the streamline rudder of 
B 573 ship, (Fig. 1). Results of numerical calculations of 
the characteristics are presented in [11]. 


The actual numerical calculations of screw propeller 
efficiency with accounting for impact of ruddder,were then 
performed. They were carried out for the real dimensions of 
the propeller and rudder installed on B 573 ship, Fig.1. 

In Fig. 5 an example pressure distribution and streamlines 
for both the propeller and rudder without any deflection is 
presented, and in Fig. 6 — for both the propeller and the rudder 
deflected by 15° (The complete set of pressure and velocity 
distributions as well as streamlines is given in [12]). 
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-1006+05 
-1.156+05 ty 
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Driving face of propeller 


Drag face of propeller 
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Fig. 6. Pressure distribution and streamlines for the set of propeller 
and rudder deflected by 15°, the advance ratio J = 0.7 
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Tab. 1. Data and results of the analysis performed by using the Fluent system. 


Screw propeller without rudder 


TIN] 


Kt [-] 


Q [Nm] 


10Kq [-] 


n, [-] 


1624159 


0.596514 


948528.9 


0.59046 


0.160868 


1365931 


0.501673 


839757.7 


0.52275 


0.305631 


1122907 


0.412416 


722395.3 


0.449691 


0.438109 


924202.1 


0.339437 


623854.7 


0.38835 


0.556719 


720384.5 


0.264579 


518421.4 


0.322718 


0.652745 


517910.9 


0.190216 


405355.8 


0.252334 


0.720215 


281523.2 


0.103397 


281425.6 


0.175188 


0.657871 


1 
2 
3 
4 
5 
6 
7 
8 


51307.51 


Screw propeller + rudder without deflecti 


TIN] 


0.018844 


Kt [-] 


146254.1 


Q [Nm] 


0.091043 


10Kq [H 


0.263667 


n, [+] 


1529143 


0.561616 


930528 


0.579254 


0.154387 


1295804 


0.475917 


802417.5 


0.499505 


0.303432 


354941.1 


0.130361 


298832.2 


0.186023 


0.781123 


Screw propeller + deflected rudder 


Jt] 


D [m] 


n [1/s] 


V [m/s] 


TIN] 


Kt [-] 


Q [Nm] 


10Kq [-] 


p [t/m*] 


0.1 
0.2 


5.9 
5.9 


1.5 
1:5 


0.885 
1.77 


1618327 
1388741 


0.594372 
0.51005 


950265.5 
839860.6 


0.591541 
0.522814 


998.66 
998.66 


na [-] 
0.159998 


0.7 


5.9 


1.5 


6.195 


380906.1 


0.139897 


315136.1 


0.196172 


998.66 


0.310697 
0.794896 


Results of the calculations of thrust, torque and efficiency of the propeller with accounting for the rudder, are presented in 


Tab. 1, and its hydrodynamic characteristics — in Fig. 7 through 10. 


Free propeller LOK q(J) 
@ Kq free propeller 
E Kq + rudder defle 
A Kq + rudder without deflection 


ected by 15 


Kt. 10Kq. eta 
10Kq 


Fig. 9. The torque coefficient Kq for the free propeller, the set of the free 
propeller and rudder without deflection and that of the free propeller 
and rudder deflected by 15° 
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Fig. 10. The efficiency n, for the free propeller, the set of the free propeller 
and rudder without deflection and that of the free propeller and rudder 
deflected by 15° 


Fig. 8. The thrust coefficient K, for the free propeller, the set of the free 
propeller and rudder without deflection and that of the free propeller 
and rudder deflected by 15° 


POLISH MARITIME RESEARCH, No 2/2010 21 


CONCLUSIONS 


As data concerning results of model tests are very scarce, 
to make comparison of the numerical analysis results with 
those of real model tests for free propeller is possible only for 
values corresponding with the advance ratio J = 0.70; relevant 
results of numerical calculations and model tests are compared 
in Tab. 2. 


| | model tests 


Tab. 2. 


numerical 
analysis 


model / analysis 


314800 


281523.2 


1.11820269 


Q 


228800 
0.564 


281425.59 
0.65787141 


0.81300354 
0.8573 104 


The comparison of results for only one value of the advance 
ratio J does not allow to assess if the remaining results are 
correct or not, and does not provide information on correctness 
of run of the characteristics. However this is the only way of 
assessing the performed analysis in view of the lack of data. 

Comparison of the efficiency, thrust and torque coefficients 
for free propeller, the set of propeller and rudder without 
deflection and that of propeller and rudder deflected by 15°, is 


performed in Tab. 3, 4 and 5. 


Tab. 3. 


Efficiency n, 


Free 
propeller 


0.160868285 


Propeller + 
rudder without 
deflection 


0.154387068 


Propeller 

+ rudder 

deflected 
by 15° 


0.159997698 


0.305630894 


0.303431839 


0.310696571 


0.65787141 


0.781122746 


Tab. 4. 


Thrust coefficient K, 


0.794895671 


Free 
propeller 


0.1 | 0.596513711 


Propeller + 
rudder without 
deflection 


0.561616382 


Propeller + 
rudder 
deflected 
by 15° 


0.594371651 


0.2 | 0.501672932 


0.47591672 


0.510050409 


0.7 | 0.103396532 


0.130361105 


Tab.5. 


Torque coefficient K, 


0.13989741 


Free 
propeller 


0.059045983 


Propeller + 
rudder without 
deflection 


0.057925422 


Propeller + 
rudder 
deflected 
by 15° 
0.05915408 


0.052274967 


0.049950537 


0.052281374 


22 


0.01751876 
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0.018602322 


0.019617241 


In numerical calculations an applied computational model, 
type and size of mesh, and quality of geometrical model of 
analysed propeller play a very important role. The applied 
computational model could differ from the real propeller 
because of very scarce data on propeller geometry, hence to 
assess a degree of similarity between model and real propeller 
is not possible. The qualitative analysis performed by means 
of CFD methods demonstrated that the working screw 
propeller accompanied with streamline rudder located behind 
it, is capable of developing a greater efficiency than the free 
propeller without rudder. 

Asa result of presently conducted investigations an impact 
of geometrical parameters of rudder on propeller efficiency 
will be determined. 
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ABSTRACT 


Fin stabilizers are very effective devices for controlling the ship roll motion against external wave-generated 
moments. Lift forces due to flow around fin with an angle of attack produce anti - roll moment. Therefore 
control of attack angle plays important role in reducing roll of ships. This paper presents results of using 
a combined neural network and PID for roll control of ship with small draught. Numerical results are given 
of around-fin flow analysis with considering free surface effect modelled by neural network and imposed 
to controlling loop. Hydraulic machinery constraints are also considered in the modelling. The obtained 
results show good performance of the controller in reducing roll amplitude in random seas. The approach 
can be used for any irregular sea conditions. 


Key words: Fin stabilizer; neural network; PID control; restoring force 


INTRODUCTION 


The aim of using roll stabilizing systems on ships depends 
on requirements of their mission such as the preventing 
against damage of cargo, increasing the effectiveness of crew 
or accuracy of positioning and handling weapons on naval 
vessels. 

Roll motion is one of the most important ship motions. 
A large number of papers and publications are available on 
modelling the ship roll motion. There are many ship mission 
parameters influenced by rolling, such as safety of cargo, crew 
and passengers, fuel combustion, weapon positioning etc. Most 
of the roll stabilizers create a moment against roll exciting 
moment to decrease roll angle amplitude. 

Fin stabilizer systems are found useful for roll damping. 
In this case, simulation of flow field around the fin and its 
performance has important role in research on decreasing the 
roll angle. 

Lifting characteristics of the fins have been reasonably 
well predicted by semi-empirical formulae by Whicker and 
Fehlner [1] and Pitt et al. [2]. Sheldahl et al. [3] presented 
also experimental data of NACA0015 foil obtained in a wind 
tunnel test. 

Hydrofoil lift characteristics are among thoroughly studied 
subjects and many papers and text books have been published 
in this regards (Newman. [4]). 

Numerical results of foil analysis including cavitation and 
2D and 3D free surface effects were published by Bal et al. [5]. 
They. presented also an improved numerical wave tank model 
for the hydrofoils by considering cavitations. 


Effects of free surface and wave generated by a strut on fins 
attached to it were also studied numerically and experimentally 
by Lee et al. [6]. 

In this paper the results of flow field analysis around fin 
stabilizer including free surface, presented by the authors [7], 
are used for modelling the fin in the controlling loop. According 
to the results, hydrofoil lift characteristics are influenced by the 
presence of free surface in case of a considered small naval 
vessel. 

There are many nonlinear and free-surface effects involved 
in modelling behaviour of fin stabilizer. The effects have 
been discussed in detail e.g. by Tristan Perez [8]. In dynamic 
modelling for controller design it is useful to include as much 
nonlinear and environmental effects as possible in order to 
achieve more accurate model without losing controllability 
of the system. 

There are many controlling methods which could be used 
to achieve a good performance of fin stabilizers, for instance: 
PID, fuzzy logic, neural network. There are many papers on 
using traditional PID controller in fin stabilizers. The work by 
S. Surendrana et al. [9] can be referred to as an example of 
using PID in controlling the fin stabilizer. They used the CFD 
for calculating fin hydrodynamic behaviour by solving flow 
around a 3D horizontal fin attached to a wall. 

In order to improve the performance of PID controller in 
fin stabilizing a combination of PID and modern controlling 
methods is sometimes used. For instance, the using of fuzzy- 
PID controller has been presented by Liang et al. [10]. 

In the present work the results of using NN-PID controllers 
in series arrangement are presented. The combined controller 
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in parallel arrangement was used for trajectory tracking of 
robotic manipulators by Wu et al. [11]. The performance of NN 
controlling algorithm was tested in four degrees-of-freedom 
by Zirilli et al. [12]. Moreover, Alarcin et al. [13] used neural 
network to control a rudder roll stabilizer. 

In the present work the controlling loop consists of: ship 
dynamic model, neural network model of fin hydrodynamics, 
PID controller and inverse neural network model of ship 
heeling moment. In order to increase accuracy of the results, 
constraints of a hydraulic mechanism used for turning fins are 
also added to the fin model. According to the results presented 
by the authors [7] lifting force of fins with two opposite and 
equal attack angles, placed on two sides of sample ship should 
not be equal. In order to obtain a desirable restoring moment 
at large roll angles, two different angles of attack should be 
used instead of equal and opposite angles. Therefore in the 
present work, instead of using two equal angles with opposite 
signs for fins, equal portions of restoring moment from each 
fin is considered a desired situation. Neural network models 
are trained by using classical Error Back Propagation (EBP) 
technique. Combined NN-PID controllers are designed and 
optimized by using the MATLAB and Simulink environment. 
This combination shows good results in reducing ship roll 
motion resulting from external disturbances. It is worth to 
mention that instead of modelling the waves, exciting force is 
modelled as a disturbance in controlling loop. 


GOVERNING EQUATION 


There are many approaches to developing the ship roll 
governing equation, such as linear, nonlinear, coupled and 
uncoupled methods. In this paper the model presented by 
T. Perez [8] is used to derive a model for controlling loop design. 
The uncoupled ship roll motion equation is as follows: 


1,4 +Bo + Co = Teontrol +Texiting (1) 
In Eq. (1), Tonua 18 the restoring moment produced by fins, 
and T iani is the roll exciting moment applied to the ship. 


There are analytical equations for calculating the lift force 
of fins in function of angle of attack, but in the present work 
numerical results accounting for the free surface effects, were 
used. 

It is important to note that the effective attack angle of fins 
may change because of roll motion and a new attack angle 
should be considered in finding the restoring moment. The 
formula for calculating the effective attack angle is given by 
Eq. (2) and (3): 

reo 
O, = arctan| — | ~ 
U 


no (2) 
U 


O stectine = a, ta 6) 


The force component which is directed perpendicularly to 
the flow velocity far ahead the foil constitutes the lifting force 
represented by the lift coefficient. Eq. (4) shows the relation 
between the restoring moment and lift coefficient used for 
calculating the moment: 


=YpAUC.r, (4) 


restoring 


MODELING 


The first step of controller design is modelling. There 
are two different elements which should be modelled in this 
problem: vessel and fin stabilizer. 
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Vessel modeling 


In this work a naval vessel model prepared by Blanke and 
Christensen [14] is used. The model was also implemented 
by Perez [8] and MARINTEK at NTNU through open source 
MATLAB tool box, called MSS. The model contains all 
coefficients for nonlinear ship manoeuvring in surge, sway, 
yaw and roll. Tab. 1 shows the ship’s main dimensions and 
properties. 


Tab. 1. Main particulars of the ship 
51.5 m 
8.6m 
2.3m 
362000 Kg 
47934*106 
23763*106 


8 m/s 

4.47 m 
1.53 m 
1.1m 
20.41 m 
3.36 m 


The model was modified in order to be combined with 
fin-stabilizer. The initial values are: surge, sway, roll and yaw 
velocities as well as roll and yaw angles. The inputs are: surge 
and sway forces as well as roll and yaw moments. For purpose of 
our task the sway and yaw moment is set to be zero and the surge 
force to be constant in order to simulate behaviour of the ship 
on straight course with the constant speed of about 8 m/s. 

The roll moment is summation of the fin-stabilizing moment 
and an exciting moment due to environmental effects. The 
exciting roll moment is simply modelled by combining the 
sinusoidal moments with different amplitudes, frequencies and 
phase angles. It should be noted that it was not an intention 
of the authors of this work to calculate a real exciting roll 
moment in waves. 

The model outputs are: surge, sway roll and yaw velocities 
as well as roll and yaw angles. Although the straight course at 
constant forward speed were only used for the simulation, the 
model is capable of being applied to simulating roll in turning 
or other manoeuvring scenarios. 


Fin-stabilizer modeling 


In the present work the fin-stabilizer consisted of two 
identical non-rectangular hydrofoils of a low aspect ratio, 
symmetrically placed on both sides of vessel. The geometrical 
properties of fins are specified in Tab. 2, and Fig. 1 shows 
schematically the shape of the fin. 


Tab. 2. Fin parameters 


NACA0015 
Tip chord [m]__ | 09 | 


Root chord [m] 
Span [m] 


Hydrodynamic effects of the fin on vessel’s behaviour are 
considered negligible. In other words it is assumed that, except of 
roll controlling moment, the fin does not have any other effect on 
vessel’s manoeuvring model. Because of the small draught of the 
vessel, the free surface effects on hydrodynamic fin coefficient, 


are important. Only the calm free surface effects are considered 
in modelling the fins hydrodynamic properties. Effects of waves 
could be later added to controlling loop as a disturbing noise. 


PP. >” o 


Ti oi 
eee ip or 


Fig. 1. Schematic shape of the fin 


These authors investigated the flat free surface effects 
on hydrofoil’s lift coefficients by using the CFD and RANS 
model presented elsewhere, [7]. A sample case of free surface 
shape around ship and fin is presented in Fig. 2. The results of 
hydrodynamic lift coefficient in function of attack angle of fin 
and roll angle of ship are presented in Fig. 3. For considering the 
fin’s hydrodynamics in control loop, a simple neural network 
model is trained to follow this diagram. The trained neural 
network model inputs are: fin’s attack angle and ship’s roll 
angle, and the output is fin lift coefficient. 


Fig. 2. Free surface profile around ship body 
and attached fin, at 8 m/s forward speed [7] 


There are some constraints imposed by the hydraulic 
machinery such as: the maximum mechanical fin angle and speed 
of changing the fin angle. As Fig. 1 shows, after the stall angle 
the fin’s lift coefficient falls down. Moreover, the presented 
NN model is not valid after stalling, therefore the stall angle is 
also added to the model as a constraint. Dynamic effects of the 
mechanical system used for fin angle changing are represented 
by asimple integration loop. The maximum mechanical fin angle 
and its maximum rate of change are assumed to be 25 degrees 
and 25 degrees per second, respectively, i.e. usual design values 
[8]. Fig. 4 shows the considered model. 

To calculate a correct effective angle of attack is also 
important. As mentioned before, the effective angle of attack 
refers to the angle between net input velocity of foil and flow. 


Fin 
Mechanical Calculate 
System Effective 


Attack ang. 


Desired 
Mechanical 
ang. 


Roll Velocity( 2 ) 


Effective 
Attack ang. 


The net velocity will be here the summation of the inflow 
velocity due to ship forward speed and that due to ship rolling 
motion. Therefore another block added to the model to calculate 
the effective angle of attack is based on Eq. (3). Fig. 5 shows 
the complete model of the fin. The input to the fin’s model is 
the desired fin angle, and the output is the fin lift coefficient 
accounting for the mechanical system dynamics and the flat 
free surface effects. Next, the fin’s roll moment is calculated by 
multiplying the lift coefficient and a calibrating gain calculated 
from Eq. (4). 

According to the results presented in Fig. 3, because of the 
free surface effects, two identical fins placed on two sides of 
the ship do not give the same roll moment even at equal angles 
of attack. Therefore to better model the fin stabilizer two fin 


models shown in Fig. 5, are used. 
Fin 


Lift coefficient 


Attack Angle [deg] 
—E— Roll 0 
—*— Roll 25 


—@— Roll -10 
—<— Roll 20 


—4— Roll 10 
—4— Roll 30 


Fig. 3. Lift coefficient in function of angle of attack 
and various ship roll angles, U = 8 m/s [7] 
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ang. 


Desired 
Mechanical 


ang. Mechanical 


ang 
Fig. 4. Block diagram of fin mechanical system 
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Cal. 
Gain 


( 3 )Roll Angle 


Fig. 5. Block diagram of fin model combined with mechanical system 
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CONTROL LOOP AND RESULTS 


The next step is forming the controlling loop. To accomplish 
this and overcome complexities of models different controlling 
methods such as PID and neural network are combined. The PID 
controller is used as a main controlling tool. For considering 
complex behaviour of roll moment in function of roll angle 
in case of large angles and forward speed, an inverse neural 
network model is trained by using the Error Back Propagation 
technique (EBP) based on the results achieved from the 
introduced ship dynamic model. The input for the trained NN 
model is the heel angle and the output is the heeling moment 
caused by the input heel angle. 

At the first step a simple controlling loop shown in Fig. 6, 
is used. The NN model of the ship heeling moment is placed in 
a closed control loop. Then the NN gains are optimized against 
positive and negative external pulses as disturbances, by using 
a genetic algorithm. By doing so, a more accurate and stable 
NN model for heeling moment is achieved. 

In the next step the PID controller was added to the NN 
as shown in Fig. 7. Then the PID gains were optimized by 
using the genetic algorithm to achieve optimized controlling 
performance in response to the same disturbances as before. 
The gains were used as the first estimation for the complete 
controlling loop containing fin models. 


Inverse 
Fin NN 
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Calculate 
Mechanical 
Attack ang. 


( 3 )Roll Velocity 


The heel moment resulting from the NN model was 
considered as the desired roll moment which should be 
interacted by the fin stabilizer. The desired roll moment should 
be divided in two portions, translated into desired mechanical 
fin angle and then passed to the fin model. To achieve this, an 
inverse NN model of fin lift coefficient, was trained. Equal 
moments acting on two fins were used, which was equivalent to 
equal portions of heeling moment. Other portions could be also 
used to achieve equal angles of attack or maximum efficiency of 
fin stabilizer. The arrangement for one fin is shown in Fig. 8. 

The output of the fin model will be an actual restoring 
moment produced by the fin stabilizer. In this step for this 
new arrangement the genetic algorithm is used to optimize 
PID gains against the same pulse disturbances. Therefore the 
steps of the complete controlling loop shown in Fig. 9, could 
be summarized as follows: 

1. the PID controller receives vessel’s roll angle and generates 

a suitable command; 

2. the PID command is translated to a desired total restoring 
moment; 

3. the desired total restoring moment is divided in two equal 
portions, one for each fin; 

4. the desired moment for each fin is translated into a lift 
coefficient by dividing it by a calibrating gain; 


Roll koment Roll ang. 


Ship Dynamic 
Model 


DeNormalizing 


Roll Moment Roll ang. 


Ship Dynamic 
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Fin 
Roll 
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Fin Model 
Fig. 8. The arrangement for translating the desired CL into the desired mechanical angle and actual fin moment 
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Fig. 9. The complete controlling loop 


5. each resulted lift coefficient is checked and constrained by 
the maximum fin lift coefficient; 

6. the lift coefficient is translated into a fin effective angle of 
attack by using the inverse NN model of fin lift coefficients; 

7. by considering ship roll velocity, a desired mechanical 
attack angle is calculated and passed to the fin model; 

8. the fin model gives an actual fin restoring moment by 
accounting for a fin mechanical system delay and free 
surface effect; 

9. the fins’ restoring moments are added to the environmental 
roll exciting moments and passed to the ship’s dynamic 
model; 

10. the ship dynamic model gives a new roll angle and roll 
velocity to be used in the next step. 


The model responses to two roll exciting moments are 
simulated. Fig. 10 shows two sample time history of the exciting 
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Fig. 10. Two sample exciting moments: case (a) - the upper figure, 
case (b) - the lower figure 
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moments. The upper figure is considered the case (a) and the 
lower figure - the case (b). 

Fig. 11 shows ship's roll angle time history in case of the 
forward speed of 8 m/s in presence of the exciting roll moment 
(a), with and without fin stabilizers. As showed, the ship roll 
angle is reduced by using the active fin models. It is important 
to note that the fin angles are not necessarily the same on both 
sides ofthe ship and each fin angle is optimized by considering 
the free surface effects. Fig. 12 shows the same results for the 
case (b). Although the complex model including the free surface 
effects was used for fin stabilizers, the designed controller was 
capable to reduce roll angle amplitude effectively. 


roll angle(deg) 


roll angle(deg) 


100 120 140 160 180 200 
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Fig. 11. Time history of roll angles for the exciting moment (a): 


the upper figure - without action of the controller, 
the lower figure - with action of the controller 
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Fig. 12. Time history of roll angles for the exciting moment (b): 
the upper figure - without action of the controller, 
the lower figure - with action of the controller 


CONCLUSION 


New fin models accounting for free surface effects were 
developed to be used in dynamic modelling of ships with fin 
stabilizers. By using the new model the fin stabilizer dynamic 
model was improved. A neural network model was trained on 
the basis of CFD simulation of hydrodynamic behaviour of 
a fin close to free surface. 

For the effective controlling of fin stabilizers combined 
PID - NN controllers were used. Two fins placed on both sides 
of the vessel were controlled in a different way to overcome 
unsymmetrical behaviour resulting from free surface effects. 
By doing so, more realistic and efficient behaviour of fins 
were achieved. 

A genetic algorithm was used to optimize behaviour of NN 
and PID controlling system. Good performance and response 
of the controlling system was achieved by using the controlling 
approach even to more complex fin models. 


NOMENCLATURE 


- Roll angle 

- Roll velocity 

- Roll acceleration 

- Ship service speed 
Immersed depth 
Submergence Froude number 
- Lift coefficient 

- Finarm 

- Pressure coefficient 

- Fin attack angle 

- Effective attack angle 
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L - Ship Length. 

B - Ship Breath 

D - Ship Draft 

A - Ship Displacement 

I, - Ship Moment of Inertia around Z axes 

I, - Ship Moment of Inertia around X axes 

KM - Ship meta-centric height from keel 

KB - Ship center of buoyancy from keel 

GM - Distance between ship’s center of gravity and meta- 
center 

LCG- Ship’s longitudinal center of gravity 

VCG- Ship’s vertical center of gravity 

PID - proportional-integral-derivative controller 

NN - Neural Network 

CFD - Computer Fluid Dynamics. 
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Entropy function application in the selection 
process of diagnostic parameters of marine 
diesel and gas turbine engines 
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Gdansk University of Technology 


ABSTRACT 


The article presents the method of analysing the diagnostic informativeness of the parameters 
characterizing gas dynamic processes observed inside the working spaces of marine diesel 
and gas turbine engines. The entropy function was used as the measure of indeterminacy 
of the identified set of engine unserviceability states. Based on numerical simulation 
experiments of the analysed gasodynamic processes, the amount of introduced diagnostic 
information was assessed and the most adequate parameters for the technical state of the 
engine were selected. These parameters compose the minimised set of diagnostic parameters 
which make it possible to assess unmistakably the technical state of the examined working spaces. 


Keywords: marine diesel and gas turbine engines; diagnostics; diagnostic informativeness; entropy function 


INTRODUCTION 


Diagnostic investigation is a basic element of the diagnostic 
process of each technical device [1]. A diagnostician, expert 
in the given branch, collects so-called primary information 
on the examined object from organoleptic examination, or by 
measuring parameters of the signals generated by the object. 

In the next step of the diagnostic activity the expert 
(diagnostician) performs a so-called measurement-based 
inference. He creates a set of symptoms, analyses it and 
formulates a preliminary diagnosis about the technical state 
of the machine based on his knowledge about past defects. 
Usually the created set of symptoms corresponds to a number of 
possible defects, therefore the diagnostician makes attempts to 
gain additional diagnostic information (symptoms), using more 
and more sophisticated (expensive and technically complicated) 
methods and means, to be able to eliminate less possible defects. 
An ideal situation is when only one possibility remains, being 
the final diagnosis characterised by a given probability [4]. 

The “bottleneck” of the diagnostic process organised in the 
above way, especially when it refers to technical objects of high 
complexity level, is the diagnostic knowledge resulting from the 
past experience of the diagnostician. This experience referring 
to the defect-symptom relations is gained from: 

e experts’ opinions, 
e experimental investigations on a real object: 

- active experiments — with the introduction of real defects 

to the machine, 

- passive experiments — with many years’ observation of 

a large number of examined machines of the same type, 
without interference into their technical state, 


e simulation experiments on specially developed computer 
codes — with possible modelling of different types of 
machine defects. 


Due to the dynamic development of the widely understood 
computer science and technology, in the next years to come 
the computer simulation methods can become very useful 
diagnostic tools. However, significant limitation in their 
practical application is the problem of experimental verification 
ofthe computer codes used for the simulation. The only way to 
confirm the credibility of the diagnostic simulation experiments 
is still the testing examination done on model rigs (or, if 
possible, on real objects), done with the aid of the traditional 
analogue technique. 


DIAGNOSTIC PARAMETER 
TOLERANCE RANGE 


Producers of marine internal combustion engines of 
both piston and turbine construction define a set of basic 
diagnostic parameters which can be determined by the user in 
sea conditions. They also give tolerance limits within which 
the values of these parameters should remain during engine 
operation. The parameters characterise the quality of engine 
operation and make it possible to assess its general technical 
state. When the value of any parameter goes beyond the given 
limits, it is a signal of inadmissible disturbance of the energy 
conversion processes taking place in the engine, and a threat 
of failure. 

In this situation the user should begin to localise the defect by 
making attempt to collect detailed data on the technical state of 
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particular functional modules. He has to analyse a large number 
of diagnostic parameters whose values, nearing the limits of 
the operating tolerance ranges, can be the symptoms of small 
technical state changes, characteristic for the “approaching” 
state of unserviceability. All this leads to a conclusion that 
precise assessment of the diagnostic tolerances determines the 
depth and quality of the formulated diagnosis [1]. 

The problem with selecting a set of adequate diagnostic 
parameters and limits of their operating tolerance is particularly 
complicated when analysing dynamic characteristics of the 
engine. An effective tool which can be used in searching an 
analytical solution to this problem is the concept of the distance 
between functions, known from the classical mathematical 
analysis [8] — Fig.1. 

The distance between two continuous functions X(t) and 
X (1) is defined as: 


t=t 
[xo -X (0) dt (1) 
t=0 

where: 

X(t) — known form of the function which maps the time- 
history of changes of the diagnostic parameter of the 
fully serviceable engine; 

X’(t) — known form of the function which maps the time- 
history of changes of the diagnostic parameter of the 
technically unserviceable engine. 


When the forms of the compared functions are not known, 
and the only information about them is given in the form of 
discrete values of the recorded dynamic time-histories, we 
can calculate the functional J(x, x’), bearing the name of 
metric, whose value for a given pair of elements {x, x’} of 
the compared time-histories is the distance between points x 
and x’. The metric defined in the above way can be used in 
two ways: 

e for quantitative assessment of the time-histories — as the 

averaged value [8]: 


t 


Nx, -x,) 2) 


q=] 


J(x,x')= 


e for qualitative assessment of the time-histories — ands the 
maximum value [8]: 


J nax (5 X') = max 


O<t<t 


' 
XX, 


6) 


where: 

x(t) — discrete value of the diagnostic parameter of the fully 
serviceable engine; 

x’(t) — discrete value of the diagnostic parameter of the 
technically unserviceable engine. 


In both cases the value of the metric will be a number 
which expresses a certain dimension of the distance between 
the compared sets. The geometrical sense of the distance is 
explained in Fig. 1. 

In order to perform a comparison analysis of a larger 
number of the recorded dynamic time-histories, given in 
different sets (units) of discrete function values, a concept 
of the referential metric is introduced, whose dimensionless 
value can be a general comparative factor for all analysed 
time-histories: 


(4) 
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If X, is the set of discrete values of the time-history of 
changes of the (reference) control parameter recorded during 
the acceleration of the analysed machine system without defects 
and X" is the time-history of changes of the same parameter for 
the system with a defect, then the value of the metric defined by 
formula (4) will be the value of the diagnostic parameter which 
characterises certain state of unserviceability, in a quantitative 
aspect. This value can be considered the measure of diagnostic 
sensitivity of the parameter. 


Mee Field J 


0 t T 


Fig. 1. Geometric interpretation of the distance between two functions 
recorded in the form of sets of discrete values 


AMOUNT OF DIAGNOSTIC 
INFORMATION 


The technical state of each module of the marine internal 
combustion engine is described by a relatively large set of 
diagnostic parameters (either directly measured or calculated), 
which could be initially determined after analysing its 
functional scheme. The analysis of all possible parameters in 
the diagnostic examination process is pointless and should be 
reduced to the optimal set of parameters which secures the 
control of the technical state of the engine and localisation of 
the identifiable (known and typical) engine defects. Therefore 
the hypothetical set of diagnostic parameters defined in the 
initial stage of the diagnostic process is to be minimised taking 
into account the following criteria: 

e maximal amount of the introduced diagnostic information 
about the unserviceability states, 
e diagnostic sensitivity of the examined engine type. 


An attempt is to be made to reach a state in which the 
finally verified set of parameters provides opportunities for 
identification of each of operational unserviceability states of 
the analysed construction unit. The ideal situation is when one 
diagnostic parameter characterises unmistakably one certain 
fault. However, as results from the performed investigations 
[3, 7], in case of marine internal combustion engines one 
diagnostic parameter corresponds, as a rule, to a number of 
possible unserviceability states. 

In order to perform a rational selection of diagnostic 
parameters, the degree of indeterminacy of the selected set of 
unserviceability states is analysed. According to Shannon, the 
amount of diagnostic information about the technical state of 
the engine which is carried by each analysed parameter can 
be assessed using: 

e unconditional entropy — as the measure of indeterminacy 
of the set of unserviceability states: 


k 
E(S,)= -5 Pp; log, Poi (5) 
i=l 
where 
S, — finite set of engine unserviceability states, 
k — number of possible unserviceability states composing 
the set S, 
P, — Probability of appearance of one of possible 


unserviceability states. 


e averaged conditional entropy — which makes it possible 
to calculate the decrease of the indeterminacy of the set 
of unserviceability states after one (successive) diagnostic 
parameter òx, is determined: 


E(S,/6x; )= Paj (ÒX; Esx, (S,) + Py OxJES (S, ) (6) 


where: 
8x; — diagnostic parameter value 


8x; =l]- òx; — value opposite to òx, 


P,,(9X;), P(x) — probability with which the parameter òx. 
reacts to engine unserviceability states, an 
the probability at which 6x; does not react, 

E, (S,), Es(S,) — conditional entropies of the unserviceability 

i set after determining the diagnostic parameter 

ôx, for the subset of states to which the 

parameter 8x òx. reacts and the subset of state 


to which 6x; does not react, respectively. 


The amount of diagnostic information about the engine 
unserviceability states S, which is carried by the parameter òx, 
can be calculated using the formula: 


Ixos, =E(S,) - ES, /8x;) (7) 


SAMPLE RESULTS OF EMPIRICAL 
INVESTIGATIONS 


The applicability of the proposed method to the process of 
selection of diagnostic parameters with the aid of the entropy 
function was verified using the results of past numerical 
experiments simulating gasodynamic processes realised in 


the working spaces of the selected marine engines: a medium 
speed diesel engine with pulsating turbo pressure charging 
system, and a three-shaft gas turbine engine with a separate 
power turbine [5, 6]. The mathematical models of processes, 
developed for this purpose, made it possible to simulate selected 
known and recognisable defects of the constructional structure 
of the examined engines, which are most frequently observed 
in the engine operation process. 


The marine diesel engine 


The states of engine unserviceability and the state of its full 
serviceability, which are used in the analysis, were obtained as 
a result of mathematical modelling of gasodynamic processes 
taking place in the working spaces of the marine engine. The 
analysed time-histories of gasodynamic parameters (temperature 
T, pressure p, and the speed v of the exhaust gas pressure peak 
amplitude displacement) were calculated for the selected control 
section of the exhaust gas outlet channel (the channel which 
connects engine cylinders with the turbo compressor turbine). 
Percent values of the referential metrics (diagnostic parameters) 
òx; which are given in Tab. 1, were assessed by comparing the 
standard time-histories of gasodynamic parameters recorded 
during the simulation of engine processes in full serviceability 
conditions, and corresponding curves recorded for the engine 
with the modelled operational E eae states Sp S,» 
Sa Saa Oe Sap Pup Sg ge Suite Sa a and S i4 


n6? ~n7? ~ n8? ~n9? ~ nlo’ ~ nll’? aip S 


Mutual relations between the finite set of engine working 
space unserviceability states S and diagnostic parameters 
òx, identifying those states can be clearly presented using 

so called diagnostic matrices — Tab. 2. It was assumed that 
if the diagnostic parameter value òx; exceeds the tolerance 
range limits by 10 percent or more (x, > > 10%) while reacting 
to the unserviceability state S; then “1” is placed in the 
diagnostic matrix entry situated at the crossing of the j-th 
row and the i-th column. If the parameter does not react to the 
unserviceability state — “0” is placed. The last matrix column 
contains the amounts of diagnostic information calculated 
using formula (7). 


Tab. 1. Metric values dx, for the modelled operational unserviceability states of the working spaces of SULZER engine 6AL20/24 type 


Operational unserviceability states 


Diagnostic 


parameter 


S,, ~ decreasing outlet valve active cross-section areas in cylinders 1, 2 and 3 by 5 %, 

S,, — decreasing outlet valve active cross-section areas in cylinders 1, 2 and 3 by 20 %, 

S „~ decreasing outlet valve active cross-section areas in cylinders 1, 2 and 3 by 30 %, 

S,, ~ decreasing outlet valve active cross-section areas in cylinder | by 5 %, 

S,; — decreasing outlet valve active cross-section areas in cylinder | by 20 %, 

S ~ decreasing outlet valve active cross-section areas in cylinder 1 by 30 %, 

S,, — changing outlet valve opening and closing angles in cylinders 1, 2 and 3 by +5° OWK, 
Sig — changing outlet valve opening and closing angles in cylinders 1, 2 and 3 by -5° OWK, 
S „~ changing outlet valve opening and closing angles in cylinder 1 by +5° OWK, 

S 1o ~ changing outlet valve opening and closing angles in cylinder 1 by -5° OWK, 

Sı ~ decreasing combustion chamber volumes in cylinder sections 1, 2 and 3 by 25%, 


Si) ~ decreasing combustion chamber volumes in cylinder sections 1, 2 and 3 by 50%, 
e decreasing combustion chamber volumes in cylinder section 1 by 25%, 

Sı, ~ decreasing combustion chamber volumes in cylinder section 1 by 50%, 

S „~ State of full operational serviceability of engine working spaces. 
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Tab. 2. Diagnostic matrix of the SULZER engine working spaces 64L20/24 type 


Operational unserviceability states 


Diagnostic 


parameter 


As results from the numerical data collected in the diagnostic 
matrix (Tab. 2), the diagnostic information maximum criterion, 
I = max, leads to the selection of the exhaust gas temperature 
as the best diagnostic parameter (I,, = 1.0). This makes it 
possible to reduce the diagnostic inference analysis to eight 
states, to which this parameter reacts, or seven remaining 
states, depending on the value which it reaches in the control 
process. 

The diagnostic parameter of similar applicability is the 
exhaust gas pressure in the outlet channel, for which the amount 
of the introduced diagnostic information about the analysed 
technical states of the engine working spaces is equal to 
[,, = 0.970. 

At the same time we can neglect the exhaust gas 
velocity which introduces the smallest amount of diagnostic 
information, I, = 0.355. This parameter reacts to all analysed 
unserviceability states by exceeding the assumed operating 
tolerance limits. 

In this case the selection of the adequate diagnostic 
parameters should be done using the criterion of control 
sensitivity of the examined engine type and ability to perform 
necessary measurements. Since the measurements of the 
working medium temperature reveal remarkable inertia, which 
forces the use of thermo-elements of an order of several tens 
of micrometers in diameter [9], it is the exhaust gas pressure 
in the outlet channel which should be selected for further 
analysis, as the measurement of this parameter seem to be most 
reasonable taking into account the marine engine diagnosing 
technology. 

Based on the results of the diagnostic investigations of the 
diesel engines used in the RP Navy, it was assumed that “k” 
engine working space unserviceability states S ,, where i= 1, ... k, 
compose a finite set of equally probable events p,, = 1/k [3, 7]. 
Therefore the formula which defines the unconditional entropy 
takes the form: 


i 1 
EG,)= -> es = log,k (8) 
i=l 
For the analysed system k = 15, hence the unconditional 
entropy: 


K(S_) = log,15 = 3.908 (9) 


The selected diagnostic parameters ST, dp and dv, collected 
in the diagnostic matrix 2, include only part of the diagnostic 
information about the technical state of the engine working 
spaces. Despite the fact that the condition which unmistakably 
determines the initial measure of indeterminacy of the analysed 
states is not met: 


E(S,) = I(T) +1(8p) + (èv) (10) 


We can name four groups of equally probable unserviceability 
states. For instance, when we analyse a set of diagnostic 
parameters (tab. 2) we can conclude that if the result: (ST, 
dp, dv} = (0, 0, 1}, appears in the diagnostic examination 
process, this may testify to the appearance of the following 
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unserviceability states: Sp S p S p S, and S ,, respectively. 
Any identical configuration of the results of the examination 
does not appear anywhere except these five cases. 

At the same time, equally probable appearance of the 
unserviceability states S S S s Sie Sup Sag Spo and Sio 
is fully defined by the result: (ST, dp, dv } = {1, 1, 1}. In 
a similar way the unserviceability state S, — {OT, dp, dv } = 
{0, 0, 1} can be identified. When the result (ST, dp, dv } = {0, 
0, 0} is obtained in the diagnostic examination process, this 
may signal full serviceability of the engine working spaces, 
provided that relation (10) is met. But this is not true, as the 
sum of the amounts of diagnostic information on the technical 
state of the engine working spaces introduced by parameters 
ST, Op, dv is equal to 2.325, and the unconditional entropy is 
equal to 3.908. To reduce the level of indeterminacy of the 
analysed unserviceability states (missing 1.583 of the amount 
of diagnostic information) we have to complete the set of 
diagnostic parameters in such a way that the condition (10) 
of the ability to make unmistakable distinction between these 
states is met. 

The missing amount of the diagnostic information about 
the technical state of the engine working spaces (additional 
diagnostic parameters) can be obtained from the harmonic 
analysis of the time-histories of pressure pulsations in the 
exhaust gas channel recorded experimentally on real objects 
(with real introduction of defects to laboratory engines) and 
systematic observations of the thermal and flow processes on 
a large number of examined objects of the same type, without 
interference into their technical state (such observations of the 
engines in operation on RP Navy vessels have been conducted 
for three years now). 


Marine gas turbine engine 


Operational conditions and small control sensitivity of 
a gas turbine engine installed in the marine power plant 
impose certain limits on the possible measurements of engine 
operation parameters. In numerous cases, complicated and 
extremely expensive measuring methods, only applicable in 
the engine test bed, are to be used. The sea conditions do not 
provide opportunities for introducing defects and changes in 
control procedures to confirm experimentally their impact on 
the quality of engine functioning and lifetime, and to assess 
diagnostic tolerances for indirect control parameters, which 
define a permissible course of the realised dynamic processes. 
The diagnostic tolerances should take into account uniqueness 
of engine production, i.e. individual characteristics of each 
engine unit, as well as different operating conditions and times 
of operation which result in different rates of aging and wear 
of elements. 

Determining the nature of mutual relations between the set 
of unserviceability states of functional engine modules and the 
set of diagnostic parameters which unmistakably identify those 
states is a basic goal of the diagnostic activity. A promising 
method of collecting so-called defect-symptom relations to 
be used in marine diagnostics of the gas turbine engines is 


modelling unserviceability states of turbine subunits with the 

aid of computer simulation [2, 5]. This method is extremely 

useful at the stage of designing the diagnostic system for 
engines used in marine conditions, as it makes it possible, in 

a relatively simple way, to: 

e determine initial static and dynamic characteristics of the 
engine, taking into account the deformation of their time- 
histories resulting from the action of external conditions, 
production differences, and irreversible aging and wear 
processes, 

e determine the set of possible unserviceability states for 
selected functional modules, 

e determine a minimal set of diagnostic parameters which 
unmistakably identify the modelled defects. 


When analysing the quality of functioning of the selected 
gas turbine engines, used on RP Navy vessels [4], the following 
unserviceability states were modelled based on the information 
about the defects recorded during their operation: 

Sı ~ pollution of compressor passages in SNC — assuming 
10% decrease in the efficiency and mass flow rate of 
the working medium compressed in the compressor 
and 2% compression decrease, 

S ~ Simultaneous pollution of the compressor passages in 
SNC and SWC — assuming 10% decrease in efficiency 

and mass flow rate of the working medium compressed 
in the compressors, along with 2% compression decrease 
for SNC and 7% compression decrease for SWC, 

— leakage in the compressor passages — assuming that the 
mass flow rate of the “lost” working medium is equal 
to: ED asad z 0.05m, wez 

— leakage in one of two air release valves behind SWC 
— assuming 10% decrease in SWC compression, 2.2% 
decrease in SNC compression, 4.9% increase of the 
mass flow rate of the working medium compressed in 
SWC and 2,8% increase of the mass flow rate of the 
working medium compressed in SNC, 

— defect of the automatic engine fuel supply system 
which results in forcing the engine acceleration process 
— the time of the working medium pressure increase, 
set at the same range behind SWC, was assumed to be 
reduced from 15 to 10 seconds. 


m3 


m5 


The object of simulation investigations was the process 
of acceleration of a three-shaft engine after introducing 
the modelled defects. The simulation aimed at determining 
the effect of changes of parameters that characterise the 
constructional structure of the analysed machine system on its 
dynamic characteristics. 

Tab. 3 collects per cent values of the metrics (diagnostic 
parameters), which were assessed by comparing the referential 
time-histories of the gasodynamic parameters recorded during 
the simulation of the process of acceleration of the fully 
serviceable engine, with corresponding time-histories for the 
engine with the modelled unserviceability states S p S, S, 
S|, and S, These values were calculated using the referential 
metric (4) in the same way as for the diesel engine, as was 
described in Section 4.1. 

This way a set of nineteen diagnostic parameters was 
obtained. This set is to be minimised using the entropy 
function: 
òT p Sp, ÖM p — temperature, pressure and mass of 

the air accumulated in the inter- 

compressor space, respectively 
éT,, — averaged temperature of the flow 
passage constructional material, 


St, — time of the air flow in the compressor 
space, 
Olaz — rotational speed of the high pressure 


rotor (reduced to the normal 
atmospheric conditions at SWC 


inlet), 

BH po OF upe — compression of the low and high 
pressure compressor, respectively 

SM, po OM pc» 5M,,,,, — 1s the mass flow rate through SNC and 
SWC, respectively 

Snipe ON ape — is the SNC and SWC efficiency, 
respectively 

OT 545 OP 5) — is, respectively, the air temperature 
and pressure behind SNC, 

ST", Sp’, — is, respectively, the air temperature 
and pressure at SWC inlet, 

ST OP a5 — is, respectively, the air temperature 
and pressure behind SWC. 


Tab. 3. Metric values dx, for the modelled unserviceability states 
of the three shaft ZORYA marine gas turbine engine UGT3000 type 


Operational unserviceability states 


Diagnostic 
parameter 


Firstly, the diagnostic matrix was created, at the same 
assumptions as for the earlier analysed diesel engine — matrix 1 
in Fig. 2. The last column in matrix 1 contains the amounts of 
diagnostic information carried by particular parameters and 
calculated using formula (7). Applying the criterion I, = max 
we can select fourteen diagnostic parameters which carry 
identical amount of information I, = 0.971. But as the first 
step, a parameter is to be selected whose measurement is the 
simplest from the constructional point of view. In this case the 
selected parameter was the pressure of the working medium 
accumulated in the compressor space, dp,,,. It is noteworthy 
that the diagnostic parameter which reacts (or not) to all 
unserviceability states carries no diagnostic information - for ST, 
and dp’,, we get, respectively: I (ôT) = 0 and I (Sp”,,) = 0. 

In the next selection steps, by proper restructuring of the 
diagnostic matrices and using available generalised calculation 
procedures [8], we can relatively easily select next parameters 
which carry the maximal amount of diagnostic information, 
provided that the first, second and next diagnostic parameter 
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were selected - matrices 2 and 3. In each case the diagnostic 
sensitivity is also analysed. 

The selection is continued until full diagnostic information 
about the technical state of the examined engine is obtained. 

Assuming that “k” engine compressor system unserviceability 
states, S, 1 = 1, ... k, compose a finite set of equally probable 
events p, = 1/k, which in cases of gas turbine engines is 
confirmed in practice [3,4], we get k = 5 for the analysed system, 
hence the unconditional entropy is: 

(11) 


E(S,) = log,5 = 2.322 
and 6m 


The selected diagnostic parameters Sp, 9T”,, a 
collected in matrix 4, carry, in total, full diagnostic information 
on the technical state of the compressor system. Then the 
condition is met which unmistakably defines the initial measure 
of indeterminacy of the analysed states: 


E(S,) = 1,(6Ppy) + 1, (T3) +1;(6Mpy) (12) 


Analysing the minimised set of diagnostic parameters we 
can conclude that if the result: {Sp,,,, ST”,,, dm,,,} = {1, 0, O}, 
appears in the diagnostic investigation process, then we can 
conclude about intensive pollution of blade passages in SNC 
(S). Any identical arrangement of the results of investigations 
does not exist anywhere except this only case. 

Simultaneous intensive pollution of blade passages in the 
two cooperating compressors (S_,) is completely defined by 
the result: {dp,,,, ST”, 5m,,,} = {1, 1, 0}. In a similar way we 
can identify the leakage in the air release valve (S,,) — (6p,,,, 


ST”, Mpu = {0, 0, 1}, and the defect of the engine fuel supply 
system (S,.) — {ðP u T” 6m,,,} = {0, 1, 1}. When the result 
{5p,4» OT ,,. dm,,,} = {0, O, O} is obtained in the diagnostic 


investigation process, it signals full serviceability of the engine. 
The modelled defect consisting in the leakage in the compressor 
passage (S ,) did not bring any substantial disturbance to engine 
operation which would result in exceeding tolerance range 
limits for the selected set of diagnostic parameters. 


CONCLUSIONS 


The method presented in the article consists in assessing 
the diagnostic informativeness of the measureable gasodynamic 
parameters of marine engines of both piston and turbine 
constructions. This assessment makes it possible to select a set 


of diagnostic parameters which precisely defines the technical 
state of constructional elements of the marine engine working 
spaces. The method provides opportunities for calculating the 
amount of diagnostic information of the selected parameters 
with the aid of the entropy function. As a consequence, it 
eliminates those diagnostic parameters which do not introduce 
valuable information on the technical state of the engine 
working spaces. 
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ABSTRACT 


In this paper the problem is discussed of supplying high power diesel engines with biofuels in the aspect 
of consequences which occurr during their operation. Attention was paid to advantages and difficulties 
associated with application of biofuels to self-ignition engines. A relevant research problem was 
characterized and research targets were proposed to make it possible to most favourably shape energy merits 
of the engines. Special attention was given to wear process of such engines and their failures resulting from 
application of biofuels. The following problems were also considered: calcium carbonate sedimentation, 
closing graphite platelets in structure of cylider bearing surface and surface defects of piston rings made 
of aluminium bronze. It was demonstrated that to assess operation of the engines is possible under the 
assumption that in a valuating approach their operation can be compared to a physical quantity represented 
numerically and measured with the use of the joule *sec unit. Conclusions resulting from the domestic and 
foreign investigations on power and reliability of such engines were also presented. 


Keywords: biofuel; operation; engine about automatic ignition; expendable 


INTRODUCTION 


In 21* century world economy has been dominated from its 
beginning by two problems: warfare against global warming-up 
and energy shortage, mainly electric one, and threat of running- low 
natural energy raw material resources, closely associated with the 
first. The themes are very closely connected to each other. Limitations 
associated with the global climate warming-up, consisting a. o. in 
reduction of emission limits of noxious gases such as NOx, SOx 
and first of all CO, have resulted in searching for alternative, pure 
and renewable energy sources [12]. However in applying renewable 
energy sources the today demand on electric power was to be taken 
into account. 

Until now in lower developed countries (e. g.Poland) electric 
power demand has been mainly covered by conventional coal-fired 
power plants operating on the basis of steam-turbine-driven electric 
generators. Unfortunately, the electric power stations appear to be 
the largest emitters of carbon dioxide, and their efficiency is rather 
not high. 

An alternative to the conventional electric power stations, as far 
as large power blocks are concerned, seems to be nuclear power 
industry, however after the Czarnobyl disaster its development has 
been effectively slowed down, at least in Eastern Europe. 

Development of nuclear power industry is associated with a high 
initial cost of electric power plants. The highest outlays should be 
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incurred during the building of nuclear reactor and before its putting 
into operation. Its investment cycle is also relatively long. Therefore 
it should be stated that power industry of the kind is rather intended 
for reach countries and its impact could be really felt at the end of 
the second decade. 

Pure energy sources are represented obviously by water - power 
plants as well as wind - power plants flourishing now in Europe. 
Experience gained at the run of century has showed that the wind- 
driven generator itself does not solve the problem in question as 
lack of wind is capable of leading to ,,black-out” and consequently 
paralysing life of the entire country as it was the case of Denmark 
and Holland. Presently we know that the wind electric power plants 
are to be bufferred at least in 50% by another power source being 
in a hot stand-by state. To this end, electric generators driven either 
by four-stroke self-ignition engines or medium- or large- power gas 
turbines, are commonly used. 

The self-ignition engine since its very beginning has been used 
for driving the electric generating sets, however it has practically 
been auxiliary, emergency, small power generators installed on 
objects requiring continuous supply or located far from electric 
networks, as well as generators on board of ships. Their output does 
not generally exceed 1 MW. 

Contemporary power industry surely expects solutions of electric 
generating sets of a much greater output — up to 20 MW and more 
- and available from one unit. The most modern diesel engines both 


four-stroke and two-stroke ones, designed on the basis of main 
engines intended for propelling large ships, are capable of coping 
with such tasks. Morover, a high efficiency of electric generating set 
driven by diesel engine supports attractivity of such solution. 

The large-power, two-stroke crosshead engine has an advantage 
which cannot be overestimated: it is relatively insensitive to fuel of 
a low quality, because of its specific design solution and adjustment 
to combusting heavy oils, of a rather not high quality, usually. This is 
especially important today in the situation of runing-low resources 
of natural liquid fuels and searching for alternative ones. For this 
reason research investigations on application of alternative fuels 
of vegetable origin, so called biofuels to self-ignition engines have 
been conducted for many years by domestic and foreign scientific 
research institutions!. 


GENERAL ASSESSMENT OF RESULTS 
OF THE INVESTIGATIONS HAVE BEEN 
CONDUCTED SO FAR 


Research work on applicability of biofuels to supplying self- 
ignition engines, of small power in particular, has been carried 
out by various scientific and scientific-industrial research 
centres. The activity has been triggered off due to growing 
awareness that [1, 2,3,8]: 

+ crude oil resources are lower and lower, 

e it is urgently necessary to limit emission of noxious 
compounds contained in exhaust gases, 

e bio-fuels are bio-degradeable. 


Itresulted in that biofuels which are methyl esters (sometimes 
ethyls) of fatty acids (FAME — Fatty Acid Metyl Esters) have 
achieved greater and greater importance in worldwide power 
industry. In Poland the methyl esters of fatty acids of rape oil 
(EMKOR) are mainly used. In this case their biodegradeability 
due to which they may be considered pro-ecological fuels, 
is of a special importance. The feature of biofuels has really 
directed interest of leading diesel engine producers to the fuels; 
it may be hence expected that investigations on consequences 
of application of such fuels will be continued. 

The conducted research [1, 13, 15, 16] deals mainly with: 
the methyl ester of fatty acids of rape oil, soybean oil, sunflower 
oil, their mixtures with diesel oil, and the methyl ester of fatty 
acids of palm oil, just recently. Investigations on application 
of animal fat waste as an alternative fuel are also carried out. 
It brings an additional advantage due to reduction of waste 
utilization cost. 

The vegetable fuels and their mixtures with diesel oil are 
characterized by a higher viscosity as compared with that of 
diesel oil, that makes engine starting-up process much more 
difficult. However engine starting-up features can be improved 
by providing heaters for biofuels characterized by high viscosity. 
The heating-up of such fuels can reduce their viscosity down 
to a level comparable to that of diesel oil. As results also from 
the invstigations, there are still important problems concerning 
the following items: excessively long combustion, reduced 
ignition lag, long injection process, slowed down evaporation 
of fuel droplets and worser fuel spraying quality. Moreover the 
fuels show high aggressiveness (chemical activity of highly 
acid reaction) against fuel piping and apparatus, that makes it 
necessary to apply special materials (stainless steel, sealings 
etc). Also, they do not have good lubricating features for 


1 


precision pairs of injection pumps and injectors, moreover they 

should have a low impact on change of physical and chemical 

features of lubricating oils in case of a leakage of the fuels to 
the oils resulting from loss of tightness of injection apparatus. 

For the above mentioned reasons further resarch is necessary 

to reach the following aims [8]: 

e improvement of starting-up features of engines from cold 
state, especially at low ambient temperatures, 

e determination of impact of vegetable fuels and their 
derivatives on: fuel piping corrosion, wear of injection 
pumps, injectors and the tribological system consisted of 
piston, rings and cylinder liner (resulting from friction 
and corrosion wear, sedimentation of lakes and calcium 
carbonate, cracks in piston heads and cylinder liners etc) 
as well as wear of crankshaft bearings, 

e determination of impact of vegetable fuels and their 
derivatives, in case of their leaking to lubricating oils during 
operation of engines, on physical and chemical features of 
the oils. 


The investigations have been performed so far, on 
applicability of the vegetable oils and their esters in question 
as alternative fuels for diesel engines, have yielded promising 
results. 

As a matter of fact from results of investigations on 
combustion of biofuels in engines the following was concluded: 
a decrease of their output, increase of specific fuel consumption 
(especially at higher engine speeds) and a decrease of efficiency, 
but simultaneously — decrease of exhaust smokiness, as 
compared with the case of diesel oil combustion in the engines. 
Exhaust gas temperature of biofuel-fed engine is lower as 
compared with that produced during diesel oil combustion. 

As stated from the performed investigations, the supplying 
of diesel engines with biofuels results in: 

e greater sedimentation of contaminations on surfaces of 
combustion chamber elements as well as on surfaces of 
spraying nozzles of injectors, 

e worsened operational parameters of the engines in question 
(lower output, higher specific fuel oil consumption, lower 
efficiency), sometimes higher smokiness of exhaust gas at 
lower rotational speeds of the engines, 

e increased emisssion of noxiuos exhaust-gas components 
such as: CO and C,H „ resulting in a lower content of CO,, 
but also lower content of NO, in exhaust gas at higher 
engine loads. 


The higher content of hydrocarbons (C,H, ) in exhaust gas 
of the engines combusting biofuels results mainly from that the 
flame ionization detectors of C,H „ analyzer do not seperate, 
from the hydrocarbons, aldehydes which are in a much greater 
amount in biofuels than in diesel oils. 

Along with growing engine wear the differentiation of 
emission of the above mentioned compounds increases. When 
analyzing the presented research results it should be observed 
that they are promising as compared with the results obtained 
from investigations of gas and steam turbines used in electric 
power stations. The quantative data given in Fig. | and 2 and 
in Tab. 1 support the statement. 

As illustrated by the presented data, in spite of the very 
promising prospects of the application of self-ignition engines 
to driving the electric generating sets it should be however 


Already in 1900 Rudolf Diesel anticipating that in the future crude oil and coal tar may be in shortage, 


presented, during Paris exhibition, an engine combusting arachis oil. For many years such fuels have not 
been used as the fuel obtained from crude oil was (and is still) less expensive 
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stated that the results obtained during the investigations have 
been performed so far on combustion processes in biofuel-fed 
engines, could not fully reveal their operational merits, as they 
have dealt only with energy aspects. No complex investigations 
concerning durability of main tribological units of the engines 
have been conducted with accounting for technical diagnostics, 
endoscopy based in particular. Proposals of such investigations 
are presented in this paper. 


CO, [g/kWh] 
MDE -diesel engine 
™ GCC -gas turbine of combined cycle En 
M GT -gas turbine 


ST 


— steam turbine 


Gas Fluid fuels Coal 


Fig. 1. Specific emission of CO, contained in exhaust gases for various 
types of electric power stations and mineral fuels [14] 
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Fig. 2. Comparison of efficiency of typical driving units used in electric 
power stations [14] 


Summing up, further development of drives based on self- 
ignition engines should be expected in power industry, and the 
crucial aspect deciding on reaching the success will be their 
reliability and reliability of electric power station as a whole. 


PROBLEM OF WEAR OF BIOFUEL-FED 
ENGINES AND THEIR TESTS 


Initial remarks 


The promising results ofthe investigations associated with 
determination of energy merits of biofuels proposed for self- 
ignition engines, have made that also research on intensity of 
degradation process of structures of such engines have been 
undertaken in order to determine a guaranteed time to failure 
of their particular constructional units. The results of diagnostic 
tests performed by MAN Diesel Co, have been obtained so far, 
showed that the energy parameters of ,,bio-engine” are kept 
in the range of contractual values, but values of the structural 
parameters of its working space change significantly due to 
calcium carbonate contamination — Fig. 3 [15, 16]. From 
the performed microscopic tests of a cylinder liner it was 
additionally concluded that the closing of graphite platelets 
in its cristallic structure took place. (Fig. 4), that may within 
a short time result in growing rate of wear of cylinder bearing 
surface - and consequently - worsening engine performance 
and efficiency of realized heat-and-fluid-flow processes [15, 
16]. Many defects of sliding coating (aluminium bronze) of 
piston sealing rings of cylinder systems (Fig. 5) confirm the 
observation. 


Formation of calcium carbonate sediments on 
surfaces of combustion chamber elements 


During a routine overhaul of internal spaces of one of MAN 
Diesel biofuel — fed engines a large amount of light-beige 
contaminations (,,chalk-like” hard sediments) were found on 
piston head surface as well as cylinder head from the side of 
combustion chamber, in all cylinder systems of the engine. 
Similar sediments were also found in internal spaces of the 
waste-heat boiler [15]. As preliminarily stated the cause of 
its occurence was an accumulated excess of non-used (non- 
combusted) additions to lubricating oil, a.o. detergent, in the 
form of acid calcium carbonate Ca(HCO,), which, in high 


Tab. 1. Comparison of features of selected technological methods for electric power production [14] 


Initial 
cost/kW 


Investment 


Technology period 


Piston engines small very short 


Operational 


cost Fuel price 


Regulation 
risk 


low high 


Coal large long 


medium medium high 


Gas turbine - 


combined cycle shon 


medium 


low high low 


Nuclear very large long 


medium low high 


Hydro 


large long 


very low null high 


Wind small short 


very low null medium 


Fuel cell small very short 


medium high low 


Photovoltaic | very large very short 


38 POLISH MARITIME RESEARCH, No 2/2010 


very low null low 


Fig. 3. View of cylinder space of the engine working on a biofuel produced on the basis of crude (not processed)palm oil. 
a) Calcium carbonate (CaCO,) layer sedimented on piston head. 
b) Contamination layer sedimented on the head from the side of combustion chamber in the vicinity of outlet valve and fuel injector 


200 um 


——y 


Fig. 4. Enlarged view of cylinder bearing surface of the engine working on a biofuel produced on the basis of crude (not processed)palm oil. 
a) Cylinder bearing surface — opened graphite platelets of good tribological qualities. 
b) Cylinder bearing surface — partly closed graphite platelets of worser tribological qualities 


Fig. 5. Ring part of piston — the engine working on a biofuel produced on the basis of crude (not processed)palm oil. 
a) Sliding surface of piston sealing ring — defects of aluminium bronze layer in the vicinity of piston-ring joint. 
b) Sliding surface of piston sealing ring — defects of aluminium bronze layer 
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temperatures, releases calcium carbonate CaCO, water vapour 
H,O w and carbon dioxide CO,: 


Ca(HCO,), > CaCO,+H,O+CO, (1) 


As deemed, it was caused by mutual interaction (an 
inappropriate combination) of base number of lubricating oil 
(TBN40) and a supplied fuel not containing sulfur. On the 
other side there is no commercial production of a lubricating 
oil of a base number smaller than TBN40. As separated hard 
particles of calcium carbonate sediments disturbe oil film 
as well as stable work of piston rings, they contribute in 
accelerated friction and errosion wear process of the cylinder 
bearing surface. Knowing that an inceased charge of cylinder 
oil could not improve the situation, but on the contrary it could 
accelerate rate of CaCO, sedimenation, one decided to limit the 
charge of lubricating oil (supplied to cylinder bearing surface 
through two lubricators) to the value of 1.05 g/kWh. According 
to MAN Diesel Co statement, such value should be generally 
determined on the basis of the relation of F x S%, where the 
empirical coefficient F takes values from 0.26 to 0.34 (and S% 
- percentage content of sulfur in fuel oil). Positive results of the 
lubrication at the lower coefficient have been demonstrated by 
tests and such values are recommended today. Generally, it can 
be concluded that the cylinder oil charge should not be taken 
lower than 0.8 g/kWh. From the producer’s experience have 
been gained so far it results that CL-DX405 oil of a reduced 
base number TBN40, reduced acid calcium carbonate content, 
at maintained high washing capacity, is the best adjusted 
cylinder oil for biofuel -fed engines (as well as for engines fed 
with fuel oil having sulfur content smaller than 1.5%). 

One may expect very soon that producers of lubricating 
means for low-speed engines will undertake intensive research 
on a new composition of lubricating oil of a TBN value close 
to 20 mgKOH/g. It results not only from demands put by 
biofuels but first of all from new legal regulations dealing 
with sea environment protection (Anex VI to Marpol 73/78 
Convention) in which some navigation regions were specified 
(a.o. the Baltic Sea, selected regions of the North Sea, region 
of La Manche Channel etc) where only use of fuel oils of sulfur 
content smaller than 1.5 %, with its target value lowered to 
0.5%, is permitted. 


To other 
cylinders 
A 


Because of the so small amount of lubricating oil applied to 
cylinder liner, it must be very precisely spread over the entire 
cylinder bearing surface in order to ensure a high quality of all 
realized functions of lubricating oil, namely: 

e limitation of adhesive wear by eliminating direct metal-to- 
metal contact of faying surfaces of piston rings and cylinder 
bearing surface, 

e limitation of corrosion wear by neutralizing acid combusion 
products, first of all sulfuric acid, in the form of calcium 
carbonate, 

e ensurance of stable (durable) oil film between piston ring 
and cylinder bearing surface, 

e limitation of sedimenting rate of contaminations in engine 
working spaces (combustion chambers), 

e ensurance of proper sealing the rings in piston grooves. 


Therefore a perspective, effective way to further control 
of oil consumption for lubrication of cylinder liners is the 
supplying of oil in an amount proportional to engine load 
(i.e. fuel oil charge delivered to cylinders) and proportional to 
percentage content of sulfur in fuel oil. In older solutions the 
functions have been realized by using a mechanical system 
connected with engine load indicator, today — by specialty 
electronic systems, e.g. the Alfa—lubrication system for MAN 
Diesel engines — Fig. 6, intended for controlling value of 
cylinder oil charge. 

However two basic conditions must be fulfilled: 

1. the cylinder oil charge cannot be smaller than a value 
which limits tribological wear of friction surfaces to 
a minimum, 

2. amounts of alkalic anticorrosive additions and detergents 
should be limited to a minimum necessary to neutralize 
formed acids and keep piston rings clean (ensurance of 
a sufficient washing capacity). 


In the case of the engine fed with sulfur - free biofuel the 
cylinder oil charge should be proportional to engine’s load 
with taking into account measurement results of contamination 
of the lubrication oil by Fe-particles. This makes it possible 
to estimate rate of wear of cylinder bearing surface and on 
this basis to correct an optimum charge from lubricator. As 
determined from service measurements of wear of cylinder 


To other 
cylinders 


Crankshaft position 
Engine load 


Fig. 6. Schematic diagram of the Alfa-lubrication system for controlling cylinder oil charge in MAN Diesel engines [14] 
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liners of MAN Diesel two-stroke engines, the mean rate of wear 
of cylinder liner was not greater than 0,1 mm/1000 h [15,16] 
(the data concerned the heavy-oil-fed engines). 


Closing graphite platelets in outer layer 
of cylinder bearing surface 


As concluded from results of diagnostic tests have been 
conducted so far on large power diesel engines, lubrication 
conditions can be improved to a certain extent by allowing 
controlled corrosion process of cylinder bearing surface to 
happen [15]. Local corrosion cells occurring in the form of 
microscopic ,,pocketlets” filled with lubricating oil support 
hydrodynamical lubrication process. Complete elimination of 
erosion and corrosion processes may lead to a loss of honing 
structure and polishing effect of cylinder bearing surface 
(Fig. 7). 

It usually results in closing the graphite platelets within 
cristallic structure of outer layer of cylinder bearing surface, 
that disturbes the building process of essential oil film. 
Consequently, it leads to increasing consumption of lubricating 
oil, scuffing piston rings in cylinder liner, and even seizing 
effects in PRC (Piston-Rings-Cylinder) system. 


P Oep 


Fig. 7. Polished (mirror) cylinder bearing surface 


4. Wear-affected operation of engines 


Operation of diesel engines consists in the transforming 
and transferring of supplied energy [4, 5]. It becomes worser 
and worser along with growing wear of the engines. To assess 
their operation is hence deemed necessary. 

Doing so, one should take into account that the transformation 
of energy in the form of heat or work, which occurs in working 
spaces of every diesel engine, should be such as to make the 
work to be done within a given time, as large as possible, or to 
do the work as fast as possible. In practice it is also important 
to release, during combustion process, as much heat as possible, 
and to reduce waste heat to a possible minimum. If due to 
wear of the engine to obtain such most favourable energy 
transformation is not possible, then the operation of the engine 
is considered incorrect and it is assumed that the engine is then 
in the state of partial serviceability [4, 5]. 

As the operation of such engines consists in the transformation 
of the energy E into the form of the work L and heat Q it can 
be hence generally interpreted as follows: 

t 
D= fear Q) 
where: i 
D — engine operation 
E — transformed (obtained) energy which makes realization 
of the task Z possible, 
t — duration time of transforming (using) the energy E 


In the most general case, E(t) may be considered the 
internal combustion gas energy which makes it possible to do 
the useful work L (1) by the engine. The energy is obviously 
obtained from the chemical energy contained in air-fuel mixture 
during fuel combustion in working spaces of every engine. 
Transformation of the chemical energy during operation of 
the engine results in the thermal and mechanical loading of 
the engine, that consequently results in growing wear of the 
engine and energy losses associated with the wear [4, 10, 
11]. To conclude about serviceability of a particular engine is 
possible after calculating value of its operation from Eq. (2), 
which, in the proposed interpretation, is equated to a physical 
quantity having the measure unit called joule-second. In order 
to determine the area of operation (D) it is obviously necessary 
to know a functional relation of energy and time, i.e. E = f(t). 
The so-understood operation determined by Eq. (2), can be 
presented in the coordinate frame ,,E—t”, hence also in the 
form of a diagram called diagram of operation [4, 5]. In view 
of that D = f(E, t) the operation of devices can be presented in 
the coordinate frame ,,D—E-*” [4, 5]. 

In the case of an arbitrary self-ignition engine, the 
transformation of the chemical energy (contained in air-fuel 
mixture produced in combustion chambers) into thermal one 
and then mechanical one, makes it possible to generate the 
crankshaft torque (M,) at a given rotational speed (n) [5]. 
Therefore the engine operation described by Eq. (2) may be 
interpreted as the energy transformation into the form of the 
useful work L, and expressed by the formula: 


D, = [L.@dt =27 fncom, (t)tdt (3) 
0 0 


As results from the presented considerations, it is 
reasonable to analyze not only the work L, but also the 
operation (D) of such engines, in this case understood as the 
energy transformation in the engine, which leads to achieving 
a necessary amount of the useful work (L,) within a definite 
time (t). It makes it possible to determine if the possible 
engine operation (D,,) in given conditions is at least equal 
to the demanded operation (D,,) necessary for realization of 
a given task Z. 

The task (Z) to realization of which the combustion engine 
is adjusted in the phase of its designing and manufacturing, can 
be realized only when the following inequality is satisfied: 


Dyu > Dy (4) 
i.e. when: 


tu 2 ty, and simultaneously E,, > Ewy 


where: 

tu — possible operational time of engine 

ty — demanded operational time of engine 

E,, — energy which can be transformed by engine 

E,, — energy demanded (necessary) to realize the task Z (energy 
which must be transformed to realize the task Z). 


It means that by analyzing energy merits of combustion 
engines in a general sense with the use of Eq. (2), one can 
consider serviceability of the engines in the following variants, 
acounting for their wearing: 

da a 
= ty» and simultaneously E, > E 
w» and simultaneously E,, = E 


and simultaneously E,, = Ey 
w 
w 
w 


(5) 


> ty» and simultaneously E,, > E 
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In the particular case when in assessing the engine operation 
itis important the definite work L, to be done by the engine, then 
the relation (5) should be presented in the following form: 


tu = tw» and simultaneously L y = Low 
tu = tw» and simultaneously L u > Lew 
tu > tw» and simultaneously L y = Ly 
í 6 
tu > tw» and simultaneously L y? Low (6) 
where 
tu — possible operational time of engine 
ty — demanded operational time of engine 
Lu — useful work which can be done by engine 
L w — useful work which is demanded (necessary) to realize 


the task Z (useful work which must be done to realize 
the task Z). 


As a rule it is important to do the useful work L, within 
a definite time. Then to assess the engine operation it is 
necessary to know how fast the work can be done. A physical 
quantity which contains such information is the effective power 
N, as it determines, in case of any machine, how fast a given 
work will be (or can be) done within a given time. Hence in 
the case of assessing the engine operation when it is important 


the engine to develop a definite power (N_), the relation (6) can 
be presented as follows: 

tu = tw» and simultaneously N y = Nay 

tu = tw» and simultaneously N y > Nay 

tu > tw» and simultaneously N y = Ny 

tu > tw» and simultaneously N y > Nay (7) 
where 
tu — possible operational time of engine 
ty — demanded operational time of engine 
Nou — effective power which can be developed by engine 
Nay — effective power demanded (necessary) to realize the 


task Z (effective power which must be generated to 
realize the task Z). 


In the case when the inequality is satisfied: 
D, < DW (8) 


then it should be concluded that the engine in question is 
damaged and not capable of realizing the assumed task Z. 

In the case of considering the situations (6) and (7) the 
inequality (8) should be expressed as follows: 


La = Loy or Nos > New (9) 


Therefore in order to determine if a given engine is 
serviceable to realize a given task, its operation should be 
identified [10, 11]. As the operation depends on a degree 
of engine’s wear, its identification must be connected with 
identification of engine’s wear, e.g. in a way presented above 
in this paper. 


RECAPITULATION. FINAL REMARKS 
AND CONCLUSIONS 


Features of diesel oil and considered alternative fuels are 
very different. Therefore it should be expected that self-ignition 
engines must be subjected to a special adjustment, or even 
modification, depending on a kind of fuel delivered to them. 
To adjust the engine fuel system to features of vegetable fuel 
oils of and their mixtures with diesel oil is necessary because 
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of their high viscosity as compared with that of diesel oil. One 

of the possible ways of decreasing the viscosity of biofuels and 

improving their spraying quality is to apply a heating-up system 

to the fuel oils, that makes it necessary to supplement the fuel 

system, not only with the fuel - oil heater but also a viscometer 

for measuring fuel oil viscosity as well as a viscostat - for 

controlling flow of heating medium through the heater and 

maintaining this way fuel oil viscosity on a required level. 
Biodegradability of the vegetable fuel oils is a very positive 

quality. However to be capable of serving as fully useful fuel 

cin they should: 
ensure similar energy qualities as compared with diesel 
oil, 

e do not cause an excessive wear of engine and consequently 
reduce life time of engine, 

e be a source of only low quantities of environmentally 
noxious substances (CO,), 

e ensure a low specific fuel oil consumption by engine, 

e facilitate recycling of raw and packing materials associated 
with their production. 


Choice of an alternative fuel oil obtained from fatty acid 
esters of vegetable oils and being their mixture with diesel oil 
should be based on research results dealing not only with energy 
and ecology aspects but also with durability and reliability 
of crucial structural units of the engines, namely: injection, 
crankshaft and piston-cylinder systems. 

In order to be able to select a kind of biofuel, out of those 
have been tested so far, additional tests on durability of biofuel- 
fed self-ignition engines should be performed. 

For application of biofuels it is necessary to achieve 
credible information concerning not only the energy aspects 
and potential pollution hazard to the environment by noxious 
substances, but also the developing of wear processes within 
the engine, i.e. its reliability and durability. In order to achieve 
a comprehensive knowledge dealing with the problems 
technical diagnostics should be used [6, 9]. 

Features of diesel oil and the considered alternative fuel 
oils influence in a different way not only energy merits of self- 
ignition engines but also their wear and operation. 

The operation of engines may be understood in different 
ways. In this paper it is understood as the process of generating, 
by the engines, the energy E within a definite time t. It has 
been equated to a physical quantity expressed by a numerical 
value measured by the unit called joule-second [J - sec]. The 
so-understood operation becomes worser and worser along 
with progressing wear of such engines. It means that value 
of the operation within a given time decreases as a result 
of lowering amount of energy generated by the engine. 
The so-interpreted operation depends on a technical state of the 
engine characterized simultaneously by the energy transformed 
by it and the energy generating time. 

The engine operation presented in this version is of such 
virtue that it can be investigated by performing measurements 
of quantities which characterize energy and can be expressed 
numerically as a value measured by the joule-second unit 
[J - sec] (Eq. 2 and 3), as well as presented in a graphical form 
as an area of operation. 

Depending on practical needs and possibile assessment 
of operation of the engines, their operation can be valuated 
as aresult of consideration of the useful work (6) which 
can be done by the engine, or the effective power (7) which 
can be developed by the engine. As engine’s power contains 
information on how fast the work can be done, hence the 
effective engine power can be considered an energy rate of its 
operation (engine operation rate). 


The so-interpreted operation, in spite of that it was 
formulated with a view of self-ignition engines, may deal also 
with spark-iginition ones. Similar interpretation of operation 
can be applied also to gas turbines and other power devices. 
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Energy losses in the hydraulic rotational motor 
definitions and relations for evaluation of the 
efficiency of motor and hydrostatic drive 


Zygmunt Paszota, Prof. 
Gdansk University of Technology 


Abstract 


The evaluation methods of energy losses and efficiency of the hydraulic rotational motors 

for the hydrostatic drives, used so far in the scientific research and in the industrial practice, 

give wrong results because the parameters that the losses and efficiencies are a function 
of are themselves dependent on those losses. 

The aim of the paper is to define the motor operating parameters, developed powers, 

energy losses and efficiencies and also to show the respective relations. Conclusions are 

drawn, based on the analyses of presented definitions and relations, on the motor energy 


investigations. 


Keywords: hydrostatic drive; hydraulic motor; energy efficiency 


INTRODUCTION 


Evaluation of the energy behaviour of a hydraulic motor 
is an evaluation of its overall efficiency ny = f (n,,, Mẹ V) 
i.e. evaluation of the overall efficiency n, as a function of 
motor shaft speed n,, and load M,, and of the working fluid 
viscosity v. This is also assessment of the value and proportions 
of the motor mechanical, volumetric and pressure losses 
deciding of the motor mechanical Ny» volumetric ny, and 
pressure n),,, efficiency, where the product ny = Nym Nyy vp 
determines the motor overall efficiency n,,. The energy losses 
and the corresponding efficiencies nyn» Ny, and Miya should be 
determined as a function of parameters having a direct impact 
on the particular losses and efficiencies. 

Designers and makers of rotational hydraulic motors and 
hydrostatic systems have not had so far a tool to determine 
their energy behaviour in the (0<04(Owm.o OSM (Munad 
field of change of the hydraulic motor shaft speed and load 
coefficients and in the v, <v < Vaa field of change of the 
working fluid viscosity. 

The rotational motor producers make erroneous routine 
evaluations of the following energy efficiencies and work 
parameters: 

— the motor overall efficiency n,, as a function of the shaft 
speed n,, and motor pressure decrease Ap,, (e.g. [11+14]) 

— the motor overall efficiency n,, as a product of the volumetric 
efficiency y,,, and the so called „mechanical — hydraulic 
efficiency” n,,,,,, determined all the three as a function of 

the same parameters (e.g. [11, 12]) 

— motor shaft speed n, as a function of the motor capacity Q,, 
and volumetric efficiency n,,,, determined in turn as a function 

of the motor pressure decrease Ap,, (e.g. [11, 13]) 
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— motor shaft torque M,, as a function of pressure decrease 
Ap,, and the so called „mechanical — hydraulic efficiency” 
Mun, OF the motor (e.g. [11+14]) 

— motor capacity Q,, as a function of the shaft speed n,, and 
volumetric efficiency Ny, determined in turn as a function 
of the motor pressure decrease Ap,, (e.g. [11+13]) 

— motor shaft usefull power P,,, as a function of the motor 
capacity Q,, and pressure decrease Ap,, and as a function 
of the motor overall efficiency n,, determined in turn as 
a function of the motor shaft speed n, and pressure decrease 


Ap,, (e.g. [11+13]). 


The hydraulic motor researchers evaluate in a wrong way 
the losses arising in the motor: 

— the motor torque M,,, of mechanical losses as a function 
of the motor pressure decrease Ap,, and shaft speed n,, 

— asum of torque M,,, of mechanical losses and the so called 
„torque of pressure losses” (resulting from the pressure 
losses Ap,, in the motor) — as a function of motor pressure 
decrease Ap,, and the shaft speed n,, 

— the intensity Q,,, of volumetric losses in the motor as 
a function of the motor pressure decrease Ap,, (or as 
a function of the motor shaft torque M,,) and as a function 
of the motor shaft speed n, 


The evaluation methods of the energy losses and efficiency 
of the rotational hydraulic motors, used so far in the scientific 
research and in the industrial practice give wrong results 
because the parameters that the losses and efficiencies are 
a function of are themselves dependent of those losses. 

There are very few informations of the motor makers 
presenting properly the motor overall efficiency ny = f (n,,, 


M,,) as a function of the motor shaft speed n,, and torque M,, 
at a specified fluid viscosity v and presenting the impact of 
viscosity v on the overall efficiency n,, (e.g. [10]). 

It is a common deficiency that no information is given 
about the dependence of the motor mechanical, volumetric 
and pressure losses on the kinematic viscosity v of the working 
fluid used in the hydrostatic drive system. 

The fundamental reason of the erroneous evaluations are 
commonly accepted views on the research methodology and 
on the method of determining the energy losses in pumps and 
in hydraulic motors . That method is based, among others, on 
the traditional reading of the energy balance of a hydrostatic 
drive system from the Sankey diagram [1+9]. The present 
unsatisfactory state is also effect of using simplified evaluations 
of the relations of particular losses to the motor or pump 
working parameters and to the working fluid viscosity. 

Therefore, the aim of this paper is to define the work 
parameters, developed powers, losses and energy efficiency 
of a rotational hydraulic motor and also demonstrating their 
complex interdependence. The analysis of those definitions 
and relations will be a basis of conclusions regarding the 
investigations of motor energy characteristics. 


ROTATIONAL HYDRAULIC MOTOR 
— WORK PARAMETERS, POWERS, 
ENERGY LOSSES, ENERGY EFFICIENCY 
— DEFINITIONS AND RELATIONS 


— Motor shaft rotational (angular) speed n,, (@,,) varies in 
the (0 < Ou (Oumas OSMy( My.) field of the hydrostatic 
drive system operation. The instantaneous n,, (@,,) value is 
required by the machine (device) driven by the motor. The 
instantaneous value of the n, (@,,) speed is independent 
of the instantaneous value of the M,, torque loading the 
motor shaft and also independent of the energy losses in 
the hydraulic motor and in the hydrostatic drive system. 

— M,, torque loading the motor shaft varies in the 
(0 < Ou (Ouna OS My£M y ,,, field of the hydrostatic drive 
system operation. The instantaneous value of M,, torque 
is required by the motor driven machine (device). The 
instantaneous value of M,, torque is independent of the 
instantaneous value of the required motor shaft speed n,, 
(@,,) and also of the energy losses in a hydraulic motor and 
in the hydrostatic drive system. 

— Working fluid (hydraulic oil, oil-water emulsion) kinematic 
viscosity v changes in the vin <v < Vna range. The 
instantaneous value v of the viscosity of fluid flux reaching 
the hydraulic motor is independent of the motor and of the 
energy losses in the motor. 

— Motor useful power P,,, required on the motor shaft by 
the driven machine (device), is a product of the required 
M,, torque loading the motor shaft and the required shaft 
angular speed œ; 


Pin = My Oy 7 2H My 1, (1) 


The motor useful power P,,, is independent of the energy 
losses in the hydraulic motor and the hydrostatic drive 
system. 

— M „torque of mechanical losses in the motor, occurring in 
the ,,shaft — working chambers” assembly, is a function of 
the required M,, torque loading the motor shaft and of the 
required shaft rotational speed n,,. The n,, speed influences 
the inertia forces of „shaft — working chambers” assembly 
elements and in effect the friction losses in the piston, 
satellite and vane motors. The M,,,, torque of losses is to 
a some extent also a function of the working fluid viscosity v. 


The impact of fluid viscosity on the mechanical losses in 
the ,,shaft — working chambers” assembly occurs mainly in 
the piston motors with fluid in the motor casing: 


Mun EA te V) (2) 


M 
Power AP, of mechanical losses in the motor, occuring 
in the ,,shaft — working chambers” assembly, is a product 
of the Mn torque of mechanical losses and shaft angular 
speed ©; 


AP j= My, 8, = 20M, 1 (3) 


Mm M 


M, indicated torque in the motor working chambers (at the 
point of conversion of the working fluid pressure energy 
into mechanical energy of the ,,shaft — working chambers” 
assembly), required by the motor from the driving working 
fluid, must be greater than the M,, torque loading the 
motor shaft [required by the driven machine (device)] 
because of the necessity of balancing also the M,, torque 
of mechanical losses in the ,,shaft — working chambers” 
assembly. The M,,, torque is equal to the sum of shaft torque 
M,, and M „n torque of mechanical losses. The indicated 
torque M,,, requires a value of the product of decrease 
Ap,,, Of the indicated pressure in working chambers and 
the theoretical motor capacity q,,, per one shaft revolution 
(theoretical motor working volume V,,,) in accordance with 
the expression: 


AP; 
Pu M = Mag =M„ +My, (4) 
The M,,, torque indicated in the motor working chambers is 
nota function of the decrease Ap, and of theoretical motor 
capacity q,,, per one shaft revolution. 
For evaluation of the M,,, torque indicated in the motor 
working chambers a formula can be used relating the 
M,, torque loading the motor shaft with the known motor 
mechanical efficiency ny, [formula (11)]: 
AP mi Ie _ _ My 
Pulm = My, = (5) 
211 Mm 


with: nym = £ Mum My) = £ M,e Mp V) 


i.e. a formula where mechanical efficiency n,,,, is defined 
as a function of parameters influencing the M,,, torque 
of mechanical losses in the ,,shaft — working chambers” 
assembly and as a function of the M,, torque loading the 
motor shaft. 

Pressure decrease Ap „indicated in the motor working 
chambers is a function of the required M,,, torque indicated 
in the chambers and theoretical capacity q,,, per one shaft 


revolution: 

AP wi = 211 Mu = 21My t Mm) (6) 

Ami Am: 

Therefore, the pressure decrease Ap, , indicated in the motor 
working chambers (with determined theoretical capacity 
qu, Per one shaft revolution) is a function of the required 
M,, torque loading the motor shaft and the M torque 
of mechanical losses in the ,,shaft — working chambers” 
assembly. The pressure decrease Ap, is indirectly a function 
of the shaft rotational speed n,, and a function of the working 
fluid viscosity v, which have an impact (apart from M,,) on 
the M „m torque of mechanical losses: 


APyi = f (Mye Mum) = f My Dy V) O 


For evaluation of the decrease Ap,,, of pressure indicated in 
the motor working chambers (with determined theoretical 


Mm? 


Mm: 
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capacity q,,, per one shaft revolution) a formula can be used 
relating the M,, torque loading the motor shaft with the 
known motor mechanical efficiency n,,,, [formula (11)]: 
211M,, 
AP yi = (8) 
qu Num 


with: nyum = f£ Mum Mi) = £ (M,e nw V) 


i.e. a formula where mechanical efficiency n,,,, is defined 
as a function of parameters influencing the M,,,, torque 
of mechanical losses in the „shaft — working chambers” 
assembly and as a function of the M,, torque loading the 
motor shaft. 

P „indicated power in the motor working chambers is 
required by the motor from the driving working fluid at the 
point of conversion of the working fluid pressure energy 
into mechanical energy of the „shaft — working chambers” 
assembly. The P „ is equal to the product of the M,,, torque 
indicated in the chambers and the shaft angular speed œ,- 
The power P,,, indicated in the working chamber is a sum of 
useful power P,, [required on the motor shaft by the driven 
machine (device)] and the power AP „m of mechanical losses 
in the „shaft — working chambers” assembly: 


AD... 
EME = Api Am: Om = Pmi = 
211 (9) 


= Mmi On = (Mu + Mum) Om = Puu + AP vim 


The P „; power indicated in the motor working chambers is 
not a function of the decrease Ap, , of pressure indicated in 
the chambers and of theoretical motor capacity q,,, per one 
shaft revolution. 
For evaluation of the P,,, power indicated in the motor 
working chambers a formula can be used relating the motor 
shaft useful power P,,, with the known motor mechanical 
efficiency Nym [formula (11)]: 

P 

Pu == (10) 

Num 

with: Nym = £ (Mum My = £ M,e Dye V) 


i.e. a formula, where mechanical efficiency n,,,, is defined 
as a function of parameters influencing the M„ torque 
of mechanical losses in the „shaft — working chambers” 
assembly and as a function of the M,, torque loading the 
motor shaft. 

The motor mechanical efficiency n,,,, is a ratio of useful 
power P,,, on the shaft [required by the motor driven 
machine (device)] to the power P,,, indicated in the motor 
working chambers (required by the motor of the driving 
fluid at the point of conversion (change) of the working 
fluid pressure energy into the mechanical energy of the 
„Shaft — working chambers” assembly). The n,,, efficiency 
can be also determined by the ratio of motor shaft torque 
M,, to the torque M,,, indicated in the working chambers: 


Mm? 


Mm’ 


fice ie Mu 
Mm Pui Piia + AP vim a 1) 
Mu Oy _ Mu Mu 
(Mu Mum) nu Mu+Mym Mmi 


The motor mechanical efficiency ny, is a function of torque 
My, Of mechanical losses in the ,,shaft— working chambers” 
assembly and the shaft loading torque M,,. Therefore, the 
Nym efficiency is a function of torque M,, and shaft rotational 
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speed n,, and a function of working fluid viscosity v, which 
influences (apart from M,,) the torque M,,,, of mechanical 
losses in the „shaft — working chambers” assembly: 


Num E fM Mi) E f (Myv Nv v) 


because: M im = f (M,,, D V). 


(12) 


Mm’ 


The motor mechanical efficiency n,,,, can be also evaluated 
from the formula: 
211 M,, 


n m = 13 
úl AP mi Um = 


However, mechanical efficiency n,,,, is not a function of the 
decrease Ap,,, of pressure indicated in the motor working 
chambers and of the motor theoretical capacity q,,, per one 
shaft revolution. 

Intensity Q „ofthe motor volumetric losses in the working 
chambers takes into account internal volumetric losses 
(between the chamber inlet channel and chamber outlet 
channel) and external volumetric losses (between chambers 
and casing and then led out of the casing). The intensity 
Q,,, of motor volumetric losses (with determined theoretical 
capacity q,,, per one shaft revolution) is a function of the 
decrease Ap,,, of pressure indicated in the chambers and, 
to some extent, of the shaft rotational speed ny as well as 
working fluid viscosity v: 


Or CAP yá ye V) (14) 
The intensity Q,,, of volumetric losses in the motor working 
chambers is a complex function of torque M,, and motor 
shaft speed n,, and also working fluid viscosity v, i.e. of 
parameters independent of the motor and motor losses. 
The decrease Ap,,, of pressure indicated in the chambers, 
influencing directly Q,,, [formula (14)], is a function 
[formula (6)] of the shaft torque M,, and of torque Mym 
of mechanical losses in the „shaft — working chambers” 
assembly and, in turn, torque M,,,, of the losses [formula 
(2)] is a function of torque M,, and of motor shaft speed ny 
and also of the working fluid viscosity v. At the same time, 
the impact of shaft speed n,, and working fluid viscosity 
v on the intensity Q,,, of the volumetric losses in working 
chambers differs from the impact ofn, and v on the torque 
M ym Of mechanical losses in the ,,shaft — working chambers” 
assembly. A direct evaluation of the relation of intensity 
Q,,, of volumetric losses in the motor working chambers 
to the motor shaft torque M,, and speed n,, and also to the 
working fluid viscosity v would be unjustified and wrong 
because it would be under a complex impact of the torque 
M uim Of mechanical losses. 
Power AP,,, of the motor volumetric losses, in the working 
chambers, is a product of decrease Ap,,, of pressure 
indicated in the chambers and intensity Q,,, of volumetric 
losses in the chambers (on the assumption that the external 
volumetric losses are small and negligible from the energy 
point of view): 


AP = AP vii Oy (15) 


Motor capacity Q,, required by the motor from the driving 
fluid must be greater than the product q,,,n,, [theoretical 
capacity q,,, per one shaft resolution and motor shaft 
rotational speed n,, required by the motor driven machine 
(device)] because of the necessity of balancing also the 
intensity Q,,, of volumetric losses in the motor working 
chambers. Capacity Q,, is equal to the sum of intensity q,,, 
n,, and intensity Q: 


Qu = dm Ou + Quv (16) 


= Mu Oy + Mun Ou + 


Evaluation of the motor capacity Q,, (with determined 
theoretical capacity q,,, per one shaft revolution) can 
be performed with a formula including the motor shaft 
rotational speed n,, required by the motor driven machine 
and known motor volumetric efficiency ny, [formula 

Q3)1: n 
Qu = Amt M (17) 

Nv 


with: ny, = f (Quo Dy) = f (APs Bap V) 


i.e. a formula, where the volumetric efficiency ny, is 
determined as a function of parameters influencing the 
intensity Q,,, of volumetric losses in the working chambers 
and as a function of the motor rotational speed n,,- 

Power P , of the working fluid absorbed by the motor 
in working chambers is required by the motor from the 
driving fluid as a difference between power p,,,, Qu of the 
fluid inflowing in the working chambers from the inlet 
channel and power p,,,; Qu of the fluid outflowing from 
the chambers to the outlet channel. Assuming that the 
external volumetric losses are small and negligible from the 
energy point of view, it may be accepted that the intensity 
of the outflowing flux is equal to the intensity Q,, of the 
inflowing flux. Therefore, power P, may be determined 
as a product of the decrease Ap,,, of pressure indicated in 
working chambers and the motor capacity Q,,. Power P yei 
must be greater than power P „ indicated in the chambers 
(required by the motor from the driving fluid at the point of 
conversion (change) of pressure energy of working fluid into 
the mechanical energy of the „shaft — working chambers” 
assembly) because of the necessity of balancing also the 
power AP „ of volumetric losses in the chambers. Power 
P ua; 18 equal to the sum of power P , and power AP.,, : 


are Pm Qu — Pumzi Q= 
= AP O AP Gaeta Om) = 08) 
5 APyi Am Pu Oi ie Pui t AP w 


Power P ,., of the working fluid consumed by the motor in 
working chambers is a sum of useful power P,,, [required 
on the motor shaft by driven machine (device)], power AP ,, 
of mechanical losses in the „shaft — working chambers” 
assembly and power AP „ of volumetric losses in the motor 
working chambers: 


Pia = E + AP im + AP (19) 


After replacing in the equation (19) the useful power P,,,, 
power AP,, of mechanical losses and power AP,, of 
volumetric losses by the expressions relating those powers 
to the parameters and losses deciding of their values, 
a picture can be obtained of the impact of those parameters 
and losses on power P,,,, consumed in the working 
chambers: 


Puci = Mu Om + Mum Om + AP oi Qu, = 
211 (My +Myim) 
dm 

Qu 


= 20 (Myu +M yum fo, + Qs 
Am: 


Power Pu; of working fluid consumed by the motor in 
working chambers can be evaluated by means of a formula 


Qu, = (20) 


expressing the ratio of power P,,, indicated in the working 
chambers to a known volumetric efficiency n,,, of the motor 
[formula (23)]: 
Pp. 
A (21) 
Mv 


with: hig f (Qie Ny) = f (AP yi Ny? v) 


i.e. a formula where volumetric efficiency n,,, is determined 
as a function of parameters influencing the intensity Q,,, of 
volumetric losses in the working chambers and as a function 
of the motor rotational speed n,,. 

Power P y; of the working fluid consumed by motor in the 
working chambers can be evaluated also from the known 
useful power P,,, on the motor shaft, known mechanical 
efficiency Nym [formula (11)] and known volumetric 
efficiency n,,, of the motor [formula (23)]: 


Pu 
me = a (22) 
Mmm Nv 
with: nym = f Mya» My) = £ M,e Mp V) 


and: Ny = f (Qu My) =f (Ap, Nhs v) 


i.e. a formula where mechanical efficiency n, is determined 
as a function of parameters influencing the torque M ym of 
mechanical losses in the ,,shaft — working chambers” 
assembly and as a function of torque M,, loading the motor 
shaft, and the volumetric efficiency Ny, is determined as 
a function of parameters influencing the intensity Q,, of 
volumetric losses in the working chambers and of the motor 
rotational speed n, 

Volumetric efficiency n,,, of the motor is a ratio of power 
P,,, indicated in the motor working chambers to power P yei 
of working fluid consumed by the motor in the chambers: 


Pui Pu 
Tv = = —___—_ = 
Puci Pui + AP 
= AP Mi Ami Am = (23) 
AP ui Am: Am + AP i Qu, 
qm My Ami AM 


The motor volumetric efficiency n,,, (with determined 
theoretical capacity q„ per one shaft revolution) is a function 
of intensity Q,,, of volumetric losses in the motor and of the 
motor shaft rotational speed n,,. Therefore, efficiency ny, 
is a function of the decrease Ap,,, of pressure indicated in 
working chambers and a function of the rotational speed n,, 
as well as a function of the working fluid viscosity v (which 
have an impact on the intensity Q,,, of volumetric losses) 
and also directly a function of rotational speed n: 


Muy = f (Quv My) =F (AP Be VY) QA) 
because: Qu, =f (APui My» V). 


Losses Ap, of working fluid pressure in the motor channels 
are a sum of two pressure losses i.e. loss Ap of pressure 
in the inlet channel (between the motor inlet point and 
working chambers) and loss Ap,,,, of pressure in the outlet 
channel (between the working chambers and the motor 
outlet point). Losses Ap,,, are a function of motor capacity 
Q, and of working fluid viscosity v: 


AP sip = AP spt + AP tps (Qs V) 05) 


Losses Ap, y, of working fluid pressure in the motor channels 
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are a complex function of the motor shaft speed n,, and 
torque M „and the working fluid viscosity v, i.e. parameters 
independent of the motor and of losses in it. The motor 
capacity Q,,, having a direct impact on AP up [formula (14)], 
is a function [formula (16)] of the shaft rotational speed n, 
and of intensity of volumetric losses Q,,, in the working 
chambers. The decrease Ap,,, of pressure indicated in the 
working chambers, having a direct impact on Q,,, [formula 
(14)], is a function [formula (6)] of the shaft torque M,, and 
of torque M ym of mechanical losses in the „shaft — working 
chambers” assembly, and the torque M,,,, of mechanical 
losses [formula (2)] is in turn a function of the motor shaft 
torque M,, and speed n,, and of the viscosity v of working 
fluid. At the same time, the impact of working fluid viscosity 
v on the losses Ap, of fluid pressure in the channels differs 
from the impact of viscosity v on the intensity Q,,, of 
volumetric losses in the working chambers and from the 
impact of v on the torque M,,, of mechanical losses in the 
„Shaft — working chambers” assembly. Also the impact of 
shaft speed n,, on the intensity Q,,, of volumetric losses in 
the working chambers differs from the impact of n,, on the 
torque M,, of mechanical losses in the ,,shaft — working 
chambers” assembly. Direct evaluation of the dependence 
of pressure losses AP up of the working fluid in the motor 
channels on the motor shaft speed n,, and torque M,, and 
on the viscosity v of working fluid would be unjustified and 
wrong, because it would be under a complex impact of the 
intensity Q,,, of volumetric losses in working chambers and 
of torque M „m of mechanical losses in the „shaft — working 
chambers” assembly. 

Power AP, of pressure losses in the motor, in the motor 
channels, with the assumption that external volumetric 
losses are small and negligible from the energy point of 
view, is a product of the pressure losses AP tp in the channels 
and the motor capacity Q: 


AP Api Qu (26) 


Decrease Ap,, of pressure in the motor (with determined 
theoretical capacity q,,, per one shaft revolution), required 
by the motor from the driving working fluid, must be greater 
than decrease Ap,,, of pressure indicated in the working 
chambers (required by the torque M,,, indicated in the 
chambers) because of the necessity of balancing also the 
losses AP of pressure in the motor channels. Decrease 
Ap,, is equal to the sum of the indicated decrease Ap,,, and 
losses AP vip: 


APy = AP ui + AP ap (27) 


Replacing in the equation (27) the decrease Ap,,, of pressure 
indicated in the working chambers with expression (6), we 
obtain the dependence of decrease Ap,, of pressure in the 
motor on the required torque M,, loading the motor shaft 
and on the torque M,,,, of mechanical losses in the „shaft 
— working chambers” assembly and also on the pressure 
losses Ap,,, in the motor channels: 


211 (My + Mym) 
qm 


The decrease Ap, of pressure in the motor can be evaluated 
by means ofa formula expressing the ratio of decrease Ap, y; 
of pressure indicated in the working chambers to the known 
pressure efficiency gs of the motor [formula (37)]: 


AP = + AP mip (28) 


AD,,. 
Mp 
with: ie f (APY AP yi) =f (Qu APye Y) 
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i.e. formula where the motor pressure efficiency nyp 1 
defined as a function of parameters influencing the oe 
AD i of working fluid pressure in the channels and as 
a function of the decrease Ap, ,, of pressure indicated in the 
working chambers. 

Decrease Ap,, of pressure in the motor can be evaluated 
also from a known torque M,, loading the motor shaft, 
from a known mechanical efficiency n,,,, [formula (11)] and 
a known motor pressure efficiency Mai [formula (37)]: 


211M 
Api a (30) 
qm Num Tp 
with: Nym = £ Mym My = £ (M,e Dap V) 
and: Thy, = f (AP Apa) = f (Qu Pup V) 
i.e. formula, where the mechanical efficiency n,,,, is defined 


as a function of parameters influencing the torque Mym 
of mechanical losses in the ,,shaft — working chambers” 
assembly and a function of torque M,, loading the motor 
shaft and the pressure efficiency n,,, is defined as a function 
of parameters influencing the losses Ap, , pOf pressure in the 
channels and as a function of decrease Ap,,, of pressure 
indicated in the motor working chambers. 

Working fluid power P „consumed by the motor must 
be greater than power P; (consumed by the motor 
in the working chambers) because of the necessity of 
balancing also the power AP,, of pressure losses in the 
motor channels. Power P, is equal to the sum of power 
P ya; and power AP,,, of the losses. Power P, is a product 
of the decrease Ap,, of pressure in the motor and motor 
capacity Q,;: 


Pir = APy Qu 
Bro Qu 


The working fluid power P,,, consumed by the motor is 
a sum of useful power P,,, (required on the motor shaft by 
the driven machine (device)], power AP,, of mechanical 
losses in the ,,shaft — working chambers” assembly, power 
AP,,, of volumetric losses in the working chambers and 
power AP w of pressure losses in the motor channels: 


P,.=P,, FAP y PAP CHAP, 0) 


After replacing in equation (32) the useful power Pw. 
and power AP „m of mechanical losses, power AP,,, of 
volumetric losses and power AP,,, of pressure losses ‘with 
formulae expressing dependence on parameters and losses 
deciding of the values of those powers, we can obtain 
a picture of the impact of parameters and losses on the 
consumed power P, ,. 


= (AP y+ Ay.) Qu = 
T AP vip Qu = Puci T AP p 


(31) 


me = My Oy + Mym On + APmi Qw + AP op Qu = 


= My Oy +Myim Oy + 21M y +Myng + 
Amt 
+ AP vip Gam +Qu)= (33) 


= 2T1(My, + Main Yn, qe Que), 
Ami 


+ AP op Gar Mu + Ola) 


The expression describing the working fluid power P. 
consumed by the motor can be also obtained from the 


product of decrease Ap,, of pressure in the motor [formula 
(28)] and motor capacity Q,, [formula (16)]: 


Pue = ApyQu= 


_ 2I (My Msi) 
Ame 


Expressions (33) and (34) are equivalent. 

Evaluation of the working fluid power P,,, consumed by 
the motor can be performed with a formula expressing the 
ratio of working fluid power P „„ consumed by the motor 
in the working chambers to the motor pressure efficiency 
Np [formula (37)]: 


(34) 
+ Apu (Am My + Quy) 


Pus = (35) 
j Nimp 


with: yy, = f (AP APs) = f (Que APy> V) 


i.e. a formula where the pressure efficiency n,,, is defined 
as a function of parameters influencing the losses Ap,,, of 
working fluid pressure in the channels and as a function of 
decrease Ap,,, of pressure indicated in the motor working 
chambers. 

The working fluid power P ,, consumed by the motor can be 
evaluated from the known useful power P,,, on the motor 
shaft, a known mechanical efficiency Nym [formula (11), 
known volumetric efficiency n,,, [formula (23)] and known 
pressure efficiency Ny, [formula (37)] of the motor: 


P 
Pu, = —— (36) 
Num Tv mp 
with: My, = f Mum Mw f (M,, ty» V) 
Muy z F (Quv su a f(Ap,,. 0 M?’ v) 
and: My, = f (APup AP, )= f (Que APys V)- 


In formula (36), the mechanical efficiency n,,,, is defined 
as a function of parameters influencing the torque M,,,, 
of mechanical losses in the „shaft — working chambers” 
assembly and as a function of torque M,, loading the motor 
shaft. The volumetric efficiency n,,, is defined as a function 
of parameters influencing the intensity Q,,, of volumetric 
losses in the working chambers and as a function of the 
motor shaft rotational speed n,,. The pressure efficiency 
Thi is defined as a function of parameters influencing the 
pressure losses AD of working fluid in the channels and 
as a function of the decrease Ap,,, of pressure indicated in 
the motor working chambers. 

The motor pressure efficiency n,,, is a ratio of the working 
fluid power P; consumed by the motor in the working 
chambers to the power P „, consumed by the motor: 


-Pwa Pua _ 
TM Phe Pages + APA 
= AP ui Qu z Apat. _ AP ui 2 
(APui + APup) Qm APmi +APMp  APy 


Therefore, the pressure efficiency can by presented as a ratio 
of the decrease Ap, of pressure indicated in the working 
chambers to the decrease Ap,, of pressure in the motor. The 
motor pressure efficiency Ny, is a function of losses AP, of 
the working fluid pressure in motor channel and decrease 
Ap,,, of pressure indicated in the motor working chambers. 
Therefore, efficiency ny »is a function of the motor capacity 
Q,, and a function of the working fluid viscosity v (which 
influence the losses AP up of working fluid pressure in the 


channels) as well as a function of the decrease Ap,,, of 
pressure in the motor working chambers: 


Nop = f APuip> Api =f (Qu APyi V) 
because: AP = £(Q , V). 


(38) 


— The motor overall efficiency n, is a ratio of the useful power 
P,,, on the motor shaft required by the driven machine 
(device) to the power P ,, consumed by the motor: 


_ Pm Mic _ 2M Mu 1m 
Pue  APm Qu Apy Qu 


Replacing in formula (39) the power P,,, consumed by 
the motor with equations describing its dependence on 
the useful power P,,, and on the powers AP ym» AP,, and 
AP „ of the energy losses i in the motor (with a a “determined 
theoretical capacity q Per one shaft revolution), we obtain 
the expressions describing the motor overall efficiency ny 
as a function of the losses: 


(39) 


with reference to equation (32): 


= Piu Z Piu 40 
k Pue Pra + AP yim + AP aay + AP wip ' i 
with reference to equation (33): 
2I My ny, (41) 
Ny = Q 
211 (My + Muns + J+ AP Mp (Aum Am + Quy) 
Mt 
with reference to equation (34): 
ae 2I My ny, (42) 
e 
211 +M 
ES $ 20. (Am Mx + Quy ) 
Mt 


Expressions (41) and (42) are equivalent. 


The motor overall efficiency n,, is therefore a function 
of torque Mym of mechanical losses in the „shaft — working 
chambers” assembly, intensity Q,,, of volumetric losses in the 

working chambers and losses AD vy of pressure in the motor 
channels. The nų efficiency is also a function of the motor 
shaft torque M,, and speed n,, required by the driven machine 
(device): 


M = f (Mim Quv AP ip? My ny) 


Torque M,,,, = £ (M,,, n,, V) of mechanical losses from 
friction of elements in the motor „shaft — working chambers” 
assembly is a function of the required torque M,, loading the 
motor shaft. In the piston, satellite and vane motors, torque 
Mum Of the losses is also a function of the required speed n,, 
influencing the inertia forces in the „shaft — working chambers” 
assembly and in effect the forces of friction between those 
elements. In the piston motors in particular, with working fluid 
in the casing, the torque M,,_ of the losses is also a function of 
the fluid viscosity v, which influences the friction between the 
„Shaft — working chambers” assembly elements and the fluid. 

Intensity Qu, = f (Ap,,, Dw V) of volumetric losses in the 
motor working chambers is a function of the decrease Ap,,, 
of pressure indicated in the chambers and, to some extent, 
a function of motor shaft rotational speed n,, and also a function 
of working fluid viscosity v. 

Losses Ap,,, = f (Qp V) of the working fluid pressure in the 
motor channels are a function of the motor capacity Q, and of 
the working fluid viscosity v. 


(43) 
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In order to evaluate the dependence of overall efficiency n,, 
of the motor (with determined theoretical capacity q,,, per one 
shaft revolution) on parameters independent of the losses in the 
motor, i.e. evaluate n,, as a function of the required motor shaft 
torque M „and speed n,, and also as a function of the working 
fluid viscosity v, a product of the motor mechanical efficiency 
Nm Volumetric efficiency ny, and pressure efficiency Ny, 
must be used: 


M =f (My, Dy v) = Mim Muy Mo (44) 


with: Nym 7 £ Mum My = £ M,e Dye V) 


because: My, =f (Mo, np V) [equation (12)] 
Muy = f (Que ny) = f (AP > Dyp v) 
because: po = f (Ap,,,. Dw V) [equation (24)] 
and: = f (APup APy) = f (Quy APup V) 
because e =f (Q,, v) [equation (38)]. 


In the above equation describing the overall efficiency 
Mm» Mechanical efficiency Ny, is defined as a function of 
parameters influencing the torque M,,, of mechanical losses in 
the „shaft — working chambers” assembly and as a function of 
torque M,, loading the motor shaft. Volumetric efficiency ny, 
is defined as a function of parameters influencing the intensity 
Qu, of volumetric losses in the working chambers and also 
as a function of the motor shaft rotational speed n,,. Pressure 
efficiency ną is defined as a function of parameters influencing 
the losses AD, of working fluid pressure in the channels and 
a function of decrease Ap,,, of pressure indicated in the motor 
working chambers. 

After replacing in equation (44) n,,, with expression (13), 
Nyy With expression (23) and nyp with expression (37), we 
obtain a formula describing the motor efficiency ną as a ratio 
of useful power P „ to power P „, consumed by the motor, i.e. 
a formula confirming correctness ‘of the expressions describing 


Num Muy and Np: 


MT f(Mu: nu: V) = Nimm vv Mp = 
_ 2TIMy qmm AP — 20M nm _ Pu a 
Apidae Qu APm  QOuAPm Pre 


In equation (45), the required motor capacity Q,, [equation 
(17)] is a function: 


Que AmMy 
Nv 
with: My, = £ (Que My) = £ (AP, yo My V) 
because: Qu, =f (Ap,,. My V) 


and the required pressure decrease Ap,, (equation (30)] is 
a function: 


2IIM 
APy = ~ 
qm N mm Amp 
with: Nim = L (My, My) = £(M,,, nw V) 


because: Mym = f (Mpp yy y) 


and: Mira 7 £ (APup AP) = f (Que APup Y) 
because: AP vip = f (Que V). 


After replacing in equation (44) n,,,, with formula (11), 
Nyy With formula (23) and Tias with formula (37), we obtain 
an expression describing the efficiency n,, as a product of 
individual efficiencies described by losses and parameters 
deciding of their values and where at the same time Ap, and 
Q,, are functions of the losses. 
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Nu = f (Mu: nms V) = Nim viv ap = 
- Mu dM APs 
M m +Mym qm My + Quy AP i + AP vip 


(46) 


where: Mum = f (M,e n,,, V) [equation (2)] 
Qus = f (AP Nw V) [equation (14)] 


AP ip T f (Qu v) [equation (25)] 
Wi: Aas My t vn) [equation (6)] 
Mt 
and: Qu = am By + Qu, [equation (16)]. 


Decrease Ap,,, of pressure indicated in the working chambers 
[equation (6)] is a function of the loading torque M „and torque 
M uim of mechanical losses in the „shaft — working chambers” 
assembly. Capacity Q,, in the motor channels [equation (16)] 
is a function of motor shaft speed n,, and the intensity Q,,, of 
volumetric losses in the working chambers. 

Formula (46) shows a direct dependence of the torque 
Mum of mechanical losses in the „shaft — working chambers” 
assembly on the torque M,, and on the motor shaft rotational 
speed n,, as well as on the working fluid viscosity v. 

Formula (46) presents a complex dependence of the intensity 
Qu, of volumetric losses in the working chambers on the shaft 
loading torque M,, and on the torque M „m of mechanical losses 
in the „shaft — working chambers” assembly (decrease Ap,,, of 
pressure indicated in the working chambers depends on M,, 
and M,,,, and has direct impact on Q,,,) and also on the shaft 
speed n,, (influencing in diversified way the torque Mym of 
mechanical losses and intensity Q,,, of volumetric losses). The 
intensity Q,,, of volumetric losses depends on diversified impact 
of the working fluid viscosity v: indirectly by impact of v on 
the torque M,,,, of mechanical losses in the „shaft — working 
chambers” assembly and directly by impact of v on intensity 
Q,,, Of losses in the chambers. 

Formula (46) presents also a complex dependence of losses 
Ap,,, Of working fluid pressure in the channels on the shaft 
rotational speed n,, and on intensity Q,,, of volumetric losses 
in the working chambers. The intensity Q,,, of losses influences 
the motor capacity Q,, [equation (16)] and at the same time Q,,, 
depends in a complex way on the shaft loading torque M,, and 
on the torque M,,, of mechanical losses in the ,,shaft— working 
chambers” assembly. Pressure losses Ap,,, in the channel are 
also dependent on the diversified impact of the working fluid 
viscosity v: indirectly by impact of v on the torque M,,,, of 
mechanical losses in the „shaft — working chambers” assembly 
and by impact of v on the intensity Q,,, of volumetric losses in 
the working chambers and directly by impact of v on the losses 
Apy, Of pressure in the channels. 


ANALYSIS OF THE PRESENTED 
DEFINITIONS AND RELATIONS 


1. The power P,,, consumed by the motor is a sum of motor 
shaft useful power P,,, and powers of three different energy 
losses in the motor. The losses occur in series increasing 
power stream in the opposite direction to the direction 
of power flow. In effect, the power stream in the motor 
increases from the shaft useful power P,,, to the working 
fluid power P,,, consumed by the motor: 


Pue = Ps ! AP ym ! AP w ! AP up 
Mechanical losses (and power AP,,_) occur in the „shaft 
— working chambers” assembly, volumetric losses (and 


power AP, , ) occur in the working chambers, pressure losses 
(and power AB) occur in the motor channels. 


. Figure 1 presents a diagram of the direction of increasing 
power stream in a hydraulic motor. Direction of the increase 
of power stream is opposite to the direction of power flow 
in the motor. The diagram replaces the Sankey diagram 
of distribution of power flowing in a power transmission 
system. The use of Sankey diagram for description of 
power stream in the power transmission systems is a basic 
cause of errors in evaluation of losses in the power flow. 
The Sankey diagram suggests determination of losses in 
a system as a function of input parameters of that system. 
The suggestion can be noticed in the method of hydraulic 
motor investigations and in the related evaluations of losses 
and the motor energy efficiency. But the input parameters 
of the system depend on the losses. 
. Torque M,,,, of mechanical losses in the „shaft — working 
chambers” assembly should be evaluated as a function M, ,,, 
=f(M,,, 2w V), Le. as a function of the motor shaft loading 
torque M,, and the shaft rotational speed n,, and also as 
a function of a working fluid viscosity v. 
The picture of torque M,, of mechanical losses in the 
motor, presented in literature and in the industrial practice as 
a function M um = f(Ap,,, Ny, V), 1.e. as a direct dependence 
on the decrease Ap,, of pressure in the motor, is incorrect 
because it bears the impact of the mechanical losses and also 
of volumetric losses in the working chambers and pressure 
losses in the motor channels. 
. Intensity Q,,, of volumetric losses in the motor working 
chambers should be evaluated as a function Q,,, = f (AP 
np V), 1.e. a function of decrease Ap, , of pressure indicated 
in the working chambers and of motor shaft rotational speed 
n,, and also as a function of the working fluid viscosity v. 
The picture of intensity Q,,, of volumetric losses, presented 
in literature and in industrial practice as a function Q, = 
f(Ap,,, My V), 1.e. as a direct dependence on the decrease 
Ap,, of pressure in the motor, is incorrect, because it bears 
the impact of the pressure losses in the motor channels. 
Similarly, the picture of intensity Q,,, of volumetric losses 
as a function Q,,, =f(M,,, My» V); i.e. as a direct dependence 
on the motor shaft loading torque M,,, is also incorrect, 
because it bears the impact of the mechanical losses in the 
„shaft — working chambers” assembly. 
. Pressure losses Ap,,, in the motor channels should be 
evaluated as a function APwp = £(Q,, V), i. e. a function of 
motor capacity Q,, and of working fluid viscosity v. 
The picture of pressure losses Ap,,,, presented sometimes 
in literature and in the industrial practice as a function Ap, , 
= f(n,,, V), Le. as a direct dependence on the motor shaft 
rotational speed n,,, is incorrect because it bears the impact 
of mechanical losses in the „shaft — working chambers” 
assembly and the impact of volumetric losses in the motor 
working chambers. 
. The torque M,, of mechanical losses in the ,, shaft 
— working chambers” assembly and the so called ,,torque” 
of pressure losses AD in the motor channels cannot make 
up a „sum” and also that „sum” cannot be evaluated as 
directly dependent on the same chosen parameters (which is 
practiced in literature and in industry), because those losses 
are of different character and are dependent on different 
parameters [M um = f(M,, n,, V) AD. = f (Q Y]. 
. The impact of working fluid viscosity v on: 
— torque Mym of mechanical losses in the ,,shaft — working 
chambers” motor assembly 
— intensity Q,,, of volumetric losses in the motor working 
chambers 
— pressure losses Ap,,, in the motor channels 
is differentiated. 


The dependence of particular kinds of losses on the 
working fluid viscosity v should be presented in expressions 
describing dependence of those losses on other parameters 
which have a direct impact on them [M n = f(M,,, ny, V), 
Qu, = FAP ye Dye V), APiyp = f (Que VI 

The overall efficiency of the motor (with determined 
theoretical capacity q,,, per one shaft revolution), in the 
(0<04( Ox,» 0<My(My,. ) range of hydraulic motor 
shaft speed and torque coefficients and in the v in SV S Vna 
range of working fluid viscosity, must be evaluated only 
as a function ny = f (M,,, n,,, V)» ie. as a function of the 
required motor shaft torque M,,, required shaft rotational 
speed n,, and as a function of the working fluid viscosity 
v. Torque M,, and speed n,, are parameters required by the 
motor driven machine (device), independent of the motor 
and of the motor losses. The working fluid viscosity v is 
also independent of the motor and of the motor losses. At 
the same time those parameters (M,,, n,,, V), in a direct or 
indirect way, have an impact on the motor mechanical, 
volumetric and pressure losses and also on the internal 
parameters deciding of the losses: on pressure decrease 
Ap, indicated in the working chambers and deciding of the 
capacity Q,,, of volumetric losses in the working chambers 
as well as on the motor capacity Q,, directly deciding of the 
pressure losses AP uy in the motor channels. 

The motor overall efficiency n,, as a function of the 
motor shaft torque M,, and speed n,, and as a function of 
the working fluid viscosity v, is a product of the motor 
mechanical efficiency Ny» volumetric efficiency ny, and 
pressure efficiency ht 


P 
Mu = £(M,, typ V) = “a = Nym Uv np 


Mc 
Each of the three efficiencies, as a factor of the product 
describing the overall efficiency, is evaluated as a function 
of parameters having a direct impact on the respective losses 
and as a function of a parameter to which those losses are 
added”. 


10. The mechanical efficiency: 


11. 


P 
=M= My =f (Munum Y) 
Pui My + Mm 


must be evaluated as a function of parameters which 
have a direct impact on the torque Mym = f (M,,, Mp, V) 
of mechanical losses in the „shaft — working chambers” 
assembly, i.e. a function of the required motor shaft torque 
M,, and a function of the required shaft rotational speed 
ny as well as a function of the working fluid viscosity v. 
At the same time, the mechanical efficiency ny, is directly 
a function of the shaft torque M,, because the torque 
Mum Of mechanical losses is ,added” to the torque M,» 
causing decrease of power transmission efficiency in the 
assembly. 


The volumetric efficiency: 


Num 


Pus qmm 
Mm = = oo Ef An) 
= Puci Umm + Quy iia 


must be evaluated as a function of parameters which 
have a direct impact on the intensity Q,,, of volumetric 
losses in the working chambers, i.e. a function of pressure 
decrease Ap,,, indicated in the chambers and as a function 
of the required motor shaft rotational speed n,, as well as 
a function of the working fluid viscosity v. At the same 
time, the volumetric efficiency ną, is directly a function 
of the shaft rotational speed n,,, because the intensity 
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R= Pumu; viscosity » of the working fluid 


direction of power flow 


inlet channel (Apyo1) 


= Far Pua =(Pui Puz) WAPu A R uci tARip 


"shaft — working 
_chambers' assembly 


Fuci=Apui Qu = 


shaft speed ny (cy) 
shaft torque My 


FA PQs t AR on 


direction of power flow 


\working chambers 


a e 


outlet channel (Aprpa) / Q, =p! 


direction of power flow 


Fuz=Pw2Qyi outlet pressure py,*0 


Fig. 1. Diagram of the direction of increasing power stream in a rotational hydraulic motor; direction of the increase of power stream is opposite to the 


direction of power flow in the motor. Power stream increases from the motor useful power P „ required on the motor shaft by the driven machine (device) 


to power P „ consumed and required by the motor from the working fluid. The increase of power stream is an effect of the power of losses in the motor: 


power AP „ of mechanical losses in the ,,shaft— working chambers” assembly, power AP,,, of volumetric losses in the working chambers and power AP 


of pressure losses in the motor channels. Powers AP,,,, AP,,, and AP „ of the losses are functions of the output parameters of the motor assembly where the 


losses occur and diversified functions of the working fluid viscosity v: power AP,,, of mechanical losses is a function of torque M,, and shaft speed n,, (w,) 


required of the motor by the driven machine (device) and a function of the working fluid viscosity v, power AP „ of volumetric losses is a function of the 
decrease Ap,, of pressure indicated in working chambers (torque M,,, indicated in the chambers) and of the shaft rotational speed n,, as well as a function 
of the working fluid viscosity v, power AP ti of pressure losses is a function of motor capacity Q „and of the working fluid viscosity v. "Power P „ indicated in 
the working chambers: P = Py, + AP p» power P „; Of the working fluid consumed in the working chambers: P ga = Pip t AP m + AP vip power P u 9f the 
working fluid consumed by the motor: Py, = Py, + AP y, + AP ig, 7 AP y T he diagram replaces the Sankey diagram of distribution of power in transmission 
systems, causing incorrect loss evaluation during the hydraulic motor energy investigations. 


Qu, Of volumetric losses is „added” to the product of 12. The pressure efficiency: 


theoretical capacity q,,, per one shaft revolution and speed E Puci _ APyi -f(A 

n,,, causing decrease of power transmission efficiency in Amp = P. A A =f (APui Quo V) 

the chambers. Me Pmi + AP op 

If we wish to present the motor volumetric efficiency ny, must be evaluated as a function of parameters which have 

as a factor in the product Nym Muy Mp describing the motor a direct impact on the pressure losses Ap,,, in the motor 

overall efficiency n,,, i.e. to present Nyy aS a complex channels, i.e. as a function of the motor capacity Q,, and 

dependence on the parameters (M,,, Nw V) describing the a function of the working fluid viscosity v. At the same 

overall efficiency n,, and as a dependence on the mechanical time, the pressure efficiency Tmp is directly a function 

losses in the motor, the intensity Qu, = f (AP, y» My V) of of the pressure decrease Ap,,, indicated in the motor 

volumetric losses in the chambers should be determined working chambers, because pressure losses AP up in the 

with: motor channels are ,,added” to the pressure decrease AP ue 
211 My En Mum) causing decrease of power transmission efficiency in the 

BD ai E Aá>>-->= A channels. 

Amr If we wish to present the motor pressure efficiency ny, as 
and with torque M,,,, of mechanical losses in the „shaft a factor in the product Thhtm Ty Tp describing the motor 
— working chambers" assembly as a function of M,, =f overall efficiency n,, i.e. to present Nyu, aS a complex 
(My Mp Y): dependence on the parameters (Mp 2y v) describing the 
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13. 


overall efficiency n,, and as a dependence on the mechanical 
and volumetric losses in the motor, the pressure losses 
AP yp = £ (Qu, V) in the channels should be determined with 
Qu = am By + Que the intensity Qu, = f (Ap,,,, ny, V) of 
volumetric losses in the chambers should be determined 
with: 


_ 211 (Mu + Mum) 
' dm 


and with torque M,, of mechanical losses in the „shaft 
— working chambers” assembly determined as a function 
of M,,,, = f (Mp My, V)- 

Therefore, the picture of the characteristics of overall 
efficiency as a product nym Ny, Tiy ofefficiencies correctly 
described by the characteristics of mechanical efficiency 
Ny 7 £ (M,e n,,, V), volumetric efficiency ny, = (AP, y; Dw V) 
and pressure efficiency he = £(AP y; Qu V) is complex. 


CONCLUSIONS 


APm 


The methods of investigation of the rotational hydraulic 

motor losses and energy efficiency, used in the scientific 

research and in industrial practice, give incorrect results 
because: 

— losses and efficiencies are evaluated as functions of 
parameters which depend on those losses or which have 
no direct impact on the losses, 

— the mechanical, volumetric and pressure losses and the 
corresponding efficiencies are presented as directly 
dependent on the same parameters, although each of 
those losses is a function of different parameters and is 
a different function of the working fluid viscosity v. 


In the investigations of the hydraulic motor (pump and 
a hydrostatic transmission system) losses and energy 
efficiency it is necessary to use as a guide the diagram of 
the direction of increasing power stream from the hydraulic 
motor shaft to the pump shaft. 


The complex method of evaluation of the motor overall 
efficiency ny = f (M,,, M,, V) as a product Nym Nyy, eee of 
three efficiencies correctly described by the characteristics 
of mechanical efficiency Nym = f (Mp n,,, V), volumetric 
efficiency ny, =F(AP, y; Dyp V) and pressure efficiency ny, = 
f(Ap,,,, Que V) should be replaced by a method of evaluation 
of the motor efficiency based on the defined coefficients k, 
of the motor and the motor driving system energy losses. 
The proposed motor efficiency evaluation is performed as 
a part of the energy efficiency evaluation of the hydrostatic 
driving system where the motor is used. 


The evaluation method of the hydraulic motor (and also 
of the pump and of the hydrostatic driving system) energy 
efficiency is based on the mathematical models of losses 
where each type of losses is a function of parameters 
influencing directly the losses and independent of those 
losses. 

Evaluated are the loss coefficients k, relating the hydraulic 
motor (pump and system) mechanical, volumetric and 
pressure losses to the reference values: driving system 
nominal pressure p,, pump theoretical capacity Q,,, motor 
theoretical rotational speed n „and theoretical torque M „e 
The loss coefficient k, are determined at the working fluid 
reference viscosity v . Also the impact is determined of the 


viscosity ratio v/v, (viscosity changing in the V in SY S Vn 
range) on the value of loss coefficients k.. 

The method allows to evaluate the values and proportions 
of mechanical, volumetric and pressure losses in the motor 
(pump, driving system) and their dependence on the fluid 
viscosity v. 

The knowledge of coefficients k, of the mechanical, 
volumetric and pressure losses allows to obtain, by applying 
a numerical method, a picture of the overall efficiency 
Ny =f (@,,M,,) of the motor (pump and system) in the 
(0 < Oy (Oy, > 0<Myul Muna) motor operating field and 
for the selected ratio v/v, of fluid viscosity. 
Simultaneously, the (0 < Oy (Ommas 0 My (My, ) motor 
(pump and system) operating field is determined for 
a selected v/v, ratio of the working fluid viscosity to the 
reference viscosity. 

It is assumed that the method is precise and simple in use. 
It reduces the necessary laboratory investigations of pumps 
and hydraulic motors. Il allows to seek energy saving 
displacement machine designs. It also allows to evaluate the 
drive energy efficiency and to seek energy saving structures 
of hydrostatic transmission systems. 


ax 
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ABSTRACT 


The article assesses the applicability of new refrigerants, being possible substitution for Freon 22, in marine 

refrigeration and air-conditioning systems. A collection of technical data on the physical properties of the 

new refrigerants is presented. Suggestions are made about certain energetic and ecological profits which 

can be gained when using energetically efficient refrigerants with low GWP values. The results of the TEWI 
comparison analysis for selected refrigerants are also included. 
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INTRODUCTION 


The applicability of the refrigerants R41, R152a, R218, 
R227ea, RC318, R600a, R290, and R1270 in cooling 
installations was analysed based on the assessment of the 
use of these substances in the heat pump cycle, for devices 
operating in the heating mode. The medium used in the 
reference heat pump installation was R22, due to its favourable 
thermodynamic properties, along with low energy consumption 
of the devices filled with this medium. Indeed, many substances 
which are used as substitutes for Freon 22 reveal higher energy 
consumption and GWP. The main criterion used for assessing 
the applicability of the substances, referred to as refrigerants 
— refrigerating media, is their influence on the environment. 
The Carnot cycle is considered here the reference for real cycles 
realised using the media being homogeneous substances, such 
as Freon R22 and its substitutes. 


ECOLOGICAL CRITERIA FOR THE 
SELECTION OF REFRIGERANTS 


Present tendencies in manufacturing of the cooling devices 
and installations are closely connected with the use of new 
effective refrigerants. These refrigerants should also meet the 
requirements concerning atmosphere protection [4]. Ecological 
problems with the use of older refrigerants have led to the 
situation in which the chemical industry offers a large and 
still growing number of new substances, with a guaranty of 
their favourable thermal, flow and ecological characteristics. 
Unfortunately, we cannot be sure that no unexpected and 


unfavourable phenomena will be observed in the future, some 
time after these substances start being practically used [4]. In 
the Annex to the standard PN-EN 378 [3] the majority of the 
new refrigerants are taken into account, including criteria of 
the assessment of their influence on the environment. 

We all are aware that the Earth’s atmosphere should be 
protected. It was particularly evident for the signatories of the 
“Montreal Protocol” of 1987, as a consequence of the increasing 
area of holes in the ozone layer in previous years. Consequently, 
the signatories agreed on substantial reduction of the emission 
of CFCs until 2000. Poland signed and ratified the Montreal 
Protocol in 1990, and in 1996 the Act No. 176 was adopted 
according to which Republic of Poland signs the London and 
Copenhagen Annex to the Montreal Protocol. Initially, the 
negative effect of the halogen derivative compounds to the 
Earth’s atmosphere was described by two coefficients. The first 
is the Ozone Depletion Potential, ODP, which characterizes the 
effect of a substance on the intensity of ozone decomposition, 
while the second is the newly introduced GWP (Global 
Warming Potential) factor, calculated in the relation to CO, 
for which it is equal to 1 (Tab. 1). The carbon dioxide reveals 
a relatively low potential for creating the greenhouse effect, 
but its amounts emitted to the atmosphere are so huge that they 
contribute in about 50% to this effect. On the other hand, despite 
their small volumetric proportions in the atmosphere, vestigial 
gases, such as methane and halogen derivative compounds of 
CFC and HFC type, contribute in nearly 20% to the creation 
of this phenomenon [4]. Unfortunately for many Freon 12 and 
Freon 22 substitutes, which are expected to be safe for the 
ozone layer, reveal high GWP values, close to those revealed 
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by the withdrawn Freons. The GWP value is estimated for the 
time of 100 years, as the stability of the existence of the given 
substance in the atmosphere is to be taken into account. 

The sole information about GWP of the refrigerants turned 
out insufficient, as it did not take into account differences in the 
amount of driving energy consumed by the devices working 
with different substances. That is why the new factor includes 


the total level of CO, “emission” taking also into account 
refrigerant losses and the energy used for driving the device. 
The energy consumed for driving the cooling devices leads 
to the emission of CO,, the amount of which depends on the 
energy source. For instance, the combustion of different fuels 
generates different volumes of CO,. Also when the energy is 
generated without combustion, we will obtain different values 


Tab. 1. ODP and GWP factors describing ecological threat of selected refrigerants to the environment [3] 


ODP value 
(R11 = 1) 


GWP value (C0, = 1, 


Chemical formula 100-years period) 


Name of substance Refrigerant 


Chlorodifluoromethane 


Fluoromethane 


1,1-difluroethane 
1,1,1,2,3,3,3- 
Heptafluoropropane 


Octofluoropropane 


Octofluorocyclobutane 


2-metylopropane 


Propane 


Propylene 


Refrigerant 


44.096 
- 187.67 
-42.09 
492.08 


Molecular weight [g/mol] 


Freezing point [°C] 


Boiling point [°C] 
Density of liquid (25°C) [kg/m*] 


Density of vapour (25°C. 
1.013bar) [kg/m?] 


Critical temperature [°C] 


1.8314 


96.675 
4.2471 
218.5 


Critical pressure [MPa] 


Critical density [kg/m*] 


Heat of vaporisation (1.013bar) 
[kJ/kg] 


Thermal conductivity of liquid 
(25°C) [mW/mK] 


Thermal conductivity of vapour 
(1.013bar) [mW/mK] 


233.75 329.92 | 105.20 | 131.42 | 116.75 


83.479 97.975 | 44.753 | 59.402 | 65.136 


7.0475 9.4093 | 8.1634 | 10.415 | 10.212 


c 

1.1748 
Cy 

0.82591 


= € 
Cc. =1.8 1.135 


C.= 
1.16910 
0.80403 


Specific heat of liquid (25°C) 
[kJ/kgK] 


C= 
1.2568 
C= 


v= 
0.69086 


vo Cas 
neces 0.79750 


C,= (8) € 
0.76620 | 0.76641 | 0.74488 

Cy = Cy Cy 
0.64870 0.70035 | 0.69390 


166.92 | 239.46 | 362.88 


c C 
Specific heat of vapour 0.60629 0.99657 


(1.013bar) [kJ/kgK] Cy Cy 
0.49043 | 0.85355 


164.39 | 77.997 


Cy 
0.81640 
163.16 


Viscosity of liquid (25°C) [Pa s] 
Viscosity of vapour (1.01 3bar) 


[1Pa s] 


Cp/Cv (vapour) 
[25*C. 1.013bar] 


9.6967 | 7.0552 | 9.4093 | 9.7823 | 9.9046 | 10.230 


1.1847 | 1.2927 | 1.1549 | 1.0655 | 1.0150 | 1.0642 
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of the factor describing the volume of the produced carbon 
dioxide [kg CO,/1kWh]. Its assumed values are the following: 
burned coal z = 1.12; fuel oil z = 0.94; natural gas z = 0.57; 
wind, water and nuclear power plants z = 0.00. For Poland it 
is assumed that z ~ 0.7 — 0.8 [4]. 

Analysing different solutions and refrigerants used for 
instance in marine cooling installations we can assess the TEWI 
(Total Equivalent Warming Impact) factor for them and select 
those revealing its lowest values. 

Cooling installations and/or heat pumps contribute to the 
creation of the greenhouse effect, mostly in an indirect way via 
CO, emission connected with the production of the energy used 
by their driving systems. Since this energy is most frequently 
generated in Europe in electric power plants in which mineral 
fuel is burned, each kilowatt hour of the electric power is 
connected with the approximate emission of 0.4 + 1.2 kg of 
CO, to the atmosphere, depending on a country. 

In this way, the consumption of the driving energy by the 
cooling system contributes to the creation of the greenhouse 
effect during the entire time of system operation. Consequently, 
we should attempt not only to use substances revealing low 
GWP, which reduces part connected with the emission, but also 
tend to obtain the highest possible thermodynamic efficiency. 
An important feature of the present pro-ecological solutions is 
the use of highly effective compressors, fans, pumps, or other 
subsystems of the cooling systems or heat pumps. 

For simplicity purposes, it was assumed in the present 
analysis that the installations are identically filled with the 
medium and the level of leakage is in constant proportion to the 
installation filling. In practice, leakage flows of the refrigerant 
are different for different devices and substances, and the risk 
of their appearance is especially high in cooling installations 
which are spread over a large area. 

A high worldwide effort is made to reduce the emission of 
greenhouse gases, as a result of which some legal regulations 
were adopted. Since 2007, on the EU territory the Directive 
[6] on some fluorinated greenhouse gases is in force, which 
imposes severe limits to be complied with in the cooling and 
air conditioning engineering. (Directive no 842/2006 of May 
17, 2006, Polish text published in the EU Official Journal 
No. L161/ of June 14, 2006). 

The Total Equivalent Warming Impact (TEWI) factor is 
calculated from the formula [3]: 


TEWI=GWP-L-n+GWP-m-(1-)+n-E-z 


where: 
GWP - Global Warming Potential, in relation to CO, [-] 


L - amount of the substance emitted to the atmosphere 
[kg/year] 

n - devices operating time [years] 

m  - amount of the refrigerant in the installation [kg] 

f - dimensionless number assessing the recovery ratio [-] 

E - energy consumed by the device in one year (E = tp N), 
[kWh] 

Z - CO, emission divided by the scaling energy unit 
[kgCO,/kWh] 

N - heat pump driving power [kW] 

tp - annual time of operation [h]. 


Sample TEWI calculation for a heat pump 


For comparison purposes, the TEWI factor for the heat 
pump installation was calculated taking into account the 
following refrigerants: R22, R41, R152a, R218, R227ea, 
RC318, R600a, R290, and R1270. The operating parameters 
and technical data of the installation were the following: mass 


of refrigerant in the installation m = 5kg, annual emission of 
the substance L= 1kg, operating time of the installation n= 15 
years, annual time of operation tg = 6000 h, medium recovery 
ratio f = 0.5, heating power of the heat pump = 25kW, 


Q 
z = 0.94 kg CO9/kWh. The assumed identical filling of the 
installation and similar emission in case of each refrigerant 
are some simplification of the real situation. 

The calculations of a single-stage cycle with SLHE 
(suction-liquid line heat exchanger) were performed based on 
the assumptions collected in Tab. 3. The analysed refrigeration 
cycle is shown in Fig. 1. 


Fig. 1. Single-stage cycle with SLHE in the lgp-h diagram: 
i= evaporation temperature, tg- condensation temperature, 
Atzpg — superheating of vapour in the evaporator, At¿— extra cooling 
of liquid in extra cooler, t _- temperature of vapour at compressor 5 
~ suction side 


Tab. 3. Initial calculation data 


Data of calculation 


Evaporation temperature ty [°C] 


Condensation temperature tę [°C] 


Superheating of vapour in the 
evaporator Atzp [K] 


Resulting from 
SLHE efficiency 
0.80 


Sub-cooling of liquid in the SLHE 
Atg [K] 


Temperature of vapour at 
compressor suction side t, [°C] 


SLHE efficiency nwrp [-] 0.8 


20 


Isentropic efficiency of compression 


0.85 
Niz [-] 


Compressor’s mechanical efficiency 
Nm [-] 


0.9 


Tab. 4 collects the results of the TEWI calculations for the 
following refrigerants: R41, R152a, R218, R227ea, RC318, 
R600a, R290, R1270 and R22. 

Fig. 2 shows the relative TEWI factor calculated for the 
refrigerants: R41, R152a, R218, R-227ea, RC318, R600a, 
R290 and R1270, taking R22 (100%) as the reference. After 
assuming identical conditions for: the heating power, filling 
of the installation with the medium, the scale of leakage, and 
the operating time, the calculations have led to nearly 21% 
reduction of the TEWI factor for R600a compared to R22. The 
second favourable refrigerant was R152a, which revealed the 
TEWI factor lower by 19% than R22. Similar TEWI results 
were also obtained for the refrigerants R290 and R1270, for 
which TEWI decrease amounted to 18 and 17%, respectively. 
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Tab. 4. Results of heat pump cycle and TEWI calculations for refrigerants: R41, R152a, R218, R227ea, RC318, R600a, R290, R1270 and R22 


Refrigerant m [kg/s] 


qk [kJ/kg] 


R41* 


GWP*m*(1-f)| n*E*z 


1 455 479 414 481 112 


R152a 


2 100 365 774 368 224 


R218 


105 000 389 140 511 640 


R227ea 


43 500 365 434 416 184 


RC318** 


130 500 363 136 515 386 


R22 


25 500 377 465 407 215 


R600a 


45 : 356 777 356 829 


R290 


45 : 371 752 371 805 


R1270 


45 : 377 643 377 696 


* - not mentioned in the list in prPN-EN 378, ** - missing safety group in Table El in prPN-EN 378 -1 
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Fig. 2. Total heat pump TEWI factor for refrigerants: R41, R152a, R218, 
R-227ea, RC318, R600a, R290, R1270 and R22 


The last favourable result in this group was the refrigerant 
R227ea which revealed lower TEWI by 8%. The remaining 
refrigerants have higher TEWI values: R41 by 6%, R218 by 
13%, and RC318 by 14 % with respect to the referential TEWI 
factor of R22. 

Fig. 3 shows that the proportion of the energetic component 
in the TEWI factor amounts to 100% for R600a, R290 and 
R1270, 99.6% for R41, and 99.3 % for R-152a. In those 
cases the emission and recovery components do not exist or 
are negligibly small, below 1% of the total TEWI value. For 
the remaining refrigerants the emission component and the 
recovery component contribute to a much larger extent: 29.5 
% for RC318, 23,9 % for R218, and 12.2% for R227ea. 


COMBUSTIBILITY OF REFRIGERANTS 


It has been known [2] that the refrigerants R218 and 
R-227ea are not combustible under the atmospheric pressure 
and in temperatures not exceeding 80°C (except R152a). 
However, the tests [5] have also proved that even R22 
becomes combustible at the pressure exceeding 5,15 bar and in 
regular ambient temperature when it mixes with the air to the 
concentration equal to or higher than 65 per cent. Consequently, 
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Fig. 3. TEWI components for R41, R152a, R218, R-227ea, RC318, R600a, 
R290, R1270 and R22 


these substances should be used in the way protecting against 
the creation of air mixtures, during leakage examination, 
for instance. In general, we should not accept conditions 
for possible creation of such high-concentration mixtures at 
pressures exceeding the atmospheric pressure. 

The refrigerants R600a, R290 and R1270 belong to the 
group of media revealing increased combustibility (the lower 
combustibility limit is lower than 3.5% of the volumetric 
concentration in air). Moreover, R152a is a combustible 
medium, which was proved in relevant tests. At lower 
temperatures, higher pressure is required for the refrigerant to 
become combustible. That is why these refrigerants should not 
be mixed with air to examine the leakage flows. The problem 
of combustibility of the refrigerants and the resultant demand 
for relevant protection of the installation is an important, but 
fully realisable issue. It is also noteworthy that R717 (NH) is in 
common use in land cold stores, and even in marine refrigeration 
plants, despite their health threat and combustibility in certain 
mixtures with air. 

Below given is the classification of the refrigerants R41, 
R152a, R218, R-227ea, RC318, R600a, R290 and R1270 with 
the reference to the safety groups of these substances according 
to prPN-EN 378-1 [3] 


Refrigerant 


Tab. 5. Classification of selected refrigerants by safety group (Annex E to prPN-EN 378-1) 


Safety group 
(PN-EN-378) 


Fluid group 
according to 
PED 


Practical 
concentration 
limit! (kg/m3) 


ATEL/ODL2 
(kg/m3) 


Combustibility 
LFL3 


(kg/m3) 


Self-ignition 
temperature 
CC) 


1 - lower combustibility limit, 

2 -acute toxicity exposure limit ATEL or oxygen defect limit ODL, 
3 - lower combustibility limit , 
4 


- missing data (according to PN-EN 378, the standard does not refer to substances without attributed safety group) 


CONCLUSIONS 


The obtained results reveal that differences up to 35% can 
be obtained when calculating the TEWI factor for a heat pump 
with different media used by it. Therefore this factor can be 
considered an applicable and relatively simple criterion for 
assessing the performance of refrigerators or heat pumps. 
The refrigerants which occupy top ranks in this analysis are: 
R600a, R152a, R290, R1270 and R227ea; the results of their 
comparison are better than those obtained for R22. According 
to [1], R152a is a long-term substitute for R12 and R22; while 
R600a is a long-term substitute for R12, and R290 and R1270 
are long-term substitutes for R22 and R502. 

1. Unfortunately, despite their good properties for air 
conditioning, wider use of R600a, R152a, R290, and R1270 
is limited because of their high combustibility (R152a - 
group A2, remaining — group A3). These media have a very 
low TEWI factor, which suggests that they are likely to be 
used in individual design solutions. 

The problem of protection against the danger resulting from 

their high combustibility can be solved in many cases using 

relatively simple design means without excessive rise of costs. 

In this case controlling the filling of the installation is necessary 

(it should be done following the regulations of the standard 

PN-EN 378-1 Annex C, C3). When the device is situated in 

the open area, it does not create serious threat. R717 (NH,) is 

used in land cold stores and marine refrigeration plants. 

2. The use of R227ea is also a correct choice, as its TEWI 
factor is relatively low. However, high value of the 
GWP factor, equalling 2900, may result, in the future, in 
prohibition of the use of this medium [6]. 

3. Refrigerants R41, R218, RC318 are less favourable taking 
into account the TEWI values. They also have low value 
of the cycle performance, COP. Another disadvantage of 
R41 is very low critical temperature, equalling to 44.5°C. 
Moreover, there are no data on the safety group (the medium 


is not classified in the standard worked out in 2009!), which 
may suggest that no company submitted it for classification, 
or there are no relevant tests to allow this medium to be used. 
The data are also missing on the safety group for RC318. 
4. R218 and RC318 reveal high TEWI factors, they are also 
the substances with high GWP factor (GWP, = 7000, 


R218 


GWP.,.3;, = 8700), which practically eliminates them from 
wider use. 
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Problems of welding in shipbuilding - an 
analytic-numerical assessment of the thermal 
cycle in HAZ with three dimensional heat source 
models in agreement with modelling rules 


Part Il 
An analytical assessment of thermal cycle 
by used C-I-N and D-E heat sources models 


Eugeniusz Ranatowski, Prof 
University of Technology and Life Science, Bydgoszcz 


ABSTRACT 


This part is continuation of PART I. The basis of this analytic solution are the Fourier - 
Kirchhoff partial differential equation with appropriated boundary conditions. For a plate 
with optional thickness, the radiative heat transfer on both surfaces is taken into account. 
It is assumed that moving C-I-N or D-E heat sources during a very short period of time, 
generate an impulse of energy inducing an instantaneous thermal field in the plate area 
and the analytic solution is received by used Fourier transformation. These fields are 


t 


being continuously summed up to obtain resultant thermal field fr = J dT w) Finally, the temperature 
0 


fields generated by C-I-N and D-E heat sources in both stationary and moving co-ordinates systems are 
established. 


Keywords: welding; shipbuilding; welding in shipbuilding; thermal cycle; heat affected zone; heat source model 


Putting it into F-K equation the following form is 


LINEAR ANALYTICAL HEAT FLOW 


SOLUTION FOR A PLATE WITH received!: 
OPTIONAL THICKNESS AND RADIATIVE eT TT q -8(t) 
HEAT TRANSFER ON SURFACES BY C-I-N — t- t +t SS. 
HEAT SOURCE MODEL OX, OY, “OZ, A (3) 


Energy heat transport in HAZ is mainly progressed by 
thermal conduction and can be described by F-K partial 
differential equation. 


.[1—u(z—s)] exp k(@? + y2)-K, :2,]= ee 


Boundary conditions: 


Here is the way how the solution for temporary temperature =M=T = 
fields is received. The integral transformation method is being T(Xo Yo sZt =0)= Ty =0 (2) 
used. ôT 
HS power input in volume showed by eq. (12) — part I in —— = 0 when x, >%, x, > -00 (3b) 
impulse form is described as: 0 
e or 
"  m-[1-exp(K, -s)] z a) ay, =Q when y,>0, y, > -% (3c) 
-exp k(x’ + y*)—K, -z]-[l-u(z-s)] ƏT 
o FA when z,=0 (3d) 
Zo 
dy = ymax ` ÒC): O ar 
-exp k(x) +yo)—-K, -z,}[l-u(z-s)] ee when z,=g (3e) 
0 
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Graphic interpretation of the above conditions is shown 
in Fig. 1. 


Yo 
6 Xx 
= oT > 
A+ = Qt °T OT 
07 A = 
Zo dy, 


ôT 


0 A al T 
= L -— = -0 - 
Oy, 02, 


Fig. 1. Graphic interpretation of boundary conditions 


The x,, y, and z, variable transformations for equation 
(3) are executed with use the classic Fourier transformation. 
Transformation module for x, is described as”: 


K(, xX.) = exp: p: Xo) (4) 
The integration range is (— oo, + 00). The equation (3) after 
multiplying by (4) looks like: 


2 2 
“OX 20Yo 
ES pzz ip o T vax” O) 


soe (Oxy 
OZ. 2 ôx? A 


Pd, + 


(5) 


Ius) expEk()=K,: zoldx,= | = 


Transforming operation: 


dx, = 


+0 O°T 
Le 


sigo | . 


e dx) = 
Ox, 


g ipx, [t2 
—i- p: T-e? 0 
—0 


oT 


T. ed, | = — 


OX, 


i-p-Xo 


Xo =O 


—i-p-T-eP” 


Xo =- 


X¿=0 


+i-p-T-e Pr 


Xo == 


. | T-e'P dx, =0-0-0+04 


+(-p: fr . eP“ dx,) =-p° -T 


Transforming operations for: y,, Z, and t variables. 
These are obviously transformations for variables other 
than x, which simplifies the counting. 


Therefore for y, , z, variable we’ll obtain 


2 
Le JP dx, z 
0 
o? +00 : o? = (5b) 
OY abe dx) =2 GT 
0 0 
20°T ipx 
n 2 de , "dx, = 
OZ, 
pe (Se) 
reid. = - T 
OL; 


and for t variable: 
ie . ei Po dx, 
2 Ot 


Transforming operation for: 


2 O po i-pxXo _ ô= 
= Lts dx, == T (5a) 


m St) -[l-u(z-s)]- 
-exp-k(xp + yo)-K, +20] 
where: 
k-K,-Q 
x-[l—expCK, :s)] 


dmax = 


[cd 9) 


-exp k(x; + yo) -K, ze? dx, = 


= [E tmx 80O -[I-u(z=s)]- 


2: 
es - k(E+y9)-K, z 


(Se) 


and finally after first transformation the following is 
received: 


2 2 . 
a Tai OTs f7 a 80 
dys Zo aa (6) 
-[1-—u(z—s)]- exp- kE tyd- K; 2 l =T 
a ôt 


Similarly like for x, this classic Fourier transformation we 
will use for y, transformation. 
Transformation module for y, is described as: 


K(r, yo) =exp(i- p- yo) (7) 


The integration range is (— oo, + œ). The equation (6) after 
multiplying by (7) and transforming looks like: 


_p?- TH q? 2 Fe OT ans TSO), 


Oz, A-k (8) 
[1 —u(z—s)]) ex A oe _1ôT 
ak ak aA 
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Finally, transformation module for z, is described as: 
Q a 
K(r,z,) =cos(r-z,) + e -sin(t-Z,) 0) 
E 
Therefore: 


2,2, 2) mm e dumax o H(t) 
-° +q? +r a eee AA 
2 2 
A ~~ -K, 2} (10) 


On. A 
A r 2, dz, => 
-T 


and 


=p +a 44"): Taf a TN dymax T’ dmax TÒ), 


(11) 
. fexp- K,- Z) cosí- z,)+ 


0 


Let's find: 
s on 
[expe K,- Z) cos(r: z) + — 
0 


-sin(r-z,)dz, =D 
D =exp- K, s) 
—K, -cos(r-s)-A-r+r? 
(K +1 7)% 5 (12a) 
-A—Qyr:cos(r-s)-0a,K, -sin(r-s) K,-A+a, 
1 (K7 +17)-A 


.. T, are roots of : 


(12) 


-sin(r -Ss)- 


To B, r,, 


21.2 
At =O, 


cot(r - g) = —————_ 
A-r- (x +01,) 
So, the following is received: 


(12b) 


dmax * 7 Ó(t) 
-p +q +1?) T + da 2 OW) 
(p? +q* +r’): SK 


The solution of (13) is Green function corresponding to 
differential operator: 


dare) 


dt 
The Green's functions solution is as follows: 
= . . . 2 2 
T= _ a Q-K, -u(t) -exp| — P = A. A 
X-(1—exp(-K, -s) 4k 4k (14) 


-D-exp[-a-(p? q +1?)-t] 
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In order to obtain solution for T(x, y,z, t) - the reverse 
transformations must be provided. 

At first, reverse transformations with specified modules for 
“x,” and “y,” will be executed. 

"The a transformations modules are as follows: 


1 . 
R(p, Xo) oo AS) (15) 


1 
R(9,y0) og a (16) 


The integration range is (— oo, + œ) for both: “p” and “q” 
variables. 

In turn, the equation (14) after multiplying by (16) and ”p” 
retransformation looks like: 


Te AQ, uW D 
~ Vvn-A-[l-expCK, -s)] 


2 2 (17) 
q 2, 2 Xo 
exp| — a (q +r’) t- 
4-k a ) zis] 


Then, the equation (17) after multiplying by (16) and “q” 
retransformation looks like: 


T- Q-K,-k-u(t) l 
T:c,:(-exp(-K, s) (1+4-a-t:k) 


a (18) 
-D -exp a.t Et 
1+4-a-t-k 
where: 
C=C "Pp 


E] 


o variable is defined as: 


T= SER, Zo) 
i=l 
with z, retransformation module: 
RC, z) =B; C, 


and B, C, values are in agreement with (23a) and (23b). 
Finally, HS impulse temperature field is: 


E Q:K, -k-u(t) . 
mec, (1—exp(-K, -s):(1+4-a-t-k) 


op EL | 


1+4-a-t-k 
- y B,C/D; -exp(-a -r° - t) 
i=l 
Total temperature distribution from moving HS can be 
achieved by summing HS impulse results on its movement 
path: 


The last transformation for “z 


(19) 


(20) 


t 
T(t)= JaT(t’) 
Therefore, in stationary solidas system: 
Q:K, -k-u(t) | 
m-c, -[l—exp(-K, :s)-(1+4-a-(t-t')-k)] 


k. (x, -v-t'y +y le 
an|- 1+4-a-(t-t)-k 2 Be CDi 


(21) 


dT = 


(22) 


-exp(-a -r° 


(t=t') 


raf OK ku) 
o ac [lexp(-K, s) [+4-a- (t-t) k] 
k- (x, -v- Ey + yo] 
. A B,-C,-D.- 
an - 1+4-a-t-k 2 ' 
-expLa-1? -(t-t)] (23) 
where: 
B, = (23a) 
r 
2 
C 2i (23b) 


= x oye dor 2, . eh: 
OT ee ec ca UE 


K +H) Ay 
-A — Qy "1, -cosír, -s)- 0, K, sinf -s)]: + 039 
1 
P K,-:A+0, 
(SA 
T L, F,...r, are roots of : 
Na, 0, 


ctg(r, -g)= (23d) 
Ar (Qo +a) 
In moving co-ordinates system, x = x,-vt, Y = Yp Z= Z 
(Fig. 1, part I): 


0 


T=f'dt Q-K,-k-u(t) 
0 x-c,-[l—expCK,s)]-[1+4-a-(¢-t’)-k] 
AN 2 
exp EAN gy 
1+4-a-(t-t’)-k 
¿Y BCD -exp[-a: r°- (t-t')] 
i=] 
B, = (24a) 
where: ' 


C, D, r, - values are represented by same equations as for 
stationary system. 


The equations (23), (24) for assessment of the temperature 
fields in both stationary and moving coordinates systems have 
far more extended form than classical analytical solution of 
Rosenthal-Rykalin. Certainly it results from accepted heat 
source model C-I-N and characterising it parameters like: 
Q, k, K,, u(z-s). Furthermore, important elements of received 
solution are such parameters as: B, , C, , D, Their values also 
depend on such physical parameters as: A, a, a, and roots r, from 
equation (23d). However eq. (23), (24) have correct forms of 
solutions from mathematical and physical points of view but 
they have too compound mathematical form and they are hard 
for direct analytical account. Besides, analytical form of solution 
makes calculation impossible when taking into consideration 
non-linear form of thermal process under welding. So, above- 
mentioned affirmations make necessity of modification of 
analytical solution. This is possible through employment of 


hybrid analytical-numerical method. Present analytical form of 
solution defined by equations (23), (24) makes also impossible 
direct employment of numeric method for solution of these 
equations and assessment of the temperature fields. 


LINEAR ANALYTICAL HEAT FLOW 
SOLUTION FOR PLATE WITH OPTIONAL 
THICKNESS AND RADIATIVE HEAT 
TRANSFER ON SURFACES BY USED D-E 
HEAT SOURCE MODEL 


Similarly are estimated temperature fields for D-E heat 
source model but in a little complicated manner. At the 
beginning we must establish two partial differential heat flow 
equations for two different “quadrants” of another ellipsoids 
(Fig. 5, part I) at the same boundary conditions-equations 
Ga) + (3e): 

- x>0, (Fig. 5, part I) 


oT oT T Wl _ 1 oT 
La Soe (25a) 
Ax, ôy, an A a ot 
- x<0, (fig. 5, part I) 
oT oT T ¿dit Ler 
+ + (25b) 
Ox, ôy, an A a ot 


In stationary coordinate system the x,, y, and z, variable 
transformations for equations (25a), (25b) are executed in 
agreement with Fourier transformation as follows: 

a. for variable x, with transformation module K (p, x,)=e *° 
we will receive: 


- x,>0 
ala ol ə Ti of -Q-5(t) 
dy, az, tbc (26a) 
-3a (yac +z¿b") paí) 1 eT 
-exp — re p == — 
a, b*c 12 a ot 
- x,<0 
pT L "TA QW 
o êz, nabe (26b) 
g = 
46s -— 3a? a +Zob Zz pa; ,|_ 1 eT 
are 2 a Ot 


b. for variable y, with transformation module K (p, x,)= erro 
we will receive: 


- for x, >0 
= 22m 
TT oh 2V3f,-Q-8(t). 
OZ Jnc 
] — (27a) 
O Laal OT 
> Eas q) == 
axl Ap ) a at 
- for x, <0 
= 27 
pT- T21, 24/3£, -Q-5(t) 
OZ; mic 
— (27b) 
2 
zZ 1l sad Lor 
a q’b 
= ay) an|- at 
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c. for variable z, with transformation module 
K(r,z,) = cosrz, + yr MZ we will receive: 
r 

- forx,>0 

= 231, -Q-d(t). 
2 2 
TEP -g -r)+—__>— 
T Ac 


a- e a + wr) (28a) 
g 6 = 
1 oT 
cos(rz, ) + — sin(rz e © dz = ——— 
eose) Ein, aa 
- for x, <0 
Tep -g -ry 2B 80, 
Jarc 
1 59 34.2 
-exp| -— (p’a, +q"b")|- 
1 e +4 ) 
g ¿En = 
Jesz) t sinezo e “dz, = oe (28b) 
0 
During the Fourier transformation for the ,,z,” variable it 


appears in eq. (28a) and (28b) the integral: 
g 2 
Oy. 3z 
| cos(rz,)+— sin(rz,)|-exp) -=> ldz, (29) 
: Mr c 
The above integral can’t be solved immediately but an 
algorithm can be used to obtain a satisfactory approach. In order to 
obtain the solution, the function exp 3z¿/c?) may be written!: 


In particular, nlast may be small when calculating for thin 
plates with quite large z-semi axis, see Fig. 2. 
1 


0.5 


z (cm) 


Fig. 2. Graphs representing convergence of exp(-3z’/c?) and its approx 
function (30) with various nlast and c values: 1. exp(-3z’/c?), c = 3; 
2. approx(z, c, nlast), c = 3, nlast = 6; 3. approx(z, c, nlast), 

c = 3, nlast = 2; 4. approx(z, c, nlast), c = 7; 5. approx(z, c, nlast), 

c = 3, nlast = 6; 6. approx(z, c, nlast), c = 7, nlast = 


al 
372 e ) In general, a source of given c and plates with large 
exp =s Lsj NE = (30) thickness, require a high nlast parameter which makes the 
e? = | whole calculation longer. For example, analysing a plate 
[o with thickness g = 5 cm, penetrated by a “double ellipsoid 
n=1 configuration of source” with c = 2 cm, the approx function 
= approx(z, c, nlast) must be used with nlast = 13 (or larger) and is equal to: 
9 27 81 
approx = 1-—-z* + -z* -——~-2°+ —_-z* ZO 
(2:c’) se") (8-c”) (40-c”) 
81 243 729 243 
A Oe ee te a a 6D 
(80-c?) (560-c!*) (4480-09) (4480-c'8) 
729 2187 2187 6561 
1 ae = bo a ae ren s rT aaa 
(44800 -c”) (492800 -c”) (1971200-c”) (25625600-¢7°) 
So, a Fourier transformation for the z variable can be easily found: 
g 
Í cost zy 2 T 2 sin(r -z,) |: aprox(z, c, nlast)dz (32) 
0 


The integral equation (32) is computed with the use of computer symbolic calculation. The result is usually very long. This 


expression is signed as E! to perform the rest operations. 


Through execution the reverse transformation for specified modules for x, y, Zz, in the same way as previously for C-I-N 
H-S model the temporary temperature field generated by a pulsed “double ellipsoid configuration of source” is finally achieved 


as follows: 
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2 2 
qt 3/3 -exp|— Xo Jo + 


To AA 
n-Vn co, f2-a-t+a2)(2at +b?) fort +a} doit +ib; 
a 


2 2 2 
4 233 e M — +— . (33) 
rn e /12-0.-t+a?)(2at+b?) dot + za: dot+—b; 
OL 


J3 B;C;E; -expC a- t “t) 
i=l 


Using an additivity method, we may achieve the summary temperature field generated by moving heat source as follows: 
- stationary co-ordinates system: 


-£.-3-/3dt 
A =| dis 


0 molle, JOa (1) +a; 2a 0)+b 7D 


-exp|— A no o e Bd 
dat) 4-00) ob} nl, (M12-a-(t-t)+a/]-[12-a-(t-1)+b?) 


2 2 59) 
-exp| — EEE vty —— 7 -© B,C,- E, expa- -(t-t)] 
4-a-(t-t)+-al 4-a-(t-t’)+—b? i=] 
3 3 | (34) 


- moving co-ordinates system: x = X= VEY =Y 7 2 = 24 


TGx,y,2)=| qu£,-3-/3dt 


0 dro a Era [l2-a-G—t)+b2) 
a 


k+v (t-t y gidde 
wT lao sie ue a A l 
4-2 (1) a ta Gth mn Le -a{Ql2-a-(t—t)+ a2} 2-a- t-th+bD 
a 
ne 2 00 
-exp| — E —__y Y B,-C,:E,:expba-1(t-t)] 
4.a-(t-t)+-al 4-a-(t-t')+—b? al 
3 3 | (35) 
where: g 
E, af cost *Z)+ = «sin (1, * ») -approx(z, c, nlast)dz (36) 
-T 
0 
nlast (2) 
approx(z,c,nlast) =1+ > at d 67) 
n=1 


n 
i=l 
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Equations (34), (35) have a similar features as previously 
determined equations (23), (24) for assessment of the 
temperature fields in both stationary and moving coordinates 
systems. Both solutions have been given for three dimensional 
C-I-N and D-E heat source models but with different thermal 
characteristic. Especially it is obliged to differentiate low or 
high concentration of energy of heat sources. 

Remaining comments concerning manner of solutions of 
equations (34), (35) and their results are the same as presented 
previously for equation (23), (24). The employment of method 
of hybrid account exists also in this case. That is the main 
purpose of the III part of this work. For further calculations it 
is chosen the analytic-numerical method. 


CONCLUSIONS 


In this paper the new temperature evaluation during welding 
of a plate with optional thickness is presented. Furthermore, the 
radiative heat transfer on both surfaces is taken into account. 
In these calculations we used the following heat source 
models: cylindrical-involution-normal (C-I-N) and double 
ellipsoidal configuration of source (D-E). By using the Fourier 
transformation method, the temperature fields generated by 
above heat sources in both stationary and moving co-ordinates 
system are presented. The algebraic forms of analytic solutions 
presented by equations requires discretising in order to make 
possible numerical calculation. 


NOMENCLATURE 
T. T, T - the temperature transformation 
KG, x,), KC, y,), KC, z,) - transformations modules for x, , y, » Zo 
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R(t, x,), RG, yo), R(T, Z) - reverse transformations modules 
e, - volumetric heat, [J K“ cm] 
c - specific heat, [J K'kg'] 
P è 
p - mass density, [kg cm”] 


m 
[In - performs iterated multiplication of n over i= 1, .., m-1, [m] 
i=l 


B,, C, - values are in agreement with (23a) and (23b) 

nlast - natural positive number as large as is necessary to achieve 
the required approximation 

HAZ - Heat Affected Zone 

C-I-N- Cylindrical-Involution-Normal heat source model 

D-E - Double-Ellipsoidal heat source model 

HS - Heat Source. 
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ABSTRACT 


The purpose of this study is to estimate the wave height at the front face of breakwater (Refracted breakwater 
and Straight breakwater), when dredging like the submarine pit is performed in the distant offshore from outer 
breakwater. The wave field of the problem is considered to be two dimensional planes and the configuration 
of the pit region is designated by a single horizontal long-rectangular system. The numerical approach uses 
the Green function based on the boundary integral approach. The results of the present numerical works are 
illustrated by applying the normal and inclined incidence. It is shown that in the case of normal incidence, 
the ratio of wave height reduction at the front face of both types of breakwaters is approximately more than 
20% due to the effect of the submarine pit on the sea bed. Furthermore, regardless of the type of breakwater 
and the difference in incident wave angles, the ratio of wave height was shown to be reduced. 


Keywords: Submarine pit; Green function; boundary integral approach; breakwater 


INTRODUCTION 


Study of wave diffraction by submarine pits has important 
application to coastal engineering. The purpose of this study 
is to investigate the reduction effect of wave height by 
submarine pit using depth-discontinuity. The early studies on 
the interaction of incident waves with a dredged region has been 
investigated by Newman (1965), Hilay (1969), Lee and Ayer 
(1981), Miles (1982), Kirby and Dalrymple (1983), Ting and 
Raichlen (1986), and Kirby et al. (1987). In all of these studies, 
the pit region was assumed to be infinitely long and the problem 
was restricted to the horizontal and vertical coordinates. 

Later, Williams (1990) and McDougal et al. (1996) 
presented the numerical model for involving two horizontal 
planes in the case of single or multiple pits. The 3D (three 
dimensional) model, which is based on the Airy wave theory, 
was then investigated by Williams and Vazquez (1991) to 
analyze wave diffraction through a single pit. Recently, 
modeling diffraction of random waves by submarine pits has 
been investigated by Kim (2007). However, the focus in all of 
these studies were restricted to interaction of only one or two 
aspects, which occurred over pit (or trench) with discontinuity 
water depth. 

The present study is about the interaction of wave, which 
is propagated over pit with discontinuity water depth. This 
study investigates the affect of this wave on the font face of 
breakwater systems (refracted and straight breakwater) when 
a pit is dredged. The wave interaction is connected to three 


boundary problems: the interaction of discontinuity water 
depth boundary and pit boundary, discontinuity water depth 
boundary within pit and breakwater boundary, and pit boundary 
and breakwater boundary. The problem is considered in a two 
dimensional plane, and the configuration of the dredging region 
on sea bed is a single-long rectangular type. The numerical 
simulation is performed by using the solution of the boundary 
integral equation based on the Green function. 

In order to verify the present numerical model, comparisons 
made with the results of absolute and approximation solution 
presented by Koji and Mutsuo (1976) for regular wave 
diffraction at the front face of refracted breakwater without 
pit. Upon comparison, the present numerical simulation and 
the research by Koji and Mutsuo (1976) fell into relative 
agreement. 

Decreasing effects of diffracted wave field in vicinity of 
breakwater due to a pit with discontinuity water depth were 
shown. The present study can provide information to design 
the dredge line of the outer breakwater that is distant from this 
dredge line. This study can also be effectively utilized for wave 
interaction by pit and be applicable to coastal engineering. 


THEORETICAL DEVELOPMENT 


The geometry of the problem is presented in Fig. 1. The 
fluid domain can be separated into three regions: the interior 
pit region with uniform depth d, the refracted breakwater, and 
the exterior fluid region with uniform depth h. The boundary 
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Fig. 1. Definition sketch for fluid domain and boundaries 


regions express S, (pit boundary) and S, (breakwater region), 
respectively. In Fig. 1 the Cartesian coordinates are designated 
as having the origin at a corner of the breakwater. The x- and 
y- are directed in the horizontal plane and the z-axis is directed 
vertically upward from at the equilibrium water surface. 

The fluid region is taken to be inviscid and incompressible, 
and the flow irrotational, then it will be assumed that the fluid 
motion may be described in terms of a velocity potential 
p(x, y,z,t) = Re[D;(x, y, z)e "| for j = 1,2. This potential 
must satisfy Laplace’s equation, V'®;=0. It is subject to the 
usual boundary conditions on the free surface and seabed, and the 
solution of velocity potential can be given following forms: 


® (x, y, z) = Y cosh[k (z + d)] onS, (1) 
D,(x, y, z) = Y cosh[k (z + h)] onS, (2) 


where, d is the water depth within the interior pit region, and h 
is the water depth of the exterior fluid region. The subscript 1 
represents the pit region; subscript 2 represents the exterior fluid 
region, and the wave numbers k((j = 1, 2) are defined by: 


œ? = gk,tanh(k,d) on s, (3) 
œ? = gk,tanh(k,h) on S, (4) 
where: 
o - angular frequency 


g - gravitational acceleration. 


The governing equations in each of the fluid regions should 
be satisfied by the Helmholtz equation. 


1 
+ (P)+ ae E 


1 
P (P)-29,(P)-— S 
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+kiW=0:j=1,2 (5) 


Continuity of mass flux and pressure across the fluid 
interface between the interior pit region (S,) and the exterior 
region (S,) requires the following conditions to be satisfied: 


ov, = h o on Ss, (6) 
on on 
Y=, ons (7) 


Finally, the scattered component of the fluid potential in 
the S, is subject to a radiation or far-field boundary condition 
at large radial distances r, which may be written as: 


lim Vi Ei, Jo, -D,)=0:i=1,2 (8) 
to r 3 
where O, is the incident potential function and is given by: 
igH, cosh[k,(z+h 
20 cosh(k,h) 


Green’s function G(P, Q), j = 1, 2 may now be expressed 
for each region as: 


G;(P,Q)= FHP (KR) 5-1 (10) 


where, H P is the Hankel function of the first kind of order 
zero, P and Q have the coordinate of (x’, y”) and(x, y) on the 


boundary(S,, S,), respectively, and R=/(x9 + (y'-y)’. 
Applying Green’s second identity to Y, and extending to 
full regions (S, + S, + S,), boundary integral equations can be 


yielded as Eq. (11) and Eq. (12). 


ôG; l 
0 eo-a eo Z) ls 0 = 152. Un) 


o eo-ae Fe fs, =0;j=1,2 (12) 


The scattered wave has no effect on the imaginary boundary line (S,), and applying boundary condition of (6) and (7) proposed 
by Williams (1990) and McDougal et al. (1996) to Eq. (11) and (12), and Eq. (11) and (12) become: 


5 (P) += al o Te, 0-36, (P,Q): (Py as, + 


1 


¥o(P)-20,(P)-— | va 


+= ~ o 


The free surface elevation in each region can be obtained 
using Eq. (15): 


i ;j=1,2 (15) 


in which n(x, y, t) = Re {A, (x, yje“, where A is the spatial 
component of surface elevation in fluid region. Finally, the 
diffraction coefficient of regular wave, K, is defined as: 


K¿= (16) 
NUMERICAL RESULTS AND ANALYSIS 


Validation of the numerical model 


A computer program has been developed to implement 
the above theory for regular wave diffraction at the vicinity 
of breakwater by pit on the sea bed. To verify the present 
numerical results of the wave height at the front face of refracted 
breakwater, these were compared with the absolute and 
approximate solutions presented by Koji and Mutsuo (1976). 
The conditions of calculation and configuration of refracted 
breakwater in regular waves without pit by Koji and Mutsuo 
(1976) are defined as incident wave angle 0 = 45°, angle of 
refracted breakwater B = 135°, reflection coefficient K, = 0.8, 
and length of breakwater = 2L (L,+L, = 2L), respectively. 

Fig. 2 shows a comparison of the results at the front face of 
refracted breakwater for the diffraction coefficient of regular 
waves (CD line) obtained by the present numerical model 
and those of Koji and Mutsuo (1976). In a comparison from 
fig. 2, it is noticed that reasonable agreement is shown in the 
three cases for the refracted breakwater, but the diffraction 
coefficient at the front face of refracted breakwater by the 
present numerical model is a little bit higher than the results 
of Koji and Mutsuo (1976). They had performed their studies 
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Fig. 2. Comparison of wave diffractions of present study with those of 
approximate and absolute solutions at the front face of refracted breakwater 
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at the infinite refracted breakwater but the present numerical 
model has been investigated at the finite refracted breakwater, 
so in the case of finite breakwater, the wave interface has been 
occurred by phase difference and interaction between incident 
waves and reflect waves. 


Numerical examples 


Numerical examples are presented to investigate the influences 
by rectangular submarine pit at the front face of refracted and 
straight breakwaters with different incident wave angle. 

In this study, the condition of calculations are water depth 
of the vicinity of breakwater h = 7 m, water depth within 
submarine pit region h = 14 m, the refracted angle within the 
refracted breakwater B = 160°, and incident wave angle 0 = 90° 
and 0 = 135° for the regular wave, respectively. Figs 3 and 4 
present the results for the diffraction coefficient obtained by 
the present numerical model with and without submarine pit 
at the front face (CD line) of refracted and straight breakwater, 
respectively. 

In the case of normal incident wave angle 0 = 90°, wave 
height reduction of refracted breakwater due to submarine pit at 
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Fig. 3. Wave height distributions at the front face of breakwater systems 
with or without pit (incident wave angle 0 = 90°) 
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the offshore sea bed as shown in fig. 3(a) can be observed at the 
vicinity of corner side. The wave height reduction at the vicinity 
of the corner side is shown to be the largest. This means that 
the wave energy is concentrated on the corner of the refracted 
breakwater, and by using the pit, the wave energy is subdued 
thus reducing damage on the breakwater. Fig. 3(b) shows that 
the straight breakwater is reduced in the entire region. This 
reduction is 26.5% less than the numerical results for the case 
without submarine pit for the refracted breakwater, and 41.7% 
for straight breakwater. 

As for the case of inclined wave angle 0 = 135°, wave 
height reduction due to pit can be estimated 12.6% for refracted 
breakwater and 5.0% for straight breakwater. As shown in figs 
4(a) and (b), reduced wave height is expressed at the front 
face of breakwater except at the left side of each breakwater. 
Table 1 presents the effect of added pit on the reduction in 
wave height line CD. 

From these analyses, we find out that the overall reduction 
of wave height may be observed at the vicinity of breakwater 
with an appropriately set submarine pit. 
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Fig. 4. Wave height distributions at the front face of breakwater systems 
with or without pit (incident wave angle O = 135°) 


Tab. 1. The ratio of reduction in wave height 


Type ot breakwater Refracted Straight 


3 breakwater | breakwater 
Incident wave angle 


Fig. 5 presents the contours of the wave height ratios near 
the breakwaters due to regular diffraction wave with two case 
of breakwater. We can derive from the results that a dredging 
pit may provide an excellent means of protection from 
a wave attack. It can be seen that wave energy is weakened 
by discontinuity water depth (dredging). The results from 
numerical simulation for wave diffraction indicate that the 
present model is appreciable to estimate the wave field at the 
vicinity of breakwater when a pit is dredged on the offshore 
sea bed. 
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Fig. 5. Diffraction contour plots at the vicinity of breakwaters with or 
without pit 


CONCLUSIONS 


The purpose of this study was to estimate the decreasing 
effects of diffracted wave fields around breakwater when 
the dredging work (pit) is performed at the offshore sea bed. 
This study is about the composite interaction for the three 
problems: dredging boundaries, depth discontinuity of pit, 
and breakwater boundaries. The boundary conditions are 
established for those problems and applied to the boundary 
integral equation. 

Results for the incident wave conditions have been presented 
to illustrate the wave height distribution of the wave field 
near breakwater by the influence of the pit. The results 
of present numerical model have been compared with 
those from published calculations and present numerical 
simulations that show relative agreement with published 
data for diffraction. 

Through the present numerical simulations the wave height 
reduction of refracted and straight breakwater with different 
incident wave angle due to submarine pit at the offshore 
sea bed can be observed. 

It is noticed that the results from the present numerical model 
accurately provides the diffracted wave height of vicinity of 
breakwater systems, when the pit is dredged on the offshore 
sea bed, and so may be extended to apply with confidence 
in breakwater planning and design applications. 
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ABSTRACT 


This paper presents a Hybrid APF/PFC/PPF Circuit (active power filtering and power factor correction 

circuit aided by shunt passive filter) for harmonic suppression and power factor correction in a marine 

electric network. By employing the proposed hybrid circuit, marine electric power network voltage and 

current can be maintained as being sinusoidal, and the power factor is close to 1. The effectiveness of 

the proposed method that is applied to a marine electric network is demonstrated through a simulation 
experiment. The results show that the expected performances are achieved. 


Keywords: Electrical engineering, marine vehicle power systems, power system harmonics, 
power system measurements, quality control 


INTRODUCTION 


Recently, application of power electronic devices to marine 
systems has been wider and wider, which can be seen from 
the fact that the number of power electronic devices equipped 
in ships is increasing very quickly and the capacity of these 
devices in ships is becoming larger and larger. As a result, in 
a ship, electrical power energy consumed by power electronic 
devices occupies a continuously increasing proportion of the 
gross electrical energy produced by the ship. But harmonics 
produced by power electronic devices such as frequency 
converters or soft starters cause distortion of the current 
waveform and voltage waveform of the electric network, as 
well as having a remarkable influence on safety and capacity 
in the choice of electric network, which makes harmonic 
pollution in a marine electric network more and more serious 
and difficult to estimate [11, 16]. Furthermore, impedances of 
generators designed for ships are relatively higher than those 
of overland power sources, which worsens current distortion 
in marine electric networks. Another problem of marine 
electrical power quality is a low power factor. The selected 
capacity of ship equipment is usually very conservative and 
its power factor during routine operation is correspondingly 
lower since a marine electric network is a kind of low capacity 
isolated network when compared with an electric network that 
operates on shore. But the reliability of ship operation plays 
a vital role in safety of navigation, which leads to a substantial 
safety consideration in the design ofa marine electric network. 
This is evident in the rules of marine electric network power 
quality issued by the majority of shipping classifications in 
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the world, which are not less strict, or are even stricter, than 
the rules of public electric network power quality issued in 
international standards [4, 13, 14, 17]. For example, the limit 
of stable voltage adjusting rate is between +2.5% and -2.5% 
as ruled by the China Classification Society, and the limit of 
stable voltage fluctuation is 3.3% in IEC EMC 61000-3-11 
[4, 5]. As another example, the limit of total voltage harmonic 
distortion (THD, ) is 5% as shown in the national standards of 
low voltage (0.38kV) public electric network of China, while 
comparatively, the mentioned parameter ruled by DNV (Det 
Norske Veritas — Norwegian classification society) and ABS 
(the American Bureau of Shipping) is 5% too. But the fact 
is, it is more difficult to obtain an improvement of the power 
quality of a marine electric network than that of an overland 
electric network due to the characteristics of marine electric 
networks [10]. 

Hitherto, the majority of achievements in the field of product 
design for power quality improvement have been restricted 
to overland electric networks. There is still a lack of reports 
showing significant progress of power quality improvement 
technology applied to ship electric networks. Nevertheless, 
even considering the practical products which have been 
applied to overland electric networks, their applied effect has 
not been satisfactory due to the limitation of control effect and 
the accuracy of measurement of the electric network in the 
aspects of active power, reactive power, power factor and high 
order harmonics. That is, it is reasonable to design a new circuit 
for a marine electric network to depress harmonics and correct 
the power factor under the consideration of characteristics of 
a marine electric network. 


HYBRID APF/PFC/PPF CIRCUIT 


Definitions of harmonic distortion 
and power factor 


Normally, distorted power, based on Budeanu’s concept, 
is defined as [18]: 


Dis ¡$ -P?-0* (1) 

The most popular parameters describing the above- 

mentioned distortions are the total harmonic distortion 

coefficients THD, and THD,, voltage and current respectively, 
as shown in formulas (2) and (3). 


THD, == — (2) 


THD, = A (3) 


where V, is root mean square value of the fundamental 
harmonic, n is order of harmonics, and V „is harmonic voltage 
when the order is n. Normally, THD, and THD, are obtained 
when n a is 40 or 50, where n „is the highest harmonic order 
in the process of defining THD, or THD.. 

The definition of power factor is shown in formula (4), where 
total power factor 2 is the arithmetic product of displacement 
power factor à and the distortion power factor à. Formulas 
(5) and (6) show the definitions of A, and À respectively. The 
total power factor à can be described as formula (7), where 
S, is the apparent power under fundamental frequency, P is 
active power, Q is reactive power, D is distorted power and 
S is apparent power. 


A=2, Àp (4) 
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The design of harmonic suppression and power factor 
correction solutions for ship use aims to minimize the quantities 
of Q and D in equation (7) to improve the total power factor A 
in the marine electric network environment [9]. 


Design of Hybrid APF/PFC/PPF Circuit 


In this paper, the designed Hybrid APF/PFC/PPF Circuit is 
a combination of SAPF (series active power filter) and active 
PFC (power factor correction circuit) aided by a shunt passive 
filter, which is shown in Fig. 1. 

Application of this design takes into consideration that 
a marine electrical power system is very different from 
a land power system, because of its character, i.e. an isolated 
power system, where the capacity of a heavy load in a ship is 
comparable to that ofa ship electrical power station. Moreover, 
the impedance of generators designed for ships is higher than 
that of a land power source, which worsens voltage deviations 
in marine electric networks. Also a heavy load gives an electric 
network a great impact and causes a considerable voltage 
drop and frequency oscillations of the electric network at the 
moment it is switched on. At the same time, power electronic 
devices generally are used in ship systems and cause distortion 
of current, and in consequence, voltage waveforms. Taking into 
account the complexity of the considered matter, it means that 
the load variation range and variation frequentness, and many 
kinds of disturbances covering the current as well as the voltage 
character (including their harmonics) must be considered. At 
the moment to mitigate harmonics, mainly passive power filters 
(PPF), and rarely parallel active power filters (APF) in marine 
technology are used. To improve the dynamic properties of 
a marine electrical power station and to reduce the power of 
filtering circuits, a hybrid filter concept for marine electric 
networks was investigated. 

Anew type of hybrid solution for a marine electric network 
was considered as a study and research for the general case 
of a marine electric network, including voltage sensitive (in 
the meaning of voltage quality) loads as well as voltage non- 
sensitive loads. Both of them may be easily appointed in the 
marine networks, e.g. all computer-aided control, alarm and 
monitoring systems are strongly sensitive voltage loads, while 
in contrast the heaters are voltage non-sensitive loads. On the 
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Fig. 1. Configuration of Hybrid APF/PFC/PPF Circuit 
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other hand, power converters generate mainly current, but also 
voltage harmonics, so they generally require a compensation 
voltage and current harmonics. 

In Fig. 1, the SAPF and shunt PF constitute a combined 
system of series active power filter and shunt passive filter. 
Based on the following two reasons, the combined series active 
power filter is designed to be included in the Hybrid APF/PFC/ 
PPF Circuit. Firstly, the capacity of a marine power station, as 
well as the capacity of a single generator contained in a power 
station, becomes higher and higher, but due to the characteristics 
of a marine generator, source voltage harmonics are rather 
considerable when compared with those that can be found in an 
overland network. The SAPF is first of all a controlled voltage 
source, and it is connected in series between an electrical power 
supplier and a linear load. The authors have distinguished two 
kinds of load: one is the voltage sensitive load (VSL) and the 
second is the voltage non-sensitive load (VNSL). The SAPF 
can compensate voltage harmonics and improve quality of 
voltage supplied to the voltage sensitive load (VSL). That is, 
the SAPF is suitable for harmonic suppression in distorted 
supply voltages [3]. Secondary, large rated cycloconverters 
and thyristor converters are more and more frequently applied 
to marine systems whose capacity occupies an incremental 
percentage of the total capacity of electrical power. When the 
combined system of shunt passive and series active power 
filter is applied, the shunt PF acts as the main bearer to fulfil 
the task of current filtering, while the series APF does not act 
rather as a harmonic compensator but mainly as a harmonic 
isolator. It decreases source harmonics flowing into the shunt 
PF and causing resonance. Also the required rating of the 
series active filter will be reduced to less than 1%, as compared 
with the rating of a 3-phase 6-pulse thyristor converter or 
cycloconverter, that is, the required rating of series APF is much 
smaller than that of a conventional shunt APF [12]. Thereby, 
application of a combined system of shunt passive and series 
active power filter results in better filtering characteristics and 
lower initial and running costs. In the SAPF shown in Fig. 1, 
L, and C, constitute a set ofa small capacity filter which is used 
to suppress high frequency pulsating voltage and current caused 
by the high frequency switching of the converter. 

Acombined system of shunt passive and series active power 
filter is not sufficient for marine electric network application. 
There are two reasons to develop this filter structure. As 
mentioned in the first paragraph of this work, the problem 
of a low power factor may cause a serious problem in ship 
electric networks and even cause breakdown of an electrical 
power station during the starting process of bulk load. SAPF 
is able to improve the operational environment mainly for 
loads generating voltage harmonics. In contrast, a shunt APF 
is very effective in depression network current harmonics, 
usually generated by a non-linear load. Furthermore, the effect 
of network current harmonic mitigation is not so effective due 
to the inherent deficiencies of shunt PF, which, undoubtedly, 
contributes to suppressing the majority of current harmonics. 
In the authors’ opinion a solution is to include an active PFC 
(power factor correction) circuit in the Hybrid APF/PFC/ 
PPF Circuit. The active PFC circuit is expected to solve the 
problem of the low power factor in the electric network, and 
simultaneously acts as a current harmonic filter to suppress 
network current harmonics additionally. Taking that into 
consideration, the active PFC circuit is located between the 
SAPF and the shunt PF. 

While both the combined system of shunt passive, the series 
active power filter and the active PFC circuit are concerned, 
each of these works as a supplementary means of the other one, 
and they work together to suppress harmonics and compensate 
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reactive power at the same time more effectively, assuming 
that superiorities of different types of compensation devices 
are mutually complementary. It is also anticipated that, for 
a marine electric network — an independent small capacity 
electric network in a complex environment and requiring 
a quick response to harmonics suppression and reactive power 
compensation — the proposed solution will meet the demand 
of a marine electric network for high-speed and accurate 
measurement and management of harmonics and reactive 
power. Moreover, taking into consideration the possibility of 
harmonics suppression and reactive power compensation, the 
suggested solution will lead to the achievement of significant 
energy saving. 


Control algorithm of Hybrid APF/PFC Circuit 


The key technology to achieve a rapid and accurate 
effect of electric power quality compensation in ship electric 
networks, considering the application of the Hybrid APF/PFC 
Circuit, is the high operational performance of the SAPF and 
the active PFC circuit in the Hybrid APF/PFC Circuit due to 
the complicated operational environment of ship loads and 
high frequency of a wide range of variation of load operating 
mode. 

For an active power filtering or an active power factor 
correction circuit, the detection circuit of harmonics and reactive 
current is an important part that determines its compensation 
characteristics. Generally speaking, in practical application, 
there are two major types of harmonics and reactive current 
detection methods. One is based on series expansion of 
Fourier Transform [2], while the other is based on the theory 
of Instantaneous Reactive Power of a 3-phase circuit, also 
called the p-q method [1]. The former method is able to 
detect a detailed values of harmonics series or the harmonics 
of a specified frequency band, which qualifies it for fault 
diagnosis and protection of the electric network, but its real time 
operational property of current detection is not satisfactory. The 
latter method has quite good properties in real time detection 
but is not able to work effectively when the source voltage 
is distorted [19]. Nevertheless, the theory of Instantaneous 
Reactive Power of a 3-phase circuit has been the most used 
until now in the control of active power filters [15]. Of course, 
there are quite a number of other types of harmonics detection 
methods, e.g., the adaptive harmonic current detection method 
[8] and the predictive current detection method [6]. The fact is 
that all these proposed methods are still in the stage of research 
and development and lack practicability. A current detection 
method based on an improvement of the Instantaneous Reactive 
Power theory of a 3-phase circuit, which is called the i-i 
method, is being developed to effectively detect load current and 
divide the load current into the fundamental active component, 
fundamental reactive component and harmonic component with 
accuracy and real-time property even when the source voltage 
is distorted or 3-phase unbalanced [19]. Its working principle 
is expressed concisely in this paper by Fig. 2. 

In Fig. 2, voltage vector V, stands for phase A voltage of 
power source, with its instantaneous value v = V „sinot, where, 
V „ is amplitude of phase voltage. It provides a PLL (phase- 
locked loop) circuit with the phase information of sinot and 
cosot, with appropriately defined signs. In this detection and 
control system, three types of coordinates are involved, among 
which the “3-phase coordinate” and “af coordinate” are static 
coordinates and the “pq coordinate” is a rotary coordinate with 
its angular velocity œ — the angular velocity of source voltage. 
The quantities i,, i,, 1, are the detected 3-phase line currents of 
the electrical power network respectively. 


Fig. 2. Fundamental principle of the compensation current determined by i-i, method, 
a) detection of harmonic current, b) detection of the sum of harmonic current and reactive current 


In Fig. 2a, i,,, 1,,, 1,, are the 3-phase harmonic currents 
respectively and 1,,,1,,, 1 pare the 3-phase fundamental currents 
respectively. In Fig. 2b, Ly e. i „are the 3-phase fundamental 
active currents respectively, an les La, e are the 3-phase 
fundamental reactive currents respectively. 


The relationship of the variations in Fig. 2 is: 
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where C py is the inverse matrix of Cp and Cr Cop: Oz 15 the 
transposed matrix of C,,. The output of LPF (low pass filter), 
ip and iq, are the direct current value of i, andi, respectively. I, 


E P 4 
i,» 1,-are 3-phase fundamental currents respectively. Afterwards, 


the compensation current of harmonic current is: 
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Based on the aforementioned harmonic signals which are 
obtained according to the 1. method, activity of appropriate 
filters can be fulfilled. 

In the designed SAPF/PFC circuit, considering that the 
conventional source voltage provided by a ship electrical 
power station is distorted in evidence and even a little bit 
3-phase unbalanced in a degree, the 1 A method is applied to 
determine the compensation current. Furthermore, in Fig. 2b, 
the channel of current i is disconnected, then the signal 1, that 
goes through the LPF is 0, and the equation (8) becomes: 


ai 
ig AS a (11) 


1 cfp 


where only active current can be found in its output of i£, and 
ig, that is, the output of transposition of coordinate converter 
(CI) leaves fundamental active current only. The difference 
in value between the detected current (i, i,, i,) and the output 
of Cj, is the sum of fundamental reactive current ap ie 
i,,,) and harmonic current (i,,, i,,, ip). At that time, the final 
result of current calculation in Fig. 2b expresses the sum of 
harmonic current and fundamental reactive current. When 
this result acts as the instruction value of the compensation 
current, the compensation of fundamental reactive current and 
the compensation of harmonic current are able to be fulfilled 
simultaneously. Furthermore, since the method of 3-phase 
asymmetry control is adopted, the negative sequence current 
produced by the system itself counteracts the input negative 
sequence current, so that the simultaneous compensation of 
negative sequence current is realized [7]. 

The reactive and harmonic current detection algorithm 
adopted in the active PFC circuit in the Hybrid APF/PFC 
Circuit is the i-i method as shown in Fig. 2b, which is under 
consideration for its function of both power factor correction 
and harmonic current compensation, and the source voltage 
is distorted in ship electric networks. 

For the SAPF, in order to compensate harmonic voltage 
accurately and in real time, its harmonic voltage detection 
algorithm is similar to the harmonic current detection algorithm 
which is shown in Fig. 2a, where supplied voltage of v,, v,, V, 
are measured from the output of the power supplier, then v 
and v, are obtained through matrix transform of C,, and C 
and after the low pass filter, Vp and va, which are DC quantities 
of v, and v, corresponding to 3-phase fundamental voltages 
V e Vip Vep are obtained. 

After matrix transform of C,,and CZ, 3-phase fundamental 
voltages V p V,p V,, are detected. Afterwards, harmonic voltages 


Vane Ven Ven ae obtained, according to which the 3-phase instruction 
voltage of the power electronic converters, corresponding to —v.,, 
-V w “Vœ are ascertained and adopted to control the on-off state 
of the six switching units in the main circuit. 

For the shunt passive filter in the Hybrid APF/PFC/PPF 
Circuit, its parameters are designed according to the L-C 


Resonance Principle to suppress the 5 and 7" order current 
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harmonics, as well as higher order current harmonics, caused by 
the nonlinear load while it operates in rated conditions, since the 
eigenvalue of harmonics caused by the adopted nonlinear load 
is 5" and 7" order harmonics, 11" and 13" order harmonics, etc. 
The operational performance is not so satisfactory if the shunt 
passive filter operates alone, for a shunt passive filter does not 
employ any control algorithm and its parameters are fixed. 


SIMULATION EXPERIMENTS 


The simulated experimental platform of a marine electric 
network — including the marine synchronous generator, as 
well as the linear load and nonlinear load — is established. 
The parameters of the simulated power source are those of 
a synchronous generator for marine use which is made in 
Shanghai Electrical Machinery Plant of Shanghai Electric 
Group Co., Ltd. Its type is TFH-400/6 - Salient-pole 
synchronous generator. Its nominal data are: 400 kVA, 400 V 
line voltage, 50 Hz, 1000 rpm. Relevant reactance data are (pu): 
d-axis reactance (X) is 1.728, d-axis transient reactance (X) 
is 0.384, d-axis subtransient reactance (X,”) is 0.1793, q-axis 
reactance (X,) is 0.865, q-axis subtransient reactance (X, ”) 
is 0.181, leak reactance (X) is 0.113. The simulated current 
harmonic source is a 3 -phase nonlinear load which is a 3-phase 
full-controlled rectification bridge circuit followed by a resistor 
of 10 Q, and the trig angle of each of the 6 thyristors in the 
rectification bridge is 0°. The linear load is a 3-phase RL load 
with nominated power 20 kVA. In the shunt passive filter, the 
values of “L” and “C” in each phase are 2 mH and 200 microF 
respectively in the 5th single-tuned harmonic filter, 2 mH and 
100 microF respectively in the 7th single-tuned harmonic filter, 
and 0.37 mH and 175 microF respectively in the high pass filter, 
and “R” in each phase of the high pass filter is 5Q. For the PFC 
circuit and the SAPF, the capacitor in shunt connection in the 
DC side of the converter of both the PFC circuit and SAPF is 
2 microF, while each phase reactance that connects the PFC 
circuit and the network, as well as the SAPF and the network, 


is 5 mH. Furthermore, for the design of the SAPF, its low-pass- 
filtering part adopts the second filter module with its cut-off 
frequency 20 kHz. 

The registered voltage waveform and current waveform in 
the sub-switchboard of the nonlinear load are shown in Fig. 3. 
In Fig. 3, THD,, of the voltage supplied to the nonlinear load is 
4.94%, THD, of the nonlinear load harmonic current is 30.55%, 
and the network total power factor is 0.72. 

In the simulation experiment, a comparison between the 
operation of a combined system of shunt passive and series 
active power filter and the operation of the Hybrid APF/PFC/ 
PPF Circuit is carried out. In the first step of the experiment, 
only a combined system of shunt passive and series active 
power filter operates in the simulated marine electric network. 
The experimental result of the network voltage and current 
waveforms and their spectrum are shown in Fig. 4. In the second 
step of the experiment, the Hybrid APF/PFC/PPF Circuit is 
put into operation in the simulated marine electric network. 
The experimental result of the network voltage and current 
waveforms and their spectrum are shown in Fig. 5. In Fig. 4, 
THD,, of the voltage supplied to the nonlinear load is 3.89%, 
THD, of the nonlinear load harmonic current is 12.58%, and 
the network total power factor is 0.90. In Fig. 5, THD, of the 
voltage supplied to the nonlinear load is 3.13%, THD, of the 
nonlinear load harmonic current is 6.19%, and the network 
total power factor is 0.99. 

In the simulation experiment, a comparison between the 
operation of a combined system of shunt passive and series 
active power filter and the operation of the Hybrid APF/PFC/ 
PPF Circuit is carried out. In the first step of the experiment, 
only a combined system of shunt passive and series active 
power filter operates in the simulated marine electric network. 
The experimental result of the network voltage and current 
waveforms and their spectrum are shown in Fig. 4. In the second 
step of the experiment, the Hybrid APF/PFC/PPF Circuit is 
put into operation in the simulated marine electric network. 
The experimental result of the network voltage and current 
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Fig. 3. Voltage and current waveforms and spectrum of harmonic source, 
a) waveforms of network voltage and nonlinear load current, b) spectrum analysis. 
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Fig. 4. Voltage and current waveforms and spectrum obtained after the operation of a combined system of shunt passive and series active power filter, 
a) waveforms of network voltage and current, b) spectrum analysis 
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Fig. 5. Voltage and current waveforms and spectrum obtained after operation of the Hybrid APF/PFC/PPF Circuit, 
a) waveforms of network voltage and current, b) spectrum analysis 


waveforms and their spectrum are shown in Fig. 5. In Fig. 4, The effectiveness ofthe filtering circuits under consideration 
THD, of the voltage supplied to the nonlinear load is 3.89%, was evaluated not only based on their filtering properties, expressed 
THD, of the nonlinear load harmonic current is 12.58%, and by related THD values, but also taking into consideration such 
the network total power factor is 0.90. In Fig. 5, THD of the parameters as max (I,/I,), active power factor A, estimated cost of 
voltage supplied to the nonlinear load is 3.13%, THD, of the realization and the influence of Af = + 5% f (rated frequency) 
nonlinear load harmonic current is 6.19%, and the network changes on the filtering properties. It was verified that the 
total power factor is 0.99. efficiency of a passive parallel filter (PPF) alone is compromised 
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when the frequency is deviating. So to improve this situation, 
a hybrid filter construction should be applied, but at least two 
criteria must be taken into account at the same time: harmonic 
mitigation and sensitivity to frequency deviation. The solution 
based on combined SAPF/PFC/PPF configuration displays the 
best ability for power quality improvement; however, it is quite 
an expensive solution and the most susceptible to the marine 
environment itself. It requires careful design and manufacture 
of hardware parts in order to be resistant to marine conditions, 
i.e. vibration, humidity, high temperature. 


CONCLUSION 


According to the comparison of the experimental results 
of network voltage and current waveforms and their spectrum 
in the environment of a simulated marine electric network, it 
can be concluded that the combined system of shunt passive 
and series active power filter is able to improve the network 
harmonics and power factor, but it is not adequately effective 
in suppressing current harmonics and improving the power 
factor. In contrast, application of the Hybrid APF/PFC/PPF 
Circuit gives satisfactory result in both voltage and current 
harmonics suppression and power factor correction. That is, 
due to the importance and complicated nature of the problem 
of marine electrical power quality, application of the designed 
Hybrid APF/PFC/PPF Circuit to a marine electric network 
to improve the problem of harmonics and low power factor 
is a reasonable choice since it shows considerable merits 
in harmonic suppression and power factor correction in the 
environment of a marine electric network. 
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